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SYMBOLS 


State (Ns,Ns) or output (No,No) weighting matrix 
control (Nc,Nc) weighting matrix 

control gain matrix (Nc,Ns) 

ccentrol (No,Nc) or noise (No,Ng) feedforward 
matrix 

open-loop dynamics matrix (Ns, Ns) 

control distribution Batierasins, We) 

state disturbance distribution matrix (Ns,Ng) 
measurement scaling matrix (No,Ns) 

estimator gain matrix (Ns,No) 

humber of ccntrols 

humber of process noise sources 

number of states 

number of okfservations or meaSurements 

white process noise covariance matrix (Ng, Ng) 
white meas. noise covariance matrix (No,No) 
steady-state covariance matrix of control (Nc,Nc) 
control vector (Nc,1) 

COBtrO] An Pute ie aD 

state vector (Ns, 1) 

State vector derivative (Ns,1) 

estimate of state vector (Ns,1) 

derivative of estimate of state vector (Ns,1) 
State reconstruction error (Ns, 1) 


output/measurement vector (No, 1) 
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I. INTRODUCTION 
The purpose of this thesis is to describe and denon- 
Strate the modification and additions to the existing 
FORTRAN program (OPTSYSX) Which 1S used in the study, 
design, and application of Optimal Systems Control theory. 

The CPTimal SYStems control program (OPTSYS) was Origi- 
NMated in 1971 by Hall [Ref. 1] to support his research in 
EOrary-wilng aircraft control systems. The most recent 
Frogram godifications were made by Walker [Ref. 2] and Liu 
[Ref. 3] of Stanford University and are designated OPTSYS 4 
and CFISYS 5, respectively. The OPTSYS modifications made 
Ey Hoden [Ref. 4] were primarily devoted to creating a 
user-friendly interactive version (OPTSYSX) of the CPTSYS 4 
Soopers 

The goal of this thesis work was to develop a progran 
set which will operate in an interactive mode and plot the 
time response of a State Variable Control System which has 
been developed and/cr analyzed using the OPTSYSX Progran. 
Minor modifications tc the OPTSYSX Program were necessary to 
allow the user to build a data fiie of the matrizses reguired 
for the time response calculations plus additional matrices 
which cculd be used again upon reentry to the OPtTSYSX 
Program without the laborious and time-consuming task of 
reentering each data element in all of the required 
matrices. 

It is assumed that the reader/user is: familiar with the 
Easic concerts of Control Theory and Optimal Systems Design. 
The symbcl/naming conventions of Bryson [Ref. 5] are used in 
the projram operation discussion and in the exanples of the 
problems solved uSing this’ systen. A glossary of the 
symbols and abbreviations used ime this discussion is 


provided on page 8. 
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An explanation of the OPTSYSX capabilities and a program 
set overview are presented first. 

This work concludes with examples of variois types of 
problems demonstrated in the interactive mode, including a 
copy of each terminal session with the final results. A set 
of complete program listings are included in Appendices A, 
E, € and pe 
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A. GENERAL BACKGROUND 


OFTSYSX 1S an interactive, double-precision FORTRAN 
program employing modern control theory analysis techniques. 
Its extensive capabilities include the svnthesis and anal- 
ysis of filters and regulators as well as eigensystem anal- 
ysis, modal distribution, transfer function analysis and 
power spectral density computations. The modifications to 
the OPTSYS Program introduced by this thesis work have not 


affected any of the program's original capabilities. 


B. PROGRAM OVERVIEW 


CETSYSX is an extremely large and complex program with 
over 2800 lines of code. In order to use this program in 
its small (set up for 32X32 matrices) version, the user is 
reguired to extend the IBM 370/3033 virtual machine (VM) 
memory capacity beyond 720 kilobytes which is th= default VM 
memory size. A significant increase in the size of the 
OPTSYSK program would make the program too large to operate 
On a cne megabyte VM, the largest virtual memory available 
on a user's virtual machine. The high resolution plotting 
software is limited to Singie precision yariables. 
Therefore double precision library routines cannot be called 
EEOn the pPLetting program. For these reasons the task of 
obtaining the time response of a system was divided between 
Piece se rOdgranS worl >s YSki,sOPTCALGCsand OPTPLOT. An Executive 
program (OPTSYS EXEC) was written to make the interfacing of 
the three programs transparent to the user. 

Minor modifications were made to the OPTSYSX program 


including the addition of three subroutines to handle the 
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input and output of Matrix data to anderrom a data firilemon 
the user's disk. The OPTCALC program performs the double 
precision pumerical integration of the system df egquations 
over time and creates another disk data file of the state 
variakle variation with time. OPTPLOT takes this time 
response data and presents it ina graphical format on the 


TEK 618 craphical display or aS GER CA Tees penmeot. 


Co OL vote Ase 


The OPTSYS EXEC is written in the EXEC 2 lanjuage. This 
language allows the EXEC to issue almost any command that 
can ke entered in the direct mode at the terminal. Therefore 
an EXEC is the ideal controller for the "black box" type o£ 
system where the user 1S not aware of what is actually 
taking place within the program(s). The OPTSYS EXEC was 
written to complete all of the reguired interfacing between 
the three programs (OPTSYSX, OPTCALC and OPTPLOT), without 
the direct guidance or control of the user. By answering 
guestions presented on the terminal screen, the user deter- 
mines the logic flow through the EXEC while the EXEC estab- 
lishes the appropriate FILEDEFs and loads the programs 
required by the user's desires. 


D. OFTSYSX MODIFICATIONS 


Three subroutines (RDMATF, RDMAT, and WRTMAT) were added 
to the CPTSYSX Program for data  f£ile” read error “check, 
matrix input from a data file and matrix output to a data 
file, respectively. , These three subroutines provide the 
user with the opportunity to save the [F], [G], [ Hoe 
[GAMMA], ,;A]j and [|B] matrices for use 1n a subseyuent run of 
the OFTSYSX Progran- The WRIMAT subroutine also saves the 
[C] and [K]} matrices for system time response calculation 
and plotting by the OPTCALC and OPTPLOT programs. 


13 


1. RDMATF Subroutine 


The RDMATF subroutine is used to check for the exis- 
tence of a previously generated file containing matrix data. 
Seven flags may be set by this subroutine. Six of these 
flags correlate with the six matrices that the user may save 
ror reuse later in the OPTSYSX Progran. ine remaining £lag 
(IROMAT) must be set to enable the RDMAT subroutine to read 
matrix information from the data file. A READ statement of 
the forn 


READ (9,111,ERR=2 22, FND=333) A,3 


(where "11171", "222" and "333" designate line numosers for the 
FORMAT statement ana DG anew On ERROR ofr Draten On 
Pugeot~-“SITLE routines, respectively) is used for the data 
file check. The nenexistence of the File or premature 
END-OF-FILE are detected by the ERR and END cshecks which 
cause a Eranch to a routine that sets the IRDMAT flag to "0" 
apa becurns to the calling progran. 

If no error is detected during the initial read 
attempt, the variable Bis checked for the sentinel "1", 
This seccnd check is to help ensure that the file is actu- 
ally a file which contains valid matrix elements. The user 
is then presented a message which asks if he/she wants to 
use the matrices which are available. The user may respond 


with on cr three answers: 


1. Use all of the matrices. 
2. Use selected matrices. 


3. Use none of the matrices. 


If the answer is "i" or "2", the subroutine reads the matrix 


dimensions (Ns, Nc, No andNg) from the data file and 
changes the IRDMAT flag to 1 to key the RDMAT subroutine to 
read the matrix elements from the data file. If the anSwer 


is "1" all of the matrix-save flags are set to "1", 
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If the answer is "2", the Wsemmaceucamien Ehe opportt 
nity to select individual matrices for reuse while rejecting 
other matrix information. This is accomplished by setting 
individual matrix-save flags to "1" if the matrix is to be 
saved and "0" if new matrix data will be input from the 
terminal. If the answer is "3" (Use none of the matrices) 
the IRDMAT flag 1s set to 0O and the subroutine returns to 
the main program. When all actions have been conpleted, the 
flag infcrmation and the matrix dimensions are passed to the 


main program for later use. 


2- RDMAT Subroutine 


The RDMAT subroutine is used to read all of the 
matrix information in the data file and transfer the infor- 
Mation to the appropriate variables. AS previously 
discussed, The acticns of this subroutine depend on the 
status of the IRDMAT flag. If this flag had been set to 
"mo", no read operations are attempted and program flow imme- 
diately returns to the calling progran. 

When the IRIMAT flag is set to "1", the kRDMAT 
subroutine reads the matrix dimensions from the data file, 
and uses these dimensions to transfer the matrix information 
from the file to the appropriate variables. The file matrix 
dimensions are used for the read operations and are not fed 
back to the calling program, since the dimensions of some of 
the matrices which are not being reused may have changed 
from the the previous run. Similarly, using the cugpene 
matrix dimensions in the RDMAT subroutine would cause data 
read-in frotlems due to the changing number of elements in 


each matrix as the matrix dimensions vary. 
3. WRIMAT Subroutine 


The WRITMAT sukroutine 1s used to write a data file 


of the data file flags, the matrix dimensions ind selected 


V5 


matrices. When the user has completed the analysis/design of 
the system of interest, the WRTMAT subroutine asks the user 
if heyshe wants to calculate the time response of the system 
which the user just designed. If the user answers YES, the 
WRTMAT subroutine generates a data file of appropriate 
matrix information and halts execution of the OP‘SYSX 
program. Control then reverts to the OPTSYS EXEC. If the 
user answers NO, the WRIMAT subroutine returns control to 
the main frogram and normal OPTSYSX program oferation 
continues. 

The informaticn written to the data file consists of 
2"1"s (which are used as a sentinel or flag by the RDMATF 
Subroutine (aS previously explained) and in a similar manner 
ry the OPTSYS EXEC), followed by the matrix dimensions (Ns, 
lemceana Ng) and then by the [F], [GJ], | H], [GAMMA], ([(C], 
[Xk], {A4] and [B] matrices. These matrix elements are written 
Moeere OPIMAT PATA file using a 4D20.13 format is a compro- 
mise kLetween the maximum feasible accuracy of data exchange 
rFetween the double-frecision programs and the use of a 


moderate amount of the user's disk space. 


E. OFTCALC PROGRAM 





The OPTCALC program is a FORTRAN interactive 
doublée-precision system integration routine. i hese DEeG Lan 
uses the International Mathematical & Statistical Library 
(IMSL) subroutine DGEAR to perform the numerical integration 
of the system under analysis. The stiff system mode of DGEAR 
is used in order tc provide the capability to do time 
response calculations of the X-29A longitudinal axis back-up 


mode system which is an 98 X 98 stiff systen. 
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2. System Eguation Representation 


The OPTCALC program uses the state variable format such as 
xdot = [F]*x + (3 J*¥uc (2.1) 


to define the system. In this system the [F] matrix is the 
open-loop dynamics matrix (system or plant) and the [G] 
matrix is the control matrix. The variaole assignments are 
x aS the state vector and uc as the control input vector. It 
follows that xdot is the time derivative of x. 

Various forms of eguation 2.1 are used for all the 
time response calculations. The [F] matrix is modified to 
[F+S*C]i for closed-loop (regulator only) system calcula- 


tions @S in ejuation= 2.2 aac 


Xdot = .F+G=e Ve xeet | Gieue (252) 


u = (CC iXX.t ue (253) 


For this closed loop system, the [Cj matrix is the centrol 
gain cr regulator Gain matrix and uis the total input 
VECCOF., 

The combined filter and regulator systems can be 
represented by equations 2.2, 2.4, 2.5, 2.6 and 2.7. The 
[HE] matrix is the measurement scaling (observer output) 
matrix and the [K] matrix is the estimator or Kalman filter 
cbserver gain matrix. The variables xe, xedot and x are the 
state estimate vector, the derivative of the state estimate 


vector and the state reconstruction error, respectively. 


ee ee eG es ee 


1The OPISYSX sign Convention tor tice Eee is the 
negative of the “standard normally _used GChtiromss 
Therefore rF4Gee | has the correct sign for OPTSYs x matrices. 
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xdot = [F+G4C ]*x + [G]*uc rez) 


z= [ H }*x (2.4) 


xedot = [F]¥xe + [G]*u + [K]*(z - [H]*xe) (Ze) 


S 
I 


[ C }*¥xe (250) 


xX =x - xe (257) 


Pegaraons 222, 2.4, 2.5, and 2.6 Can be combized into the 


augmented matrix form of equation 2.8. 


xdot | | F G*C i | 
= | | * | 
xedot | K*H F+#+G*C-K*H | | 


Sas Ge ie 


Equation 2.8 1S an augmented eguation in whicn the Ns dimen- 
Sion has been doubled andthe state and state estimate 
vectors have been combined into one vector of 2*Ns length. 
The [G] watrix is also augmented by repeating the first Ns 
rowS again beginning at row Ns+1, making the new [G] matrix 
dimensions (2*Ns,Nc). | 

In a Similar manner, a filter only system can repre- 
sented by the same eguations with the [C] matrix set to 0. 


These equations are: 
moot f=. (Rie xX tb (CG | *uc (21) 
ZH) ex (234) 
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xedot = [PP }*xe +9 6 Jeu. hae ee x2 ) (2.5) 
% =x - xe (2a 


The comEination of equations 2.1, "2Zaateeeleoeand 2.7 invowan 
augmented [F]} matrix is Similar to eguation 2.8 but with 
fewer terms and the upper right guadrant eyual to zero. 
This filter-only augmented system equation is shown below. 


a 0 x | (2Ze55 


«| | + 
| K*H F-K*H | [xe] 


xdot 
|] xedot 


*uUC 
G 











a. System Selection 


When the CPTCALC program is run, the |Fl, |G] 
ard |C] matrices (and |H] and |K] matrices, if available) 
are presented on the terminal as a check and as a reminder 
of the characteristics of the system that has been passed 
from OPTSYSX to OPTCALC. The user is reguested to select the 


type cf system response to be calculated. 


4= Wh) 


Selection of 2, 3 or 4 forms the appropriate system matrix 
ecuationse 2.2, 25 9 crease respectively and doubles the 


length of the |G| matrax, if required. 
k. Defining Calculation ms tseandeinput s 


After the user determines the type of systen 


under study, the OPTCALC program prompts for the integration 
start and stop times and the number of data points desired. 


Lg 


The user has some centrol over the tradeoff between curve 
fidelity and computer time used by varying the number of 
data points calculated. Computer time use is normally a 
factor only on very large systems. If less thans200 points 
are calculated, the OPTPLOT program uses a curve smoothing 
function which may cause minor inaccuracies in the plotted 
curve but avoids the sharp peaks and irregular appearance 
generated by plotting straight lines between an insufficient 
number of data points. When 200 or more points are calcu- 
lated no smoothing is done. The points are connested by very 
short straight lines which has the appearance of a smooth 
Curve. 

step or ramp functions are available as control 
inputs. Only one type of function can be used for. each 
control, but the magnitude and start and stop times can be 
set as desired. 

The state and state eStimate initial conditions 
can ke set to any valve by the user. The control inputs and 
non-zero initial conditions can be used simultaneously, if 
desired. PFefore the time response calculations begin, the 
user iS given the opfortunity to make changes in any area of 
the system integraticn initial conditions that have been 


previcusly selected. 
3. system Time Response 


Eguation 2.1 is evaluated directly in the open loop 
system response calculations. The FCN subroutine was 
written to evaluate the system of eguations for the DGEAF 
IMSL suproutine. Each time the FCN subroutine is called by 
CGEAR, it updates the control inputs (uc) and then evaluates 
each state derivative by summing all the terms across that 
row of the [F]*x and [G]*uc matrices. The same FCN sukrou- 
tine is used for all system integrations. As >xplained in 
the previous section, the [fF] matrix is replaced by the 


Paso Gio Matrix for closed loop System response problems. 
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a. Systems With Filters 


Augmented equations 2.9 and 2.8 are 1sed for the 
time response evaluation of systems with filters only and 
systems with filters plus regulators, respectively. The 
augmented matrix is developed as a dummy matrix and is then 
inserted as the [FJ] matrix with the dimensions doubled 
(2*NS). The ~G] matrix is also augmented by repeating the 
first NS rowS again beginning at row Ns#t1, maxing the new 
{G] matrix dimensions (2*Ns,Nc). 

The augmented system can be evaluated by simply 
doubling the old system row and column dimension (Ns) and 
calling the DGEAR integration subroutine. Using this 
method, the existing FCN subroutine reguires no changes to 


evaluate the augmented systen. 
be OPTCALG= Cpe oUE 


The OPTCAIC program uses FILEDEF 8 for the data 
file output as well as FILEDEF 5 to read and write to the 
terminal. The output data file contains the following 
discrete information: the matrix dimensions Ns and Nc, the 
augmented matrix dimensions, the number of data points 
calculated and a flag to indicate that an augmented matrix 
waS calculated. The [C] matrix iS passed to permit the 
calculation of u the total control input to the system. The 
final portion of the data file is individual data foints of 
time, external control input (uc) and each state (x) and 
State estimate (xe). This data file provides all of the 
data required by the OPTPLOT program to make a smooth graph- 


ical response curve. 


Fo OFIPLOT PROGRAS 


The CPIFLOT program 1S a FORTRAN interactive plotting 


program uSing the Display Integrated Software System and 


Za 


Plotting Language (DISSPLA) Ey Integrated Software Systems 
Corporation. This program provides the user a high resolu- 
tion graphical display of the system's time response ani if 


desired will provide a VERSATEC pen-plot of the same graph. 


1. General Operation 


Flotting data is received from the projram OPTCALC 
via a data file on FILEDEF 8. The types of data provided in 
the file are discussed in the previous section. The program 


presents the user with a series of guestions to determine: 


The number of curves to plot. 

select the type cf variable for each curve. 

Select the variable subscript for each curv:;. 
select the number of headings and contents of each. 


4> CON) od 
@ @ @ @ 


The progran then plots the selected variables and frovides 
the user a graphical display on the Tektronics 618 (TEK 618) 
Manso. aly 

The following Main Menu is then presented which 


provides the user with the major decision points of the 


program. 
leeeseGuN NEW GRAPH CF OTHER CONTROLS ,STATES, OR ESTIMATES. 
2. REFLOT PREVIOUSLY SAVED GRAPH DATA. 
Semeur) THE CURRENT GRAPH. 
eee teen evi ORD GRAPH ON THE TEK6E 18. : 
ae oan AO, O Reser: MElAF TES OF ITHE CURVES 
PRev@OuSLY SAVED. 


The purpose of each of the selections is  self- 
explanatory, however the methods of their use may not be. Ii 
number 3 ais selected, the user is then presented the 
following Edit Menu of items to make additions, deleticns or 
corrections to the curves that are plotted on the TEK 618 


=Green: 


ae 


1. CHANGE VARIABLES OR ADD A CURVE ON THE CURRENT PLOT. 
2. DELETE CURVE FRCM CURRENT PLOT. 

3. EDIT CURVE TITLE (S). 

WY. EDIT PAGE HEADING(S). 

5. CHANGE THE Y-AXIS SCALE. 

6. CHANGE THE TIME AXIS SCALE. 

7. CHANGE THE PLOT SIZE. 

8. CHANGE THE LETTERING HEIGHT. 

9. CHANGE POSITION OF THE LEGEND. 

10. EDITING COMPLETE. 


This extensive list of modification capabiliti2s provides 
the user with the tocl to make almost any imaginable altera- 
tion to an existing eelor. 

Since the OPTPLOT program receives the time response 
data from the OPTCALC program, item 6 of the Edit Menu 
cannot be used to expand the time scale beyond the time span 
previously calculated. Therefore the time axis change 
feature can be used cnly to select a subset of the original 
daytyare 


2. WERSATEC Pen-rlots 


The VERSATEC pen-plots are provided through the 
CDISSPOP [LOrtion, “ef DISS een. In order to use the DISSPOP 
feature, a device independent plot file called a metafile 
must Fe generated. 

To view the graphical time resporse plot on the TEK 
618 terminal, the TEK618 cption of DISSPLA must be called 
within the OPTPLOT flotting program. The graphical image 
data is then sent to the TEK 618 display screen. 

A metafile is created when the COMPRS option of 
DISSFLA is called by the OPTPLOT plotting program. When the 
Flotting pregram is executed with the COMPRS option, the 
Graphical image data is sent to a metafile on the user's 
disk. The TEK618 crtion and COMPRS o;tion are mutually 
exclusive (only one can be active at a time), therefore 
graphical data cannot go to both the terminal screen and the 


metafile, concurrently. As a further complication, the 
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TEK61€ and COMPRS oftions cannot be used in alternating 
pattern, first to originate and edit each graph and then to 
add this graph to a metafile possibly containing several 
cther ograrhs. Therefore if more than one pen-plot per 
terminal session is desired, some type of capability must be 
provided to save the information required to reproduce a 
given graph. 

Khen the user attempts to leave the carrent plot 
(le. selecting items 1, 2 or 5 of the Main Menu) the progran 
asks the user to save the current graph data for later use 
in generating a retafile. This feature provides the cafa- 
bility to save any desired graph data in order to later make 
a metafile and obtain a pen-plot. When the COMPRS oftion is 
used (by selecting item 5 of the Main Menu), any number of 
graphs may be addedto the metafile up to the limit of 
available user disk space (provided graph data has’ been 
previously saved). After exiting the OPTPLOT program, the 
CPTSYS EXEC asks the user if he/she wants a hard copy of the 
metafile that had been generated during the session. If the 
user answers YES, the OPTSYS EXEC cails the DISSPOP EXEC 
Wee the VRogkc op ticn. 

When the user exits the DISSPOP EXEC, the CPISYS 


EXEC gives the user the options to: 


N 
Ik. 


QIN 
e 0 e@ 


MewmoOgmenconeto TUN OrloYoX vagain allows the user to use all 
or part of the matrices that had been saved in tne data file 
without manually reentering each element. The OPTCALC 
option could be exercised if the user wants to use the same 
system matrices again, but chanje the contrdl input oF 


initial conditions or change the type of system (open, 
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closed, filter only cr filter plus regulator) that was eval- 


uated on the previous run. 
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III. S¥STEM USE AND EXAMPLES 
This chapter contains several basic examples of the four 
types of problems which may be solved using OPTSY SX, OPTCALC 
and OPTPLOT under control of the OPTSYS EXEC. Included with 


these examples are copies of each recorded terminal session. 


A. OQOPEN-LOOP SYSTEM TIME RESPONSE 


The following ofen-loop system example was taken from 
ibeiemo, pp 5.3 = 5.7 1. 

The full terminal session is recorded below, with user 
input at the left margin in lower case letters oor numbers 


below each "?", 


record on 
BEGIN RECORDING OF TERMINAL SESSION 
Peet —0.01/0.02 19:58:26 
Opes ys 
THE OPTSYS EXEC CONTROLS A TRIO OF PROGRAMS: 
1. OPTSYSX FORTRAN (SYSTEM ANALYSTS) 
2. OPTCALC FORTRAN CALCULATE TIME ESBONSE) 
3. OFTPLOT FORTRAN {(DISSPLA PLOTTING ROUTINE) 
EACH PROGRAM PASSES INFORMATION TO THE NEXT PROGRAM 
THROUGH A DATA FILE WRITTEN TO THE USERS DISK. IN THIS 
CASE, THESE FILES ARE "OPTMAT DATA" AND “OPTPLOT DATA". 
THE SIZE OF THESE FILES VARY WITH THE SYSTEM ORDER, AND 
CAN USE ABOUT 20% OF TEE USERS DISK SPACE. THEREFORE 
ENSURE THAT SUFFICIENT DISK SPACE IS AVAILABLE. 


Seer ern wi, AN Olu ERE NITRY =2O CONTINUE - 


He ek oe eK a ok ea ko a ak ke ok ik a oe ae ok ok ao a ok ok ok aie ake ik ake ote ake aie ak ae ie ic ik i ai oie ik i oie ok kK akc ok Kok 


YOU HAVE A DATA FILE NAMED 'OPTMAT DATA" ON YOUR A 
Deoweelnun lenis FREVTOU SLY GENERATED BY THE OPTSYS 
FROGRAM AND CCNTAINS THE F, G, H, GAMMA, A AND B 
MATRICES FROM THAT RUN. 


IF YCU WOULD LIKE TO USE THESE SAME 

THIS RUN, THE OPTSYS PROGRAM WILL 

DESIRED DATA AT THE APPROPRIATE TIME, 
IF YOU TYPE (Y) ES. 


ANY OTHER INPUT WILL RESULT IN THAT FILE BEING ERASED! 


MATRICES 
READ IN 
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Ae RK I KR ea kK ae oo oe ae a ee ee eK a ko oe ok a ak ok ee eK a ok ok ok eK Kk ok 
Yy 


He HK a OK a ae ea eae oe ee ae oe a a a Ke a Ke Ke ee ok OK ie ke a ko ee ok ok ok ok 


DO YOU WANT THE NUMERICAL OUTPUT FROM OPTSYSX TIO GO 
TO YOUR TERMINAL 74 683} ORO wee ok) PF lire 


eA I A a ee ee a ee ee a Ke eR a oe a ee ek i a Ok ae oe oe oe eo ak a Kok a oe ie oe ak a ok ok 


Ss 


koko ko ok igkok kk kok kok ok kokkok dokokok kak kok kak kook ak aka ak ake aka akaak ak ake ok ak ok ak ak ak ake ak ak ak akc ake ak 


OUTPUT WIIL COME TO YOUR TERMINAL SCREEN. 
FR kk kok kk ok kk kk kkk kk kk kak kkk kkk kk kok kk kk kk ok ok de 


LOR DNC NG ESy cme 
EXECUZTICN BEGINS... 
CPTSYSX IS A COMPLETELY INTE Y 
PROGRAM. IT WLIL SOLVE NUMEROUS CONTROL PRO 
FOLLOWING TYPES OF SYSTEMS CONTROL EQUATION 
XDOT = {F}*X + {G}*U + {GAM}* (W+WO 
MEASUREMENT EQUATION-- 
{H}*X + {D}*W + V 
REGULATOR PERFORMANCE INDEX-- 
J = 1/2 * INTEGRAL (Y¥Y *{A}*Y + U *{B} *U) DT 
STATE FEEDBACK GAIN DEFINITION-- 
U = = {C} ¥*xX 
DO YOU WISH TO CONTINUE? TYPZ "YES" CR "NON, 
== DATASENERY —— 
IS SPECIFICALLY DESIG 
NTERACTIVELY, SEVERAL 
ABLE TO USERS: 


Urt "GN AND "GAMMA" MATRICES 
D FROM SEPARATE DATA FILES. 


wri wGM) AND "GAMMA" MATRICES MAY BE 
DEF{NED WITHIN SUBROUTINE "SETUP", 


CASE, THE USER SHOULD OBTAIN A COPY 
OGRAM LISTING AND EXAMINE 
LES CONTAINED ING SK Vise 


TO CONT INGE Te ee eae UO Ne 


RACTIVE OPIITMAL = ONTROL 


N SI EUS 16 
BLEMS ON THE 

Ak Soe 

) 


Nm 
I) 


teH 


TG mii <ht ei HO ob 
ro oot 


DO YOU WISH TO INPUT GHEE “FG Vee eco 
MATRICES FECHE SUB ROUTINES EDU a eee 
METHOD DES@CREEED VON id Eee REV EO ee erreere 


TYP ey SS eras NiO 


a, 


a= *4) 


Seinieeer LoYsk OPTIONS: 
OO ee ere en aero St LH OUT 
Orie SiGENSLSLEM CALCULATIONS. 
SEITEN 2a es OloleNmatneyots WiTH OPEN=LOOP 
ELGEN SYSitaM, CALCULATIONS, 
OFIFON Sse OPEN] LOGEMEIGENSYSTEN FOUND 
AND PROGRAM TERMINATES. 
(ee hie TRY POLLOWS IMMEDIATELY. } 
OFTLON 4 == MODAL Sento AON MATRICES COMPUTED 
Sree eg Daeetane abiggegro® S¥ATAESHS 
Sols C lean er LON: 172,33, OR “4. 
DO YOU DESIRE RMS VALUES OF STATE AND CONTROL? 
ee es eee NOY, 
Pee i- LLCO? LRANSTER FUNCTION OPTIONS: 
OPTION 1 -~ NO OPEN-LOCP TRANSFER FUNCTIONS COMPUTED. 
CPOE aaa CLoS,ekEoIDUES, AND ZEROS COMPUTED. 
Cel tOv see eOlulL! sOLES AND ZEROSMeONPULEDE 
Celi Oem caecCNLY POLES AND RESTDUES COMPUTED. 
SELFCT AN OPTION: 1, 2, 3, OR 4. 
NOISE TRANSFER FUNCTION OFTIONS: 
OPTION 1 --= NO NOISE TRANSFER FUNCTIONS COMPUTED. 
Cr Owe -—BECLEo, —ReosIDUES, AND ZEROS COMPUTED. 
Co ON (3a eee LES AND ZEROS COMPUTED. 
COP tis Se NLY wero Lio saNven ES DUBS COMPUTED. 
SELECT AN OFTION: 1, 2, 3, OR 4. 
COUPENSATOR TRANSFER FUNCTION OPTIONS: 
OPTION 1 -- NO COMP. TRANSFER FUNCTIONS ZS OMPUTED. 
Cee ie 2e-—weClLEo,enESIOUta, AND ZEROS COMPUTED. 
CPTION 3 -- CNLY POLES AND ZEROS COMPOTED. 
CPTTON 4 -- CNLY POLES AND RESIDULS COMPUTED. 


{NOTE: A COMPENSATOR TRANSFER FUNCTION CAN BE 
CCMPUTED ONLY IF BOTH A REGULATOR 
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Cue) 


AND FILTER ARE SYNZEES@ ZED 
AND/OR INPUT.} 


SELECT AN OFTION: 1, 2, 3, OR Sie 


WELL A FEED=FORWARD DISTRIBUTION Bia Tr x 
{"O' = MATRIX} 2h TNE Ue 


LYE ER SyEsS On “Nee 


THIS OPTION, DETERMINES. THETCRITERIA FORD DECIDING Sh teat. 
MARKCV FARAMETER IS ZERO-THE “HARKOV PARAMETER INDICATES 
THE CRDER OF THE NUMERATOR POLYNGETAURCr BACH SIR ANS 
FUNCTION. 
ALL "N" ZEROS CF THIS POLYNOMIAL ARLE PRINTED Cul AN D 
THIS TEST TELLS HOW MANY POPE A reese rer ol Ala On 
LESS THAN 10.07* {-IE} JS CONSdD Eee Zn aG. 
IHE DEFAULT VALUE OF THIS @eARAE Cee UE} Sac. 
IN CTHER WORLDS, IE — S30 ar 
if YOU DESIRE A DIFFERENT MNARRC Chl 2a, 
TYFE THE INTEGER VALUE. 
IF YCU DESIRE THE SEPAD rR eA Ue ee Ol 7 ERO} 
POWER SPECTRALSDENSITY Si PSpimeeEt ION @ 7; 
OPTICN te=- COMPUTE THE PS DPOReiEE OUTPUTS aD 72 he tae 
CONTROLS OF THE CONTROLLED SYSTEM WHEN FORCEDRE 
PROCESS AND NE ASURERINI NOISE. ae BOTH. 
REGULATOR AND A FILTER MUSI BE RESIDENT IN THE 
PROGRAM 10 UsSt Tits eGe acon.) 
OPTION 2 -- SAME AS OPTION APO ee aeOn ty eet fue 
RESIDUES OF EA CH Ris Erne Gc Trey 
USED INSIHE eS ieee yeu Viton. 
OPTICN 3 -= "NOT DES IEeae 
SELECT AN- OPTION: pec gmok S. 


THE "FU, "GH mAW) ugAM"™, "AN AND "3B" MATRICES 
FROM YOUR. PREV [ous OPTSYS KUN WERE SAVED. 


THE FOLLOWING OPTIONS ARE AVAILABLE: 
1 o.US2 Ales OF THE SAME MATRICES ASAIN. 
2. USE SELES TEDEMAIRICES AGAIN 
3: INPUT ALL NEW MATRICES. 


ENTER 1,278 On see 


NOTE: EACH SAVED MATRIX WILL BE REDISPLAYED Al 
THE PROPER INPUT (SEQUENCE INTERV 

AND YOU WIIL HAVE THE OPTION OF CHANGING 
INDIVIDUAL MATRIX ELEMENTS. 
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FIAG/PARAMETER SEITINGS FOR THIS RUN ARE AS FOLLCWS: 
ee ee SS Lie leer 2 eeiTr3 IFDFY ITE IDEBUG 


3 0 0 0 1 0 0 0 0 O 0 
ISET 1DSTAE IPSD IYU0 INORM IREG NS NC NOB NG 
0 0 0 0 0 0 4 O QO 
OP DPR@OF SYSTEM = 4 
NUMBER CF CONTROLS = 0 
NUMBER CF CBSERVATICNS = 0 
HentoreeOr FEROCESS NCISE SOURCES = 0 
ites YomeMee As Roker OP a ATR TMs. 
Ome 1.00000 0.0 Oo 
Oe =i 4 1500 =O et 0 ° 
wee 000 0 -1.43000 =0 0198.0 0.0 
00 0.0 1.00000 0.0 


DO YOU WISH TO CHANGE THe VALUE OF ANY MATRIX ELEMENT? 
TY Poe o OR NOM , 


cs MAC Re ha reve retateta se eleteves «0 «6 fe 


OPEN LOOP DYNAMI 

00L+01 Diag Hi--0 
505r+00 - 
30 


800D+01 


Owoooe 


ODdoo 
COOe 


QO. 
0 


Ger MOORE GENWNEUES: 25... .6..DE 7 (SI-F) .. 
0.0 1-6. 80767D-01: 1.22984D-01, 3.80349D-01: 


MCDAL MEASUREMENT SCALING MATRIX...H (BAR) *TI.. 
0.0 Ona0 O'. 0 0.0 
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DO YOU WISH TO OBTAIN A TIME RESPONSE 

OF THE SYSTEM YOU ARE EVALUATING? 
(Y OR N) 
NOTE: YOU MUST BE LOGGED ON AT A DUAL SCREEN 
(TBK 618) TERMINAL TO UTILIZE THIS MODE. 
THE a S¥st BM) G (CONTRO 1) by H (OBSERVABLES) GAM {NOISE), 
(cu C OST) AND B, (eo NTROL Cost) MATRICES WILL BE 

ED FOR RSENTRY HE MAIN OPTSYS PROSRAM. 


y 


AK fe eo oe te he ee ee eo 2 a ee eK oe KK he ee he eo he ee ek ee KK EK KK KO 


If YOU ARE EISSATISFIED Sian tae oe 
THUS FAR AND WOULD LIKE le Est re cus, 


=IYPE *Y* 10 52e— 


ANY OTHER INPUT 10 GeoNnTINGE 
KK Kk KK 


PP ee he ek mare re Beenie 5 Be ek Medes +k 4x oo 


IOADING CPTCALC. ... 
EXECUTION BEGINS... 


DURING THIS SECTION OF THE PROGRAM YOU WILL: 
- SELECT. THE TYPE CE Sis ie scels 2 etme Lon 
(CFEN LOOP, CLOSED LOOP, OR PULren REGULATOR! 
- PROVIDE START AND STOP TIME FOR PLOTTING CALCULATIONS 
=- SELECT THs Tigh Cl URE see ee ee ea OR a 
- PROVIDE SIART AND SICP tae ore RIVING FUNCTION (S) 
- FROVIDE DRIVING FUNCTION MAGNITUDE(S). 
CLEAR THE S@RER Nees CON PEUE 
THE F MATEUS 
0.0 1.00000 0.0 0.0 
0.0 -0.41500 =O 2000 10 0.0 
¢.80000 -1.43000 -0.01980 0.0 
0.0 0.0 1.00000 0.0 
THE G MATRIX 
0.0 
Gao 
0.0 
0.0 
THE VC NADER x 
0.0 0.0 0.0 0.0 


THE FOLLOWING PLCTTING OPTIONS ARE AVAILABLE IF THE 
REGULRED MATRICES WERE CARPCUDAT ED  fN (Oe ters 


1. OPEN LOOP TIME RESPONSE 
XDOT = {F}*X + {G}* 


2. CLOSED LOOP TIME RESPONSE 
XDOT = {F-G*C}*X + {G}*UC, U = {C}*X 


3. OPTINIZED FILTER CLOSED LOOP es) Siar eer. 


cm 


XDOT = {F}*X + {G}*UC 2 = {H} *X 
XHDOT = ch # XE Shy o + {K}*{Z - deka} 
4, FILTER + REGULATOR CLOSED LOOP SYSTEM RESPONSE. 
XDOT = {F+G*C}*X + {6} *0¢ Z = {H} *X 
XHDOT = {F}*XH #{G}U + {K{*(zZ - H*XH}, U = Cf XE 
SHUEC Tats Oe eomonaiie 


AT WHAT TIME DO YOU WANT TO START 
thi taeis RESECNSESeALCULAT TONS ? 


PPO ol Akh lee EN SECONDS. (NORFALLY 0.0) 


? 
0 
AT WHAT TIME DO YOU WANT TO STOP 
THE TIME RESPONSE CALCULATIONS? 
, INPUT STOE TIME IN SECONDS. 
Os 
THIS PROGRAM DIVIDES THE TIME INTERVAL YOU iAVE 
JUST SPECIFIED INTO UP TO 500 SMALL INTERVALS FOR 
THE INTEGRATION AND PLOTTING ROUTINES. IN ORDER 
TO SAVE COMPUTER TIME, THE NUMBER OF POINTS CAN BE 
CAN BE REDUCED WITH SOME LOSS IN CURVE FIDELITY. 
, HOW MANY PCINTS DO YOU WANT TO CALCULATE? 
500 


DOZS THE SYSTEM UTILIZE A DRIVING FUNCTION (CONTI RCL INEUT)? 
(Y)ES OR (N)O 

nH 

DOES THE SYSTEM START WITH ALL INITIAL CONDITIONS = 0.0 ? 
(Y) ES OR (N) 0? 


Nn 
5 Colt toe PNG APMCONDITION FOR X({ Ij 2? 
OD 2 
Wie eee om rte iN Men ewCONDITION FOR X{( 2) ? 
0 
: iA tees orn hienteCcONDETION FOR X( 32) 7 
0 

WHAT IS THE INITIAL CONDITION FOR X{ 4) ? 
0 

HAKE CHANGES IN THE FOLLOWING AREAS. 

(ccc NOT HEN UYPE OF SY¥sten. 10 PLOT 


or 


FILTER SORE ten: EGULATECR) 


CLCS El, 
SLART AND SIOP TIMES 


(OPEN, 


Ze 


FUNCTIONS 
INIT TAEWCON DET IONS 
CONTINUE 


DRIVING 


Be 


4. 


De 
SELECT A NUMEER BETWEEN 1 AND 5. 


RAM OPERATION. 


VIDED ONLY 


5 SPR 
PROG 
SLATES AND STATE ESTIMATES CAN BE PLOTPEre 


ALLVCCHNT ROL, 


X (3) 


X (2) 


X (1) 


U(1) 


Dee 


DOO OOOO OCOOO OOO OOOOGCe eee Kes =. OO-reee ee oe ere rec 
DDODDODOD ODI OOOO DID DODOO DODO VDOWDODOD GDODODOD OODOIODOMVODO DOGVOVCOO 
ree er He He eH HHH HHH IL EH HH HEHE HEHE CH CHC HEC HE Ih HEHE HOH HH HHH + + 
RAIA AAA AAR AQAQAA AIFARDQANARAAAARARARANRARARARAAQAAQAANRA 
NO OT ODO AAN TOOT YN FT STrVDOONMODVO ARNON AMrKnowadanMsm™MNOoowArmy) 
WOMOMNOWOM FNWMIDOWINANM™ OM MN MDDS OM ODM ONM ANT HDONONKMMEE MO? 
FANS SCOUMN FDOMNSCSHM FONMIN SWNOM™ NOWOWD NNO D QHonwWwONAOowr 
NM NEM OOM ATWO TFOMNNMYN gH SFOOW SFHHWMNAHNA AVOIO™ OODOMMN Done +zovrw 
MMM N OWN DID OW SFMN— Ot SH DODO STONOMMHYNOOMDOVMOODSPOAMNNNE KN C-Ore 
MHNOMMIINM TAKE OOM TK HODODMUMO HK MMNDOMMONNNSToOPMmormMandosalor zor 
DDT NA ASTNYI DOO WD AFNSCMH ANE DCCC KT CK KT TONVUATONNMN ts tT STAI MNHONNGN 
ee?eee f@ @ @ @ @ @ &@ @ e@ & ee &F eB we ef © © & © Ff © ef © eG e@ eg © 8G © &F © ee © F&F FF FS BC S&F FF SF ge 
BS Oe OO ae OO meena Lia ao 
;rfered’edeveeret &€ bg i | 

MOMNANCINNA QE eC CT OOK NNNNNNC CCCP Cr Cre cer OmNNC eo er Coon 
DODDOOCDOOSCSQDODOO ODOC VDOOODODOCOOCODOOOOOCOOOQOoGOOCOoO® 
(rrtrntrrrttveeretr’?tvdetttrtetvrt)ttetrret’trr)yrer’)tvprrrrt et!) -vpvererve)¢peva 
QAPRQAQQAQARQQAQAQQANAAQAQAANARAQAAQAANRARAAAQARAARAQAAARAAQAR 

NM YNDWT MEM ONS KH ODO NNDONO QOMNTNDMWIMNMYN KH FA POMNODDMNNOOAN 
T TNT OM OCWOUOWN SRY OT NOMWAONDTENTOTAMNAMNNHNT OMNDVDOM™ONOMWMINAMMMSD 
IND. OO =F DSAANENOHNDNHANNOODN +s FHAIN FS ANMOKDKHANKHNANAYNVNMMN—M—wo 
NQMANNMHMHWOW CEM TMNT HAMAD OONOHONMOOMN OOTFDNODONATHADHAOSM™NO 
MFODWO FAMNNDHNODD—- DON KOVHNDNODMNIN KO +FDOHNANMNSCKKHNTKTODNDDODRHO 
WN FHOAM NNT OTAMNNANANDONE SCM MNYNFOMYTKTNNWATHYTMNODKVNKNDOKAOOMN 
DONDE TOPE DTSCH KTH TH DOHOMNVO SRDS THANAMOYAYAMNNANANANNE OM FMYNAWNOOdAGTDMO? 
®*e2e@« ft @8@ @ @ @ 8@ 8&8 @ ee heh FhUhOHhUCc HmhUCUCU HMhUCUCFmhUC HMC HMhC—~C HMhCUCrh!lhUC MhC<C SC‘ OUCOUPtCU UCU FhUCUC HmhC HLCM GLC H]HhCUCU HMhUC HMhCUCUrFhUCUMhUC PC<CPC<C PSC HCU TLC OhUCU hUCUHhUC WC‘ SCC! 
a, aT 1 ch Oe eee top eon a ioe 
(rrrretrtitrebretveoret) eb’) pes i | 
SCS Se CONN SS SC CeO SCS OE CK NNG KEK KK OOOCOCOOC OOF KK -— OOO 
DDD CV DODD DWV OVO DVDO DOVODVD ODOC DODOVOVWVDDVOODDODDOTOOCOOCO0®O 
(rrrrprerttrrrrrer)tl bd terre trere te bb pred breereeeee ye trtirirree+ 
AQAA AQAA NARQAAQQAAANR AQARAQARQAQAAAARQAQAAAAARAAQAAAAAARAAQARAN 
DFHNADWO NWN DO =F OM AM MNINMYNOPMMNOWBOAHAMNNOSTHNVOMYNOWOSATOKNAOKDOMOOoowM 
DOOM NST OTMM OV — WOOT MM YNTOM AM CWI MOM ANOMNMM™DOONnDODNWOOWOD 
DI YHAP KHON HDONS SON TOATOM ONTAVOPMSCHHYAANATCOHnOVCDOMMNMNooreo 
OUIE™ WNT PFOWNNAMHAMNADW FT NOW T ADNMNOMNMWIMNONTKDNOMNOMNTOraAr rer MHm™ 
DNHNWO DN OONT MANS HNT NN Oem SFDOOONANAS SOMMNSCHNTOAMNESFOTOOMNDHWWO 
DOMNMNOCPMOANNDMO HO ST OT NOOM YONHNT SRO PFONMNMNNATINADM MH NHOMONM 
ANE TTT STOTT CNT TANYA AAAAANNAT OODANANE DTCC KKK THNOMOANKH HN 
Cr Pe | 
CODD DD ODOD OVD DOD ODD DODO DODDOVOVDOADDOAIOOVDVD OV DODVOODODOOVOVOO 
trerete!etd vet ¢ t & @ € € 4 pruitdes 


Se) en OO a eae aren 
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DODO DODO DODO VDD OVD DOD DOD OOV DODO DVO ODOVDDODDOOVDOCO DVWVWOOV0O 
Sa a, RO, ieee naa ene cen OS eee 
e eee eee 8 @ @ ee eeeveevvretreerteereeretreereee ee ere ee eevrerve 
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#198 $9 0-204u7s79009-9-2859000-01 2. pneuzeise 94 
e e a Ve + am « a i e D+ 0 
eG 0.0 -0.3478981D+00-0.487857 1D-01-0. 2103153 D+0 1 


Ae RE A Ke kk ee ee eo ee kak fe 2 ke ok ok oo oo ke oe ako ok ake a fe ofc of kok teat oka feof ok ok oft oe fe teak eo ok ke 


Pood ARE ee oSAt see) WITH THE RESULTS 
THUS shh ean evOuULD eS tiihE tO EXIT TO CMS, 


et PE Sie Ome kT 


(ANY OL neR INE e tee CONTI NUE) 
Be Se a I a aK aK eK 2 ke a koe oo ok ok oe oo ko ok oa ok ok oe keke oR oe oo otek ke ok 


ee. Your Fortran program is now being loaded ... 
eee Execution will soon follow ... 
Px beUITON BEGINS... 


THIS PCRTICN OF THE PROGRAM PLOTS: 

- THE STATES 

- EXTERNAL CONTROL INPUTS, 

~ FEFLBACK CONTROL INPUTS, 

~ STATE ESTIMATES AND 

- RECCNSTRUCTION ERRORS 
FROM THE DATA THAT YOU JUST CALCULATED. 
THE CAPABILITY IS ALSO AVAILABLE TO REVIER ANY 
GRAPHS THAT YOU HAD PREVIOUSLY SAVED AS DATA 
FILES ON YOUR DISK. 


CLE ARM lHE scOREEN) PO CONTINUES 


THE FOLLOWING OPTIONS ARE AVAILABLE: 


Vemeiel LHe DATAw OO WUST CALCU 
Zoe eA CURVE that YOU. PREVIO 


ENE R | ee Ree 


LATED. 
USLY SAVED. 


See wAnLABDES VS Time. 


POU PEROT UP TO .d Sy 
Po tea Ome LOL? 


HOW MANY VARIABLES DC YOU W 
4 


Peer Ore VAntVe Ee DOM YOURE OH TOePLOTOAS CURVE NUMBER 12 


1. STATE VARIABLE (IE: oErO) 
2. FFEDEACK CONTRO if. ts = Get) 
3. CCNTROL INPUT (Ee, U1 Co) 
4. STATE ESTIMATE “(OBS abet ‘tE. HAT1, XHAT2, ETC.) 
5. STATE RECONSTRUCTION R (IE., X1-XHaTi, 
X2-XBAT2, ETC) 
ENTER 1,2,2,4 OR 5 
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WHAT IS THE SOUBSC REP DO yn eo eee Vee ree Le eet 
YOU WANT TO PLOT AS IHE NUBBER 91 CURVE VS TIME? 


WHAT IS THE CURVE LABEVerOre isms rAB LE? 





NOTE: 1. 40 CHARACTERS MAX LENGTH 
2. GREEK SYMBOLS WILL BE PRINTED FOR ANY LETTERS 
ENCLOSED IN PARENTHESES. 
TE. (A) => ALPHA 
EF) => BETA 
F) => PHI 
0) => THETA 
state yl 
WHICH TYPE OF VARIAPIE DO YOU WISH TO PLOT AS CURVE NUMBER 2? 
1. STATE VARIABLE (I=. op Eto) 
2. FEEDEACK CONTECL ith. = + 
3. CCNTROL INPUT (IE Gee ETC. 
4. STATE ESTIMATE haglag® (Ba, “HAT, XHAT2, Ere 
5. STATE RECONSTRUCTION ERROR (IE., X1-XHAT1, 
K2-XHAT2, ETC) 
ENTER 1,2,3,4 OB 5 


KHAT IS THE SUBSCRIPT OF THE STATE VARIABLES Aa 
YOU WANT TC PLOT AS THE NUMBER 2 CURVE VS°iitee 


Zz 
WHAT IS THE CURVE LABEL ECR Sih oo fhe ber. 
NOTE: 1. 40 CHARACTERS MAX LENGTH 
2- GREEK SYMBOLS WILL BE PRINTED FOR ANY Leta 
ENCLOSED IN PARENTHESES. 
IE. (A) => ALPHA 
E) => BETA 
F) => PHI 
Q) => THETA 
state y2 


WHICH TYPE OF VARIABLE DO YOU WISH TO PLOT AS CURVE NOUS reee 


1. STATE VARIABLE {1E., X1, X2, ETC) 
2. FEEDEACK CONTEOL (IE. = = *X) 
3. CONTROL INP (geese, Ucn ET eo) 
4. STATE ESTIMATE (OBSERVER) (IE., XHAT1, XHAT2, ETC.) 
5. STATE RECONSTRUCTION ERROR (IE., X1-XHAT1, 
X2-XHAT2, ETC) 
ENTER 1,2,3,4 OR 5 
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RHAT 2S THiweeSonne OF THE STATE VARIABLE THAT 
: YOU WANT TC PLOT AS THE NUMBER 3 CURVE VS TIME? 
3 
WHAT 2S Toteeunve LAbeleror THIS VARIABLE? 
NOTE: 1. 40 CHARACTERS MAX LENGTH 
2- GREEK SYMBOLS WILL BE PRINTED FOR ANY LETTERS 
ENCLOSED IN PAKENTHESES. 
IE. (A) => ALPHA 
EB) => BETA 
F) => PHI 
C) => THETA 
state y3 


Wpeeetetr PE OF VARIABLE BO YOU) WISH TO PLOT AS CIJRVE NUMBER 4? 


1. STATE VARIABLE (TE, oETC) 
2. FEEDEACK CONTROL ae "§ 3*tcky 
3. CCNTEOL INPUT (Ie. +f, u2, EDC. ) 
4. STATE ESTIMATE (OBSERVSR) (IE., HATI, XHAT2, ETC.) 
SB: STRTE RECONSTEUCTICN EEROR (IE:, X1-xfati, 
X2-XHAT2, ETC) 
ENTER 1,2,3,4 OR 5 


Mone or ne SUBSO REPT OF THe STATE VARL ABLE. THAT 
Es Loum i ANIs 1 © PLOT wAoe tien DER “"S CURVE VS ITILME? 
‘ 
Hiei toeTHE CURVE LABRET FOR THIS VARIABLE? 
NCTEs 1. 40 CHAFACTERS MAX LENGTH 
Ze GREEK SYMBOLS WILE SBE PRENZED FOR ANY LETIERS 
ENCLOSED IN PARENTHESES. 
IE. (A) => ALPHA 
EF) => BETA 
F) => PHI 
Cie 2 TELA 
state y4 
wCUsNayT USE UP O85 READINGS. 
HOW MANY HEADINGS DO YOU DESIRE ON THI5 GRAPH? 
, Oyj lap Rea 
a 


WHAT IS THE DESTRED HEADING NUMBER 1? 


NOTE: 1. 4QO CHARACTERS MAX LENGTH 
Peeorl uke NoOhomntLn Be rRini ko FOR ANT 
PeTiPno eNChCoED Li PARENTHESES. 
IE. (A) => ALPHA 
B 
1 
0 
open loorf systen 


WHAT IS THE DESIRED HEADING NUR? 


NOWE: 1.2 40 CHARACTERS MAX LENGTH 
2- GREEK SYMBOLS Wee Peon weD FORE 
LETTERS ENCLOSE DASime then ts bS. 
IE. (A) => ALPHA 
B) => 
F) => 
O24 
example 2 


WHAT IS THE DESIRED AEA NG SUM ER 37 


NOTE: 1. 40 CHARACTERS MAX LENSTH 
2. GREEK SiYEBOLS WILE BE PRENTED (RO Reaney 
LETT ERS@ ENCLOSED SEN PARENTHESES. 
IE. (A) => ALPHA 
B) => BETA 
F) => PHI 
Q) => THETA 


modern ccecntrol theor 
>> USING A PRE-ALLOCATED DATASFT FOR UNIT FT17F00 1. 
TEE FOLLOWING OPTIONS ARE AVAILABLE. 
1. BEGIN NEW GRAPH OF OTHER CONTROLS, STATES, OR ESTIMATES. 
2. REPLCT EREVIOUSLY SAVED CRAPH DATA. 
3. EDIT THE CURRENT GRAPH. 
4. ELOT REVISED GRAPH ON THE TEK618. 
. QUIT AND/OR MAKE METAFILE OF THE CURVES. 
PREVIOUSLY SAVED. 
SELECT A NUMBER BETWEEN 1 AND 5. 
2 
THE GRAPH EDIT MENU 
1. CHANGE VARIABLES OR ADD A CURVE ON THE CURRENT PLOT. 
2. DELETE CURVE FROM CURRENT PLOT. 
32 EDIT CURVE TITILE(S). 
4. EDIT PASE HEADING(S). 
S. CHANGE THE Y=<AXIS S@uiee 
6. CHANGE THE TIME AXIS SCALE. 
7. CHANGE PLOT SIZE. (DEFAULT IS 8.5 X 6.0) 
8; CHANGE THE LETTERING HEIGHT. 
9. CHANGE POSITION OF THE LEGEND. 
190. EDITING COMPLETE. 
SELECT A NUMBER EETWEEN 1 AND 10. 
9 
HOW MANY INCHES IN THE X DIRECTION 
LEFT OR RIGHT), DO YOU #ANT TO MOVE 
OVE THE LEGEND BOX FROM ITS PRESENT P)ISITION 
NOTE: 1. DEFAULT PLOT SIZE IS 3.5 X 6.0 
. LEFT IS NEGATIVE 
; 3. RIGHT IS POSITIVE 
-4 
HOW MANY INCHES. IN THE Y DIRECTION 
(3p OR DOWN) YOU FANT TO MOVE 
MOVE THE PESEND. BOX FROM ITS PRESENT P)SITION 


Se) 


NOTE: FAULT PAGE Helo Oe 20x 6. 0 
WN IS ee 


iPeecO ol 1 


Z 
E 


UWJN) 
e 


DE 5 
DC tb se 
UP VE 


ry 


GRAPH EDIT MENU 


LoeCnkeA DD AeURVE ON RHE CURRENT PLOT. 
OM CURRENT PLOT. 


SH 


td ‘TORR FAHEY EIS 
PP MW bIOH t eH 


3G? 


HYP RibthtUOQ hh 
td 


nr Cy 
MHiMUHHOWO<S 
Omer this Cp 


OMNMORIt itd +H 
L 


CA Tee 


CALE. 
DEC AULS sls es 5. X- 6.0) 
ING HEIGHT. 
eine LEGEND. 


Zam HO UT” 
W" 


rHOAAAAO ria 

OM mmm mooie 
p> > ft > HH Ft 
Ee PS 
G2 G2G2G1G) 

aOR H Sr OHI 


lee ed 


— 
ty CWOONDNMNLSWYm 


SLING -C 


ELECT A NUMBER BETWEEN 1 AND 10. 


— JU) 
Oo BI 


Pee PLONE NG OPTIONS ARE AVALLABLE-: 


PW GRAPH CF OTHER CONTROLS, STATES, OR ESTIMATES. 
PREVIOUSLY SAVED GRAPH DATA. 

E CURRENT GRAPH. 

VISED GRAPH ON THE TEK618. 

D/O 

VIO 


n GIR) — 
@ ¢ @ «6 6 
tO Whitt 
ChHWhiht o& 
HOKtWa 
Jtqjrgrcy 


Pea, CETAPY ES OF Tit Cuneo. 
USLY SAVED. 


LECT A NUMBER EETWEEN 1 AND 5. 


Hen 
rj 


E FOLLOWING OPTIONS ARE AVAILABLE. 


N NEW GRAPE CF OTHER CONTROLS, STATES, O23 ESTIMATES. 
CT PREVIOUSLY SAVED GRAPH DATA. 

ENT GRAPE. 

GRAPH ON THE TEK618. 

AKE METAFILE OF THE CURVES. 


HES nee rhe WRENN | AND 5. 


In Gthyj— 
ee eee 


CURR 
Sa) 

OR M 
Y SA 


O heOrMhitoD -1 
ty 


UieuW) 
te 
Ee 
bi 


Mar << OM 
Hd 

O a< 

UO titi 

== 

ho 

> ty 


HIOmW SA byt 


QIN 
rd 
mo SO try 
t41+-2O 
Naa 

mbH 
ry 
Utity 
Ono 


WHAT FILE NAME DC YOU WANT THE CURVE DATA STORED UNDER? 
(8 CHARACTERS MAX) 
openlcor 


> BEING EEbeD UNDER 
OL VECTORS GENERA LED 
RODUCT Set 1 SoC Of OA 
ENCES 6 READS: 0 W 


LAT 


E CATA 
-- 26 
FTWARE 
GE REF 


O LOOP DATA 
I PLOT Fries 


GO, CA. 


I PEN 
a) N 2 
iB N DIE 
Ek Pau 


oie 


THE FOLLOWING OPTIONS ARE AVAILABLE: 


1. MAKE METAPILE OP PREV OUsipeoe cows UV. 
23 OU TEE. 


ENTER Vie Gne2 


—) 


WHAT FILE NAME IS THE DADS Stone ee NDE : 


Oopenloop 
THE CURVE DATA IS BEING LOA DED PROMS Pip ree or Oo>) Daath 
>> USING A PRE-ALLOCATED DA TASED SeG RU Neyer eon. 


THE FOLLOWING OPTIONS ARE AVAILABLE: 


1. MAKE METAFILE OF PREVROUSLY "SoyeD l0nE. 
Ze Ue 


ENTER 7% OR 2 


* 

2 

END OF DISSFLA 9.0 == 13197 VECTORS GENERATED IN sl PLOT eee 
FROPRIETARY SOFTWARE PR ODUCT ROR ec Ope one PG, ee 

1817 VIRTUAL STORAGE REPERENCH SO orn ADS (ao nt hee 

DASD 21 DepAcr Ew 

DASD A222) DE Wena 

LASD M20, prance iy 


We I IK ie oe oie aK ee Ae a ea ee ke KE KK aK KK a a oe ok ie oc ok oe oe oe eK KK ok OK ok ok oe ok ek aK ok ok 


DO YOU WANT A VRSTEC PLOTTER SMOOTH COPY OF THE 
LEE Diss Pia Me OIE ivi MOU JUST CREATED. 


ARIK I FE Fe RK I IK a a eae IK ake ok eK oe a Ke ak OK a OK a ok Kk ok ok eK OK ok KK ak eke kK KK ok ok KK kK ek 


x shan Sec 

CTASD 001 LINKED R7G: 7 R/W pies 

Z ene Sa = 205 

DASD 001 DETACHED 

CREAZ INGE ii: 

CREATING NEW F Dies 

PUN FILS 6680 — 20s ns COPY 001 NOHOLD 
TASD 120 DETAGEED 


EEE FR IK IE I EK I aK oe ie a a EK OK 2 EE 2 2 Oo a oe KE aK ke ok aK eK KK oe a KK a ok OK 


YOUR GRAPH (5S) CAN BE PICKED SUP AT Oli COURPUTER SCENTEre 
TES GRAPH (S) WiILESBE ADRESSE a) ee Ore J > heel ee 
Fo ko Rok ok kok kok kok ok dk kk ko Ro RK RK RK AK RK ok ak KK 


FE eo IK IK ie ae IK ok eK KK SK ie OK ie ak eK KE I OK OK i KK OK OK i KK oe KK I KK OK OK KK aie eK KK KK aK OK ok 


DO YOU FANT TO 


1. FUN OPTSYSX AGAIN 
Ze RUN THE FLOT PROGRAM USENG THE SAne eee ero 
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r ore ANG Tie eee ator bi (OPEN/CLOS ED) >} 


—_ @ 


ONCE Ei eee) ee 
tok kk Rk tok ik tok kk a tok lok kkk kk kk kk kkk ak kak kak kk ok ok ak ok ak i ok kak ok ok 3 


Ss 


Se ok Kk ek a eo eo ee ok ok Ko KR KK oe KK KK eK ee eK ok He 


HAVE A GOOD DAY!! 


FO GCG og ola oki lok oko loki kokigiok i lok lok kk ok ak ok ak ak ak ak Jka akc ak 
fee 1.05/17.26 20:08:15 
record off 
EN® RECCEDING OF TERMINAL SESSION 
The graphical output generated by this exanple follows 


@seiligure 3.1. 
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0°S2 





O'02 





asuodsay auty doot-usedo 


JAS — ANIL 
o'st ool 


AYOUHL TOULNOOD NYAAOW 
Go TIA NV AG 
WALSAS dOOT Nad O 


L°e eanbty 


T°0- 
TA GLVLS 


$Z0°0 000°0 $Z0°0- 0S0'0- S2£0°0- 001°0- 
cA ALVIS 


0S0°O 


Oc= 


0°O 


o't 
CA ALVIS 


0'2 


0°S 


0°Ot- 


0‘9= 


O'S 
PA ALVIS 


‘0 
4] 


O°ST 


0°O¢ 


Pe CICSED-LOOP SYSTEB TIME RESPONSE 


The foiiowing clesed-loop system exanple was taken fron 
Wetec, pp 5.8 - 5.19]. 

The full terminal sessicn is recorded below, with user 
input at the left margin in lower case letters or numbers 


below each tt, 


d on 
RECORDING OF 5 oe SE oo LON 
Ce yo.02 20:19:44 


O Stat 


ec 
EG 
’ 

fee 


NM t4jHoO 
ke WH 


UijPeorto lout Are CONRROLS A DREO, OF PROGRAMS : 


OEITSYSX FCRRRAN Seo ial a caaeee 
PTCALC FCRTRAN SauCuDA Le TPE RESP 
TPLOT FCRTRAN DiS ee LASevOr LING RO 


a INFORMATION TO THE NE 


HEE AE ea HE HE he HIKE ee eo He oe i Ee oi eae eae oc ee ake i ok ok ek ee ok ae a ek OK i ke ak He ok a ek OK ok ok 


YOU HAVE A DATA FILE NAMED ‘OPTMAT DATA ON YOUR A 
Peo wide WAS PREVGOUSLY GENERATED BY THE OPTSYS 
FROGRAM AND CCKTAINS THE F, G, H, GAMMA, A AND B 
Paar CES FROM THAT RUN. 


ie our wOULDS LIKE TO USE THESE SANE MATRICES FOR 
ito cUN, Ine GPRISYS  PROGRAMWeWILE READ IN THE 
DESIRED DATA AT THE APPROPRIATE TIME, 


IF YOU TYPE (Y) ES. 


ANY OTHER INPUT WILL RESULT IN THAT FILE BEING ERASED! 


He HK eK oe ok ea ok ek ok ok ok ok ea eae ake eae oe oe eo oe oe oe ok ok ic aK oe oe ae eC oe a a eK ea a Ke aK 
4 


He ke ke ok ok of ote ok ok ok kak oe teak ok he oi ok oie hee oe oe a ae te a oe fe oie ok eo ie oe ake eo ie eK i ak a eK 


LemreCum in De GrrmwUMERiCAL CUT2UT FROM OPTSYSX TO GO 
POSYOUR TERMINAL To BBD} Ghee tOr a Dalok) FILE? 


He te ake ok ke fe ke ake he oe oe ok ake ok ke eae fe eo fe aie coke oe oe eo a a ae ae of ok ak ake ak ee Kk oe ok ie a Ae oe a eK Ee OK 


SS) 
Ke ok ok ok ok of of ok oe oe oe ok ok ok ok ote ke oe ok ook fe ke ke ke ok ke fe ka oe oo oo oe ok ok ke ke ee Rae ok a oe ok eK eK IK OK I 


CUTLU tei DpeecONE 20 YOUR TERMINAL SCREZi. 
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WE AK A KK a a Ke oe KK eK kK ok ok ok kK eK ke KK ok ok ok eK ok ok ok & 


SR COMPLETELY INTERACTIVE OPTIMAL SYSTEMS CONTROL 
M FILL SOLVE NUMEROUS CONTROL PROBLEMS ON THE 
ING TYPES OF SYSTEMS CONTROL EQUATIONS: 
XDOT = {F}*X + {G}*U + {GAM} *(W+H0) 
MEASUREMENT EQUATION-- 
{H}*X + {D}*U + V 
REGULATOR PERFORMANCE INDEX-- 
1/2 * INTEGRAL (Y *{A}*Y + U *{B} ¥*U) DT 
STATE FEEDBACK GAIN DEFINITION-- 

U = -{C} *Xx 
DO YOU WISH TO CONTINUE? TYPE "YES" CR "NOM, 
--DATA ENIRY-- 
SPECIFICALLY D 


TS 
NTERACTIVELY, SEV 
ABLE TO USERS: 


N 
i 


Cy 
tl 


xX 
i 
aie 


N 


® 


bits 


SIGNED TD) READ 
RAL ALTERNATE 


oO ><O 


S 
A 
A 

wr wont AND "GAMMA" MATRICES 
D FROM SEPARATE DATA FILES. 


iE 
A 
METHOD 2--THE "RYN C!, ANDOMGATMA" Meo ES ay ee 
EXPLICITLY DEFINED WITHIN SUBROUTINE "SETUP". 


{NOTE: IN EITHER CASE, 2 1EE USE So hOULDRGET ATTN AMeOrY 
OF THE PROGRAM LISTING AND EXAMINE 
THE EXAMPLES CONTAINED IN S/o sta 


DO YOU WiSH TO CONDE NU eye ES ee heel © i 


2 
E 
1 
B 


Meo 
LA 
VA 
== ro 
care 
on 
ne 


TO YOU WISH TO INPUT THEetr® eG, AND ei GAla 
MATRICES FRCMN SUS ROUTE ED Set Ul ae oe 
METHOD DESCRIBED ON So PREV SOUS Senin 

TYPE "YES@ Ck “ere 
GENERAL OPTSYSX CPYIGNs; 


OFTION 1 ~== “SYSTEM ANAL cio Ie Ur 
CPEN“-LOOP EITGENSYSTEM CALCULATIONS. 


CPTION 2 == SYSTEM OA NAies DS Be CeO CENA LOO 
FIGENS YSTEM CALCULATIONS. 
CPTION 3 =-- CPEN-LOCP RE IGENS yo Seen ou cep 
AND PROGRAM TERMINATES. 
{NPM ADeR Ey ENTRY “FOLLOWS .2 oO Ad Poi 
CETION 4 == MODAL DISTRIBUTION awh Cl Cer ia 
WITHOUT FILTER OR REGULATOR SYNTHESIS 
CR SIE ADS STATE SANA LY Sion 
SELECIVAN OP TION: 1,2, opmeen 4 


43 





eee CU Desthe RWG VAGUS GF STATE AND CONTROL? 
TYPE YE Ss se hae NG 


CES fox LORV CLASSICAL OPTIONS : 


OFTION 1 -- CPTIMAL FILTER AND/OR REGULATOR 
=. Vito Om ie eNO BATERNAL "C# OR %K" 
MATRIX INPUT. 


Cr OMe aeCrl Weer iT ERAND/Ok REGULATOR 
Seve oto EXTERNAL “Cc 
MATRIX INPUT. 

CPTION 3 -- OPTIMAL FILTER AND/OR REGULATOR 
SUN Dae Sl SmwrL tn BXTERNAL Kt 
MATRIX INPUT. 

Cello 4 0—=—" CPLINkL FILTER aan REGULA TOR 
SYNGE polo ent Lena ERNAL "C* AND hk ™ 
MATRIX INPOT. 

onc T AN SOPILON: 1... 2, con 4. 


Peo UMN ony DOM LETERMINESTHE STEADY-STATE RESPONSE 
BOR CONSTANT DISTURBANCE? 


fe ES Me RR NO the 


CO YOU WISH TO LETERMINE THE MODAL DISTRIBUTION 
AND GAIN MATRICES? 


JUNO Ses Boers y eer G15 eh 6 ae 


Cee LCP TRANSFER FUNCTION OPTIONS: 


CPItON te NOMOPEN-LOGr TRANSFERS FUNCTIONS CONPUTED. 
CeO = eeOLro, kool DURS, AND ZEROS COMPUTED. 
Clie eeanmOuLY POLES AND ZEROS COMPUTED. 

GETIONG Sew CNTY POLES@AND RESIDUES COMPUTED. 


Selec) ANS OPTION< 1, 2, 3, OR 4. 


Neto se ANS ron FUNCTION OPTIONS: 


CET@ON 1 == NO NOISE TRANSFER FUNCTIONS CCNPUOTED. 
CEI rONe2e== FOLES Re oOteUES, ANDeZEROS COMNPUIED. 
Ge7 to Nee — SONY SPOLES AND ZSROS COMPUTED. 

CETION 4 -—= CNLY POLES AND RESIDUES COMPUTED. 


Slee lt enkeOrPthONe Ut, 2, 3, “OR@4. 


4 


CONPENSATOR TRANSFER FUNCTION OF aie. 


OFTION 1 <-- NO COMP. TRANSEER ere laons COULUI re 
CETION 2 == FODES, RESDPDUES,. Atweogen es oC ONPUSED. 
CPTION 3 == CNLY POLES SiN DeZEROoeeeMPure D. 

CFIION 4 —-=-"CNLY POLES eae yh eo rs COUP UT aye 


{NOTEs A COMPENSATOR TRANSEER, FUNCTION CAN GBs 
COMPUTED ONLY IF BOTH A REGULATOR 
AND FILTER ARE SYNTHESIZED 
AND/OR INPUG.; 


SELECT AN OGPTION;S \0 zee, CR OS. 





? 
1 
WILL A FEED-FORWARD DISTEIBUTION MATRIX 
("po — MATRIX} EE INPUT ? 
TYPE "YES" OR "NO", 
N 
THIS OPTION DETERMINES THE CRITERIA FOR DECIDING WHEN A 
MARKOV PARAMETER IS ZERO-THE MARKOV PARAMETER INDICATES 
TEE CRDE® OF THE NUMERATOR POLYNOMIAL OF EACH TRANSFER 
FUNCTION. 
ALL "N" ZEROS CF THIS POLYNOMIAL ARE PRINTED OUT AND 
THIS TEST TELLS HOW MANY EXTRA ROOTS EXIST AT Z = O. 
LESS THAN 10.0**{-IE} IS CONSIDERED ZERO. 
TEE DEFAULT VALUE OF THIS PARAMETER {IE} IS 6. 
IN CTHER WORDS, IE = 1.02-6. 
IF YOU DESIRE A DIFFERENT MARKOV CRITERIA, 
TYPE THE INTEGER VALUE. 
, iF YOU DESIRE THE DEFAULT VALUE, TYPE "0" {ZERO} 
0 


DO YOU DESIRE TO SYNTHESIZE A STABLE FILTER {CR REGULATOR 
DESTABILIZING TRE OFWGINAor Sle 


{NOTE: WORKS FOR FILTER OR REGULATOR BUT NOT FOR EOTH 
IN THE SAME RUN.} 


TYPE "YES" Gre vec: 


CO YOU DESIRE TO PRINT DHE SUES LAG RANG eGENOYS oe 
FRICR TO DECEHPOSID TON SirOrReeH SCrees Tht PROGRAMA 


TYPE “YES‘ Grier 


POWER SPECTRAL DENSITY {PSD} OPTION 1 : 


OPTICN 1 --— COMPUTE THE PSD OF THE Cee 


ND/OR T 
CONTROLS OF THE CONTROLLED 5 M WH 


eh HE 
SLlBU WEHEEN FORCED 9am 
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(ed 


PROCESS Aven een one nt NOLS: {NOTE: 
REGULATCR AND A FILTER MUST BE RESIDENT 
PROGRAM 1060s Pain SSOPTiON.} 

OPTION 2 -- SAME AS OPTION 1 ABOVE BUT ONLY PRINT THE 
RESIDUES OF EACH TRANSFER FUNCTION 
USED IN THE PSD COMPUTATION. 

Oe etwn >, — NOT DESIRED. 


SLC LAN FOr EON su 2 eOR 3. 


Bem OUSDE SIRE RECULATROR SYNTHESIS ONLY? 
iy Ere ES Oh NO, 


THE "FY, WEN, WH NGAMN", WAN AND "BE" MATRICES 
FROM 1 OU PREVIOUS OPTSYS RUN WERE SAVED. 


tote ro LLOWIENG OFTIONS ARE AVAILABLE: 
iemUoL Apis tide ont MATRICES AsALN. 
PemUon SELEC TEDSMATRIGCES AGAIN. 
3. INFUT ALL NEW MATRICES. 


DUPER Wiig 25 Ol Se 


NOTE: EACH SAVED MATRIX WILL BE REDISPLAYED 
ieee ior GR, INPUT SEOUSNCE INTERV AL 

AND YOU WILL HAVE THE OPTION OF CHANGING 
ANDIVIDUAL MATRIX ELEMENTS. 


POR OU Wise TO SAVE THE jie FROM THE LAST 
poietOwse USED IN THiS RUN 


NODZEs THE NATROIA WILL BE REDESPLAYED AT 
THE FRCPER INPUT SEQUENCE INTERV AL 

AND YOU ¥ILL HAVE THE OPTION OF CHANGING 
INDIVIDUAL MATRIX ELEMENTS. 


LE Eee Eo) ORM NO kee 


DemrouU wlon YO SAVE PHEeUAN—=MATRIX FROMTHE LAST 
Rin tOeBe USED IN THIS RUN? 


NC te. ine Unlike xk Wh eeeps REDETSPLAYED AT 
THE PROPER INPUT SEQUENCE INTERVAL 

Ripe. OUsu TEL SHAW E THE OPTION OP CHANGING 
INDIVIDUAL MATRIX ELEMENTS. 


i cee LESSOR NOM. 


DOM Ouchi SH @® SAVE DHEP PV SNATRIX FROM 2HE LAST 
PUN 2O BE USED IN THiS RUN? 


46 


H 


A 


LHE 


NCTE: THE MATRIX WILL BE REDISPLAYED AT 
THE PROPER INPUI SEQUENCE INTERVAL 
AND YOU WILL HAVE THE OPTION OF CHANGING 
INDIVIDUAL MATRIX ELEMENTS. 
TYPE "YES" OR "NON, 
Nh 
ENTER THE # OF CONTROLS {NC} OF THE CONTROL SYSTEM MODEL 
{"GU-MATRIX} .« 


ane THE toe MEASURE MENTS OR OBSERVATIONS {NO} OF THE 


; WHM=MATRIX 
\ 

ENTER THE # OF EROCESS NCISE SOURCES {NG} OF THE 
: {"GAMMA"-MATRIX}. 
0 


FLAG/PARAMETER SETTINGS FOR THES RUNSARE 2S fF OLte ee 
TOL 4i¢ iIR @SS 2M) ITF] SEE Zee ee eee eer eee 


0 0 0 0 0 0 0 0 0 0 0 
ISET IDSIAB SSE SD” AX eC Cees onc ee NOR NG 
0 0 0 0 0 1 4 1 4 0 
ORDER OF SisSity— 94 
NUMEER Cr VCON TROLS. — oa 
NOMBERK CF CBSERVATICNS = 4 
NUMEER CF EFROCESS NG€ise sourers — 9 
THE SYSTEM MAT RIX a Re tee 
0.9 1.00000 0.0 0.0 
0.0 =O dal oup =Oe Ota 0.0 
9.80000 =| eu 00 0 = (00a 98 0 Or 
0.0 0.0 1.00000 0.0 


DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
TYPE “YES “(Cree nOt 


Nn 
CEEN LOOP DY NAMI CS Wat Ro co eeeeerenene meters ee nererenen er 
One 0.1000D+01 0.0 0.0 
0.0 -0.415CD+00 -0:1110D-01 0.0 
0.98COD+01 -0.14305+01 -0.1980D-01 0.0 
0.0 0.0 0.1000D+01 0.0 
ENTER THE MEASUREMENT SCALING MATRIX {"H"-MATRIX}. 
DIMENSION = # OBSERVATIONS {NO} X # STATES {NS} 
THE ELEMENT H( 1, 1)= 


4] 


]2])7 OFF OV OV OV =I OF OV OV Or 


On) Ov Ov Or) 


) 


amd * 


ie ELENENT H( 1, 2)= 
PEP LEMENT H{ 1, 3)= 
TE SELENENT H( 1, 4y= 
Peer oehbNT H( 2, °4)= 
Pe ecLEMENT H({ 2, 2)= 
Piet LEMENT H( 2, 3)= 
Ae ELEMENT H({ 2, 4)j= 
Mab fLEUENT H{ 2, y= 
PeEeeLENENT H{ 3, 2)= 
Mic kLEMENT H{ 3, S)= 
THE ELEMENT H( 3, 4)= 
ities ELENENT H({ 4, 1)= 
ie FLEMENT H{ 4, 2)= 
Mabe ELEMENT H( 4, °3)= 


THE FLEMENT H( 4, 4)= 


URE wEASUREVENT SCALING MATRIX [{VH'"=-MAT RIX} 2a. 


1.00000 0.0 OF 0.90 
0. 0 1.00000 0.9 0.0 
0. 0 O20 1.00000 0.0 
0.0 0.0 0.0 1.00009 


POm CUSWisH TO CHANGE THE VALUEVOF ANY MATRIX ELEMENT? 
Dee MVE St aOR UNO MN. 


MEASUREMENT SCALING MATRIX... 2.222 eceee ee eH. 
OCOL+0-1 0.0 0.0 Oreo 
0.1000D+01 0.0 Ch 0 
0.0 0. 1000D+01 OO 
Ono 0.0 0.10 00D+01 
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Ov OF OV OV OV OV DVI OV OV OV Ord 


Ov OF OV OV On 


Ete 
Dae 
hes 


Ae 


Ae 


ne 


BORE 


LHe 


nee 


GEE 


nee 


oH 


THE 


THE 


oe 


OOOO 
ee 6 
OOOO 


BR TRE OULEU wes 


NSION = 
FLEMENT 


EG EME Na 


ELEMENT 


ELEMENT 


ELEMENT 


ELEMENT 


EEE Ny 


ELEMENT 


ELEMENT 


ELEMENT 


ELEMENT 


PLEAENS 


ELEMENT 


ELEMENT 


ELEMENT 


ELEMENT 


OU TEU: 


* 
A ( 


A ( 
A ( 
A ( 
A ( 
A ( 
A ( 
A ( 
A ( 
A ( 
A ( 
A ( 
A ( 
A ( 
A ( 


A ( 


OOo 


1, 


MEASURBMENT COST wha tax 
0.0 
0 


verona Ree 
RTIONS {NO 


{SV AN— TARA} eee 


0.90 9.0 
0.0 0.0 
0.0 0.0 
0.0 O62 5000 


DO YOU WISH TO CHANGE THE VALUE OF ANYVSUATRA. ELE aie 


AL NERS, 


Ww YES Ww 


CP 


WNO " z 


CUTe2UT COST MA TR U3 others es oe ee « oe eeneneeeis 6 Stet 


0.0 


Oe 


0 


U0 0.0 


4g 





Oe0 0.0 0.0 ° 
OF 0.0 0.0 529 
0.6 0.0 Oee 0.2500D+00 


meer IHRE CONTRCISDGSTREBULLON MATRIX {"*G"=MATRIX} . 


DIMENSION = # STATES {NS} X # CCNTROLS {NC} 
THE ELEMENT G( 1, 1)= 


Ms 
0 
E THE ELEMENT G( 2, 1)= 
6.27 
THE ELEMENT G( 3, 1)= 
S.8 
; THE ELEMENT G( 4, 1)= 
2 
THE CONTROL DISTRIBUTION MATRIX {"G"-MATRIX}... 
0.0 
6.27000 
S. 80000 
0.0 
DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
TYPE "YES" OR "NO", 
nN . 
ENTER THE CCNTRCL COST WEIGHTING HAT TRI Le a IX} 
DIMENSION = # CONTROLS {NC} ONTROLS {NC} 
: THE ELEMENT B( 1, 1)= 
131.3 
TEE CONTROL COST MATRIX....<.-e--Beae 
131.30000 
DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
TYPE "YES" CR "NO", 
nN 
THE CONTROL DISTRIBUTION MATRIX..+.-«+--.--Ge-~ 
0.0 
0. 62705+01 
0.98COD+01 
0.0 
TEE CONTROL COST OR cnc ene e@#e#e?e89e08908 @e8 @ ere 6 


Oe lei cet 0s 
EIGENSYSTEM OF OFTIMAL KEGULATOR......... 
G=VoOCcP OPTIMAL REG. E-VALUES...DET(SI-FtG*C).. 
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-1. 23385D+t00, 5. 54546D-015-4. 198350 -00Reetee 2 


C-LCCP KIGHT EIGENVECTOR MATRIX [eee Ce ekeie 
-1.019349D-01 2.3087177D-02 —-8. 155460 p=0 2S 1222 64D-07 
1.129691D-01 ~8.5013480D-02 =-9231733600 0 23070 74-07 
1.000090D+00 0.0 1.000000D+00 0.0 
-6.742684D-01 -3.030447D-01 -2.865351D-01 -7.7484US9D-01 
CCNTROL EILGENVECTOR MAT RIX cc eee eer. 
—~-5.464314D-03 2.109409D-02 2.713925D-02 ~1.676334D-02 
-3. /64753D+00 2.578703D+00 =3. 56230 700-0) ieee cra 
—-3.4Z1605D+01 -9.486653D+00 -4.604269D+00 -3.245261D+00 
3.764753Dt00 —-2.578703D400 1.352512 00. C27 Unt 
1.526581D+01 2.419€63D+t00 1.609198D+00 4.8415 48ND-01 


THE CETIMAL FEEDBACK GAIN CONTROL MATRIX. (eS oe 
-8.54927T=01 -3.29750D-01 9-8. 5345)-072 Bee 


LEE CLOSED LOOP EXNANIGS MATEY < cae ft - Gee 


0.0 1.000000D+00 0.0 Nee 
-5. 360337D+00 -2.4511S7D+00 -5.462116D-01 -2. 73593 1D-01 
1.421803D+00 -4 .612572D+90 -8.561786D-01 -4. 276256D-01 
0.0 0.0 1.000000D+00 0.0 
DO YOU WISH TO OBTAIN A TIME RESPONSE 
OF THE SYSTEM YOU ARE EVALUATING? 
(Y OR N) 
NOTE: YOU MUST EE LOGGED ON AT A DUAL SCREEN 
(TEK 618) TERMINAL TO UTILIZE THIS MODS. 
THE F (SYSTEM), G (CONTRO L) H SOBSERVABLES) ¢ 3AM (NOISE) , 
A (OUTPUT COST) AN (CourR L COST) MATRICES WILL BE 
SAVED FOR R ENTRY TC THE MAIN OPTSYS PROGRAM. 


_ 


eK He HH I A I ee I ea oe a a eae oe ae a oe ee oe oe ke ea Kk OK ek a aK KK ok 


IF YOU ARE DESSATi SERED Iie thE ohio ees 
ZHUS FAR ANESWOULD LiRre Orr eres, 


~TYPE *Y* TOs 


(ANY OTHER INPUT TO CONTINUE) 
He IK He KK eK a ke ee a ea ee ok eK Ko ook oe KK ok oe KK OK ok ek aK he EK KOK OK ek kok Ke ok & 
LOADING ORT Ghia 
EXECUTICN BEGINS... 


DORING THIS SECTION GPS THE 2rOCk Ae Cle eam 


- SELEGaaeaH aa CEe SYSTEM eakesPONSE TO erro 
aoe COO SED LOOP, OR FILTER/REGULATOR ) 
= “FROViDs START NI STOP TIME FOR PLOTTING CALCULATIONS 
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flee 


soe THE TY PrecCreps © Vv ING Seer ED a ee RAM 
FROVIDc START AND STOP TIMES FO kIVING Cit 154 (S) 
FROVIDE DRIVING FUNCTION MAGNITUDE ( 
CLEAR THE SCREEN TO CONTINUE 

LES eh Aor nelee 
Oe 0 1. 0OO00 0.0 OreO 
9.80000 -1.43000 IO eOn oo 0 ORO 
0.0 O20 ne UOO00 0.0 

THE G MATRIX 
Oso 
6.27000 
S.o0000 
Ou 


THE C MATRIX 


Sued 52 =0.3° 2075 -9.08534 -0.04364 
ELE ett tc ey ole 


17200000 0.0 0.0 0.0 
0.0 1.00000 0.0 0.0 
0.0 0.0 1.00000 0.0 
ae 0.0 0.0 1.00000 
THE K MATRIX 
0.0 0. 0 0.0 Ore 
0.0 0.0 AMG, 0.9 
O20 0.0 0.0 0.0 
0.0 0.0 0.0 0.0 


THE FOLLOWING PLOTTING OPTIONS ARE AVAILABLE IF THE 
ree RE PSM ATRECES WERE CALCULATED IN OPTSYSX: 


CPEN LOOP TIME RESPONSE 
XDOT = {F}*X + {G}¥*UC 

CLCSED LOOP TIME RESPONSE 
XDOT = {F-G*C}*X + {G}*¥UC, U = {[C} ¥*X 

OPTIMIZED FILTER CLOSED LOOP SYSTEM RESPONSE. 

XDOT = {FI*K + (G)¥OC, Z = {H}¥*X 
XHDOT = {F}*XH #{G}U + {K}*(Z - H¥XH} 

FILTER + REGULATOR CLOSED LOOP SYSTEM RESPINSE. 
XTOT = (F+#G¥C}4X + [(G}¥*UC Z = ts ee: 
XHDOT = {F}*XH +{G}U + {K]*{Z - H¥XH}, JU = *¥ H 

SaaneT de 2cmore de 
THE AUGMENTED F MATRIX (F+G*C) 
0.0 1.00000 0.0 0.0 
-523€034 -2.45120 -0254621 -0227359 


De 


1.42180 =H a6 1257 -0. 
0.0 a0 i 


8 -0.42763 
0 0.0 


AI WHAT TIME DO YOU WANT 105s 7 
TEE TIME RESPONSE CAFCULAT SONS 


INPUT START TIME IN SECON eee C hoary 0.0) 


0 
AT WHAT TIME DO YOU WANT TO STOP 
THE TIME RESPONSE CALCULATIONS? 
, INPUT STCP TIME IN SECONDS. 
25 
THIS FROGRAM DIVIDES THE TIME INTERVAL YOU JAVE 
JUST SPECIFIED INTO UP TO 500 SMALL INTERVALS FOR 
THE INTEGRATION AND PLOTTING ROUTINES. IN ORDER 
TO SAVE COMPUTER TIME, THE NUMBER OF POINTS CAN BE 
CAN BE REDUCED WITH SOME LOSS IN CURVE FIDELITY. 
, HOW MANY PCINTS DO YOU WANT TO CALCULATE? 
500 


DOES THE SYSTEM UTILIZE A DRIVING FUNCTION (CONIROL INPUT)? 
(Y)ES OR €N)O 


DOES THE SYSTEM START WITH ALL INITIAL CONDITIONS = 0.0 ? 
(Y)ES OR (N)O? 


WHAT IS THE ZNGTIAL CONDD ARON EOE eee 


0502 

F WHAT IS THE INITIAL CONDITION SCR ae 

0 

E WHAT IS THE INTSOAD) CON DMO ON EG Rene 

0 

: WHAT IS THE INITIAL CONDITION FOR X{ 4) ? 

Q 

REE LANGE HAST aQPEOEIBITEE Teas. 
1. SELECI] ANOTHER TYPE CE NS toa Omer es 
(OPEN, CLOSED, FILTER OR FILTER/R EGULATOR) 


Ce STAR PeANDe Ss tor Tits 
ee DRIVING FUNCTIONS 
4. INITIAL CONDITIONS 


5S 





CONTINUE 


Se 
SELECT A NUMBER BETWEEN 1 AND 5. 


Col) 


Zz Ay 


OLLOWING INFORMATIO 
N INDICATION OF PRO 


G4 eq 
FIA 
E1f4 


PLOTTIEL. 


¥ (3) 


STATES AND STATE ESTIMATES CAN BI 


AIL CCNTROLS, 


X (2) 


X (1) 


U(1) 


a 


SEEK ENE EOE NNANNANNAIMANNAMMMMMMNMANMNMNAt tsa saa sq sts 
eT UIMIC ISIC IMIS, 0 171, GHG CIMICISIIE GTC 11 GMOS TIE 
' (yt Prrvb’ovrrerbteertd tte (red (rretdoee 
PA eie A AAAARRAAARAAAMAAnRAAA cacaad ccs FATQVHAR AR AAOAALIANRA 
DIANTMNHOMOTNDODIADNS MMMM NOMONDNNONNHNOTOKNHNOAMNKMNK-DOKTMOONMNMNWYN 
MADE HM KH KHOU TDONWO FO AMNDNVNOISTANSCYNTOTADOFM OMNnAMHOM OM OMrodO 
ARNON OM NOWINODDEOMNMN OO at DN APNO OO SELL Se tT VO NIDA TOO a IMHO 
MOMNANMNCK AN SFKDONDCOl>K— TFOHNDM OVA TONDNDTANOKNDDDANANNOMNOENHAe-ne +t 
OANWI TOOT OW COM COP MEM. K—HOVMO DODOOMNNOMANFHANOMNTAMNST?TSTOIMNENADY —wo 
FHNDY APHKHAVNAMNM OM FeNe DONOO AUN NOV HODNANHOAMMMolcnwonmMmMo 
OTF ON AYE OM TO TAT DENG) CNS KH TDODTMNNNNFTMYHSSCKTHNSEMNHNTeEAPSCKANN 
*e eee © @ © © © 8 &e © we we F © F © &G es Fe eFeeseee*s*se¢# 6 © » © ¢ © © © gg» © © @ © © @ @ »@ 
Fe a 8 OO ae eae aac oe 
i | ited ' | t ' A | 
TET NCCE AAMNANANANMMANNMNMIMMM OY) + FMM) FW ts Fatt se TFTAMNOWMMNNOWONWN 
OOD VO DOD OVD VO OVD AD OONV ODO OOO VC ODVOVCOCOVDOVVDDOCOVONVWOOVOVOOCOOEO 
i’rrrrvrreterttoprrrtrbvrrrtretoroervbrerervertrbtrbvrt’?trrrrterbbrrreereb’oes 
SqjcuagaqgnaaaanaaAaeoannoaanaaaAaAaaAaAaeaeonoAnAaAaanAnaAaaaAaamaAade 
SFANANTFNNANEHRONCR CHM OAMNEK KN TNDOOlH?> FLNDrVODOMNODOTOAWIES DONDMOO 
Dr DI ODO CO SH NODHAMOE FOMDANAMTOPTNOMANSCTNEKHDTNHOCADODO-Trowowo 
OO FAINT HK DNHNONDOAMNOINNHNNONECPEAM MO tT NDEKANAM Al AM HK NODOR WOO 
NDF OMOWMIDON FANKDODNONNEKTNMNNIND OT SFAMNANNNODOEKHNNDWOKNAAMoO 
OD TH MEMO HNONDN OWN DOP OF FDODAVANOVWOMDOMMOO FPOAPLNVNESPKHONSVNOI— 
OtFtOAPFPKNANATAHVANN0 STFOTMNOWNYSFTNVNTKHYASTTAPNO TODAY ArHNOMWCTEAMAt 
OANANT MNF KADAHNHAANOPSCKTKTNOMAANANKP OCC TNAMATANONTKKNMNMNsAANNAK— 
sees. * #® @# « ® @ # %€@ © @ ® @ @#2@¢ 8 8® @ @ 8 ®@ @ @ @» @ @ @ @ @ @ @ © © © @ @ @  @hUmemhmlUCUCOhmUC HUCSDUCSFhUC HMhUCPhUCU Dhl 
ODOC VO DDO DOVO VD OO DVOSC ODVODVDDVDODVDOVVDVDODVDODDOVDOVOVDVDOVOVOOCOCOO 
‘ot (rtr)oed (reeds  orered (roee 
MTA MANNNANNANNMANNANANMSANMNHMMNSAS MONON AWM SS FST TOWMI TT AMNOMNM NWN ODOM Wo 
DOD DOV OVDOD DOVOV DO DODVDOVODVDDOVCVVVOOVOVVO DDO VDOVOVDOOOCOOCCO 
(het rrrt)rrrveubrrrvettvrr)!etrprrbvrrrrrrertrvrerbrrerererrerprverebveree_evoees 
AAR AAR ARARAAAARARAAARAARAARARAARAAARAARARAARAARARAARARQAARAAAAARARARAR 
OMDVYN SPAMYMAWIINAM OO FF TN DWN DWE WDE AND OOS FRATMHMONMHNV CDOMNTANTMNIE-HAOMO 
COMNDNOPM OVD DANAE HNOMAANNANDHOTANNTONNO FOOT NAAMES STROM MOTOMrON 
DTONANDA’NAMN DOF SFOANNS CNS OK AMMNNMNOEVHNVOANAMMNAMHDNTODOWODDAAMK— KON 
OMT ANANOMNNANONSTH ODOM MNFOMN CFTOMMNMMNK-MNOOMSTTWODDWNTA TOMNDMOODNWOMAO 
DOONSCHM™NNTODMMNHMOOKTTOMVOMNDODMMNMNVWONM FONMHBDOMNAMTADODTOVNNoM 
OTM ANOS SNOATFANDNSKNDNHOTP HAT TrOM DODNODHOOTSMNK—MOWMNNHAMNAVAENDOowWoeomM 
ATMOS TFAMNANOWINADN SEK TNNHATOMNMNE NST SClOM TOWN DOTS TVYOFTTMMNNAKSCerKwo 
ee? © © © 8 ee eh Hh FhUhHmhUr H!hUCUHhMhC(UchHMhUCUhT}H OH FH BP eC fF FS & © & © e e © F © F&F © © © © © FF &F © © F&F @ © & © @ 8 @ 
OOD DOD VDOOVD DODO OV DODO OCOD ODVODVOVDODVDVOVDVDODODOVOVCVVCOOOOO 
(rede ee oe rete prbuoed (roeueg 


sre cea si ae! ee Oe ae ae 
ess 8 @ 8 @ °* 6 @ ee 8 8 @ e @ 6 e o * 6 @ @ 


e 
SOSCSCKOSSCOOCOCO OOOO OCOCOOCCOCCOOSSCCOOSOSSoSSS OO SSOSCOOCOGCCG 


OD DD DOOD DOO OOO DODOVDODODOVDOOVDVDOOODODVOVDDODOODOOVDOVO0COO 
DWNODMOMNOMNDMNOM CMON OMNOWMNONODOM OM DNDMNDNHOMODNOMOMOMNOMNONOMOMO 
x.ee¢ © @ © @ ®@ s ¢ @ © @ ® * 2@ 8 @ e* «© 8? © @ @ e ¢ 9» ® ®* @@¢@ ¢@ 8 8 ° #°® ° °® s © © @ 8 ® e 8 @ @ e 69 
OOF AAMMN ST FUWIDWOE LL DDODQNNC OT CANAMMT FMMOVM MS DODHANOCOTEKANNMNMHAAWY 


Me fe sk ok ok ok ok ok ok ok ok oko ek ka ke i eo ok EK Ko a 2 a 2 KK KK KK OK IK OK ie aK ER IK ORK OK OK OK 
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IF YOU ARE VISSATI SFIED V2 ie Eee eS 
THUS FAR AND WOULD LIKE 20mm? Gees, 


“TYPE 'y' TO EXIT- 
(ANY OTHER INPUT TO CONTINUE) 


He AK Ae Ne ok ee oe ie oe ke ee ok ok oe ok ok oe oie ok ok ok ok ok ok ok ok ok ak KK ok oie oe ok oe a ok ek ok kk ok oie Oo kak ok ok ok 
F (120) R/O 
c (121) 70 
EF (122) 70 


eee Your Fortrah program is now being loaded ... 
ees Execution will soon Follow ... 
EXECUTION BEGINS... 


Det eee eh el 


( 
tH 
tdhitd esta 
erPH ON PmdtwWHy 


TES 
CONTROL INPUTS, 
SoMa INPUTS, 


al 
CLlON BRRORS 

AT YOU JUST CALCULATED. 
S 

H 


ALSO AVAILABLE TO REVIEW ANY 
Dee PREVIOUSLY SAVED AS DATA 


AR THE SCREEN TO CONTINUE. 


+4 Od 
ORR POWMDH O 


YY SHES mUwUHhhy 


Cr 


THE FOLLOWING OPTIONS ARE AVAILABLE: 


1. PLCT THE SQDATARYOU JUST (CALC UGR ED. 
2- PLOT A CURVE THAT YOU PREVIOUSLY SAVED. 


ENTER 71 OR 2 





YOU MAY PLO] UP TO]-S SYStEM VARTA) Poe oe fe. 
HOW MANY VARIABLES DC YOU Wish fo, PLoT: 


4 


VHICH TYPE OF VARIABLE DO YOU WESH TO PLOT AS CURVE NUMNB Ee 


1. STATE VARIABLE (TE. X1, X2, ETC) 
2. FEEDBACK CONTRCL ues ~C*X) 
3. CCNTROL INPUT (IE., U1 ETC. -) 
4. STATE ESTIMATE Esa eats (IE., XHAT1, XHAI2, ETC.) 
5. STATE RECONSTRUCTION ERROR (IE., X1-XHAT1, 
X2-XHAT2, ETC.) 
ENTER 1,2,3,4 O8 5 


WHEAT IS THE SUSBS@ RPT Cr THE STALE. Vane eee 
YOU WANT TO PLOT AS THE NUMBER 1 CURVE VS TIME? 


DO 


WHAT IS THE CURVE LABEL FOR THIS VARIABLE? 


NOTE: 1. 40 CHARACTERS MAX LENGTH 

2. GREEK SYMBOLS WILL BE PRINTED FOR ANY LETTERS 
ENCLOSED IN PARENTHESES. 
Hives ALPHA 


A) => 

E) => BETA 
F) => PHI 

Q) => THETA 


state yl 


JVteovweure Or VARTAPIE DONYOU WISH £O PLOT AS CURVE NUMBER 2? 


1. STATE VARIABLE (TE; oETO) 
2. FEEDBACK CCNTRCIL IE. ae = 
3. CCNTROL INPUT (TEs ; “STC. " 
4. STATE ESTIMATE SERVER) (TE., XHAT1, ZHAT2, ETC.) 
5. STATE RECONSTRU TION ERROR (I1E., X1-XHAT1, 
X2-XHAT2, ETC) 
ENTER 1,2,3,4 OR 5 


iio Lie SepSse rip i Ore ine STATE VARIABLE THAT 
Ee MOU! A Nie. PrOn AS LEE NUEBER 2 °CURVE VS TIME? 
g 
Herve ho LHESCURVE LABE PerOGk GeiS VARIABLE? 
mete: 1. 40 CHARACTERS MAX EENGCTH 
Pee one iio lees Ve bora PRONG ED POR ANY LEITERS 
PNCLOSEDL IN PAREN HG on. 
IE. (A) => ALPHA 
E) => BETA 
F) => PHI 
¢) => THETA 
state y2 


NEee@ne) YP OF VARTAPIE DO YOU WISH TO PLOT AS CJ RVE NUMBER 37 


1. STATE VARIABLE {TEs g K1, K2e ETC) 
2. FEEDBACK CONTRC a x) 
PRCONTRCL IN DUT (TES Ot, ie. ET es) 
4W. STATE ESTIMATE (OBSERVER) (IE., XHAT1, XHAI2, ETC.) 
5. STATE RECONSTRUCTION ERROR (1E., X1-xfatT1, 
X2-XHAT2, ETC) 
ENTER 1,2,3,4 OR 5 


Whee USM HE SSUESCROP TeORRIHE STATE VAR IASLE THAT 
YOU WAN LCE OMS SEtEeonbme 3 CURVE VS TIME? 


Gen 


WHAT IS THE CURVE LABEL FOR THIS VARIABLE? 


NOTE: 1. 40 CHARACTERS MAX LENGTH 7 
PeeGrauk cYMNeEGoo WiLL BE PRINTED POR ANY LETTERS 


6 


ENCLOSED IN PARENTHESES. 
IE. A) => ALPHA 
EE)? => BETA 


State y3 


WHICH TYPE OF VARIABIE DOVYOU WiStefo PEOT AS Ci LVE Ne ee 


1. STATE VARIABLE (TEs ETC) 
2. FEEDBACK CONTRCIL (Tf. 6: = =C 
3. CCNTROL INPUT ites sf ETC. -) 
4. STATE ESTIMATE SRSLY ERY (Iz., XHAT1, XHAI2, ETC.) 
S. STATE RECONSTRUCTION ERROR (IE., X1-XHAT1, 
X2-XHAT2, ETC) 


ENTER 1,2,3,4 OR 5 


WHAT IS THE SUBSCRIPT OF Pees stave VARIABLES eae 
YOU WANT TC PLOT AS TRE NUMBER G CURVE VS" Tair 


WHAT IS THE CURVE LABEL FOR Tats VAREASLE: 


NOTE: 1. 40 CHARACTERS MAX LENGTH 
2- GREEK SYMBOLS WILL BE PRINTED FOR ANY DEPRue 
a INS ie 


iE 
F 
C 


VVVV 
ITO OO 


Wo 


State y4 


YCU MAY USE UP TO)3 HEADINGS. 
HOW MANY HEADINGS @RO YOUMDESTRE Ot -fHts Geka 


Op ie O Rees 


WHAT IS THE DESIRED HEADING NUP Eee: 
NOTE: 1. 40 CHARACTER S MAXELENGAT 


2. GREEK SYMBOLS WILL BE PRINTED FOR ANY LETTERS 
ENCLOSED IN PARENTHESES. 
TEL (A) => ALPHA 
EB) => BETA 
F) => PHI 
Q) => THETA 


closed lcop system 
VPHAT IS THE DESTRED HEADINGS NUMNBE reece 
NCTE: 1. 40 CHARACTERS MAX LENGTH 
2- GREEK SYMBOLS WiLL BE PREN@ SD eOr ANY Letiaee 
ENCLOSED” TNE AEN eo or 
i Ewen LPHA 
ie ETA 
i HI 
S Eas 


> A 
2B 
> P 
Pa 





example 3 
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WHAT IS THE DESIRED HEADING NUMBER 37 


NOTE; 1. 40 CHARACTERS MAX LENGTH 
Zs REEK SYMBOLS WILL BE PRINTED FOR ANY LETTERS 
ate PN EAR ENA HE SES. 


MH hy G 


° A) => ALPHA 
Epo eb era 
B= > ee oe 
QC) => THETA 





modern control theor 
>> USING A FRE-~ALLOCATED DATASET FOR UNIT FT17F001. 


tipo FOLLOWING OPTICNS ARE AVALIABLE. 


BEGIN NEW GRAPH OF OTHER CONTROLS, STATES, OR ESTIMATES. 
REFLOT PREVIOUSLY SAVED GRAPH DATS. 

EDIT THE CURRENT GRAPH. 

PLCT REVISED GRAPH ON THE TEK618. 

QUIT AND/OR MAKE METAFILE OF THE CURVES. 

PREVICUSLY SAVED. 


SELECT A NUMBER EETWEEN 1 AND 5. 
5 


in Who 
ee eeee 


DO YOU WANT TO SAVE THE CURRENT GRAPE DATA [IO 
Beeuoeel LATER LO GENERATE A METAFILE? 
Y OF N 
N@ we. of METAFILE FES RE@ULRED FOR SMOCTH VERSATEC PLOTS. 
THERE WILL BE AN OPPOR‘'UNITY TO GENEFATE A METAFILE 
JUST BEFORE EXITING THIS PROGRAM. 
ue 
WHAT FILE NAME DO YOU WANT THE CURVE DATA STORED UNDER? 
F (@ CHARACTERS MAX) 
ciosedlp 
THE CURVE TATA 15 BEGNG EFEEED SUNDER CLOSEDLP DATA 
END CF DISSELA 9.0 -= 16300 VECTORS GENERATED IN 1 PLOT FRAMES 
PwoentetARY SOFTWARE PRODUCT OF ISSCO, SAN DLEGO, - 
loosen SrORAGE REFERENCES: 6 READS: UO WRITES. 


THE FOLLOWING OPTIONS ARE AVAILABLE: 


leeks METAPILE “OF sPREVIOUSLY SAVED CURVE. 
ree leds 


Bellon 10 hoe 


= 


Wea ethene 1S THE DATASSTORED UNDER? 
elosie dir 


ii CURVE DATA I 


Sees LOADED FROM FIL 
>> USING A PRE-~ALLOCA DATAS 


Ere L Os 
Ei PORSUgNT? FT 


acne 
CATED 
TEE FOLLOWING OPTIONS ARE AVAILABLE: 


1. MAKE METAFILE OF PREVIOUSLY SAVED CURVE. 
2 ee OUT. 


BNO ER eOR 2 


Si 


N 1 PLOT FRAMES 
DIEGO, CA. 
ceo 


T A VRSTEC PLOTTER Mee coin (COPY eOre is 
PLA Mee OR ND YOURGUCy CREAT ED? 


US 


Me A i eae cake ee oe oe ea a aK i oe a oc ae oe oe a oe eo oe ok ea oe eo a ie ea oe Kk 


E442 3) 1 LINKED 5/03 F/W BY MVS; R/O BY 0085P 
SASD oh4*6 DETAC 


BED 
CREATING NEWT TE: 
CREATING NEW FILE: 
FUN FIiL<e 6749 TO MVS COPY 007 NOHOLD 
LASD 120 DETACHED 


AK SK ae Ke ee oe ae a ae ae a ae ae aK oe oe a ee oe a ae ae oe aK oe ee eae ee a Ko oe oe oe ok a a oe ok 


YOUR GRAPH (S) CAN BE PICKED UP AT THE COMPUTER CENTER. 
THE GRAPH(S) WIIL BE ADDRESSED TO “POP (USER ID)". 
MK KM KE a a oe ee oe ee ke eke ae oe i oe Ko oe oe ok ke ok i oe i oe ak oe oe eo oi ok ok oo ek ee ok ok ok 


HK Ke eo ee a ea a ee ee i Ke a a oe a i eae ae ek eae oe a ie aKa ie i oe eo i oe i oe oe ae ok 


DO YOU WANT TO 


1. RUN OPTSYSX AGAIN 
Z. RUN THE PLCT PROGRAM USING THE SAME MATRICES? 
: (20 BLOT ANOTHER TYPE OF SYSTEM (OPEN/CLOSED) ) 


ENTER 1, 2 OR 3 
eK HK ke ie oe a eae oo oe ee oe ee eK oe oe ee oe oe ae eae i oe ek eo ok ok OK ok 


5 


Me A ea ee ea oe eae ae i oe ee a oe ea Oe ee ie oe aK ee ie eK a a oo oe i He oe a ok aK eK a Kok 


HAVE A GOOD DAY?! 


fo Ra lok kk gk kkk ik kk kok diakok kkk dkoak kak ak kak ok ak ak dak fk kok kc of ak akc ok ok kook kak ake oe ok 
KR: 729.597 15. (oor oa04 
record off 
END RECORDING OF TERMINAL 3s tooo 
The graphical output generated by this example follcws 


as figure 3.2. 


a 
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OTO' O- St0°0- 


020°0- ¢20 08b0'0—- 


0O¢cO' O0- ¢$L0'0- 


920'0= 090°0- 


TIME 5c 


Closed-loop Yime Response 


Jeez 


Figure 


C. FILTER CLOSED-LOOP Sinv LaATeGh 


The following filter simulation was taken from [Ref. 7 
pp. 332 — 334]. 

In its present cenfiguration, OPTSYSX program sequencing 
requires the input of a [CC] matrix or design of an optinal 
regulator (if a [{G) matrix has been provided), PE1iCr "ae 
initiating the oftimal estimator synthesis or user provided 
IK] matrix evaluaticn. In order to comply with bcualt=a 
program sequencing conventions, and circumvent prcegram 
difficulties which may not be specified in the particular 
system model, optimal filter synthesis may be accomplished 
ty entering the identity matrix [I] in those program input 
sequences reguiring the entry of an output cost (weighting) 
Ho tex 

The full terminal session is recorded below, with user 
input at the left margin in lower case letters or numbers 


telow each "2", 


record on 

EFEGIN RECORDING OF TERMINAL SESSION 
R: T=0.0170.02 2: S540 

optsys 


THE OPTSYS EXEC CONTROES A ERTO OF 2Rece ae: 


1. OPISYS Xi shee TRaAn Stoel iN Pee eee 

2e OPTCALC FCRTRAN CALCULATE TIME aches 

35s OPTPLOT Fert DISSPLA PLOTTING ROUTINE) 
FACH PROGRAM PASSES INF ORMATUON (OS teee tet eG ha 
THROUGH A DATA FILE WRIIZTEN TO Tel esces Dist. 
CASE, THESE FILES ARE “OQPIMNAT DATA AND Ofte LO ae 
THE SIZE OF THESE FILES VARY WITH THE SYSTEM ORDER, AND 
CAN USE_ABOUT. 20%.0F THESUSERS@DI0 SK Sse ACr ae eh Orn 
ENSURE THAT SU FPIGIENT DiSh SP yer Se eer ore. 

- TYPE "EY TO BX ah Y COs Se he Oe ae ae 


WEE AK AK IK ae ae He ee Ha eo ae ie ie ae oie a a a a ek ee oie ie Ae ie a EE Oe OC Ke Ke ae eK eK ie I i eK He ae IK OK OK 


D *OPTMAT DATA ON YOUR A 
GENERA TEDSey THE Olt ows 
EF, G, H, GAMMA, A AND B 


YOU HAVE A DATA FILE NAUS 
DISK THAT WAS FPREVIOUSLY 
FROGRAM AND Gc TH 


MATRICES FROM THAT RU 
IF YCU WOULD LIkF TO USE THESE SAME MATRICES FOR 
THIS RUN, THE CPTSYS PROGRAM WILL READ IN THE 
DESIRED DATA AT THE APPROPRIAIE TIME, 


€é 1 


IF YOU TYPE (Y¥) ES. 


eon Onan LNP OD hibpensocui. EN THAT FILE BEING ERASED! 
ROR ak kok ak ok kkk dk oko kk kkk dk ok kkk kk ke keke kak ok ke ake ke fe ooo ofa stake kkk 


We RK a ae eK i ee a ok oe a oo ee oe a ok oe ok ae oo ke oo ae ke eo ek ok ee ek oo eo ok ke ke 


DC YOU WANT THE NUMERICAL OUTPUT FROM OPTSYSX TO GO 
TO YOUR TERMINAL _ cacramne LOA Dhow) FL LE? 


WK FE IK Re Fe eo I a a ok ok oe a Re ie oi ok a ke a ok ok ok 2 ok ok ok ok ok ke ak i eo ok ok ok ok ok kok ook ok 


Ss 


BO OR RG GK ok ko loko ok doko kk kkk a a kk gk dk kk ok koko ak ke keke oak ote ake ok 
OUTPUT WILL COME TO YOUR TERMINAL SCREEN. 
FO OR OG doi ka lok kk GR OR Jak kkk kok koko kkk ak dike ake kak dick ok akc ok oie ake ok akkak ak ake ok 


HOADENG CPISYS. 


EXECUTICN BEGINS... 

OPTSYSX IS A COMPLETELY INTERACTIVE OPTIMAL SYSTEMS CONTROL 
EROGRAM. IT WILL SOLVE NUMEROUS CONTROL PROBLEMS ON THE 
FOLLOWING TYPES OF SYSTEMS CONTROL EQUATIONS: 

XDOT = [{F}*X + {G}*U + {GAM}* (W+W0) 


MEASUREMENT EQUATION-- 


Z = {H}*X + {D}*U + V 
REGULATOR PERFORMANCE INDEX-- 
J = 1/2 * INTEGRAL (Y *{A}*Y + U *{B}¥*U) DT 
STATE FEEDEACK GAIN DEFINITION-- 
U = - {C}*xX 
DO YOU WISH TO CONTINUE? TYPE "YES" CR "NOM, 
--DATA ENTRY-- 
ALTHOUGH OPTSYSX IS SP ECIFICALLY DESIGNED TJ READ 
ALL MATRIX LATA INTERACTIVELY, SEVERAL ALTERNATE 
MFTHODS ARE AVAILABLE TO HeERes” 
METHCD 1--THE "F"."G" AND "GAMMA" MATRICES 
MAY BE READ FROM SEPARATE DATA FILES. 
METHOD 2--THE "F"."G", AND "GAMMA" MATRICES MAY BE 
EXPLICITLY DEFINED’ WITHIN SUBROUTINE "SETUP", 
(NOTE: IN EITHER CASE, THE USER SHOULD OBTAIN A COPY 
OF THE EROGRAM LISTING AND EXAMINE 
THE EXAYPLES CONTAINED IN S/E "SETJP".}3 
DO YOU KISH TO CONTINUE? TYPE "YES" OR "NO", 
y 


BO YOU WISHWBO INPUT THE “F", “G", AND “GAMNMA™ 
TRE Pebomeureloesvl ROU LINE NSSTUP" IAN THE 


62 


Cader) 


METHOD DESCRIBED ON THe PREwEOUSe oC REEN ? 
DYP.E Uist! OR NO 
GENERAL OBIS Yoo erren 5: 


OFTION 1 -- SYSTEM ANAUYSi Saino. 
CPEN-LOOP EIGENSYSRere CALCULA et Ss. 
OFTION 2 -- SYSTEM ANALYSIS eee N-1TOOLr 
FIGENS YST EM CALCULATIONS. 
CFTION 3 -- CPEN-LOOP EIGENSYSTEM FOUND 
AND PRCGRKAM TERMINATES. 
{"FV-MATRIX ENTRY FODLOWS IN MEDIATED 
CEIION 4 -=— MODAL DISTRIBUTION MATRICES COMPU, 
WITHOUT FILTER OR REGULATOR SYNTHS 
CR STEADY-STATE SANA LY SiS. 
SELECT AN OPT LON «2 Ota 


DO YOU DESIRE KMS VALUES OF STATE MAD Dees roe 
TIPE VYES." sCR NOs. 


OPISYSx LOR/ Class feCAr Or Tiere. 


OFTION 1 == CPITIMA TL Shao re MD / OR eRECULA TOE 
SYNTHESIS Wen NO EXTERNAL YC (OR 
MATRIX INPUT. 


CETITON 2 --— CPURUA Tees ND/OR REGULATOR 
SYNTH Eolseneeren EXTERNAL: “em 
MATRIX INPUT. 

GETION 3 -- OPTIMA PilLIbR AND/JOR RECULATOE 
SYNTHESIS Bieta EXTERNAL “Kt 


MATRIX INPUT. 
CETION 4 -= €PITTHAT Filter eanD 

SYNTH E SISeBerH FAT 

MATRIX INPUT. 


SELECT AN OPTION ee Oe 


ZOR REGULATOR 
PR NAcc! AN Das 


DO YOU WISH TO TEVERMINES DESL EsAD Yao Pale RE Ob ene 
YOR A CONSTANT DISTUREANCE? 


YPEMWYES" CRAUNOSS 


[TO YOU WESH TO LETTER NESTE Sep Ate Da Shere re On 
AND GAIN MATRICES? 


TYPE “YES Sere Ou 


OPEN-LOCE TRANSFER StU NG ACN ero. 
CPTION 1 -- NO OPEN-LOOP TRANSFER FUNCTIONS COMPUTED. 
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cond 


SON 2) —— PORES rec unubs, AND ZEROS COMPUTED. 
TEL TON 3° -— CNTY SPORES SAND ZEROS COMPUTED. 
CPTION 4 -~ CNLY POLES AND RESIDUES COMPUTED. 

SE ECT AWN OD. TON tee tie? 3, OR 4s 


NCISE TRANSFER FUNCTION OPTIONS: 


Vi EON] 1) =— NC NOISE TRANSEER FUNCTIONS 3 OMPUTED. 
C7 TION 2 == PCLES, RESIDUES, AND ZEROS COMPUTED. 
CEE e-—= CNLY POLES AND ZEROS COMPUTED. 

CET 1eN Yu>— CNLY POWERS AND RESLDUES COMPUTED. 


Sie ben “OP TION<- 15 2, 3, ‘OR 4s 


Sour oNonson TRANSFER FUNCITON OPTIONS: 


CE tN ee emeOU Poe TRANSFER FUNGTIONS © OMPUTED. 
CEeeiOn Zia POLES mentollUbo, AND ZEROS COMPUTED. 
SE tONe Sea CNTY POLES ANE Zenoss COMPUTED. 
elton 4 eae y POLES SAND RESIDUES COMPUTED. 
{NOTE: A COMPENSATOR TRANSFER FUNCTION CAN BE 
CCMPUTED ONLY IF BOTH A REGULATOR 
ANU FILTER ARE SYNTHESIZED 


ANI/OR INPUT.} 
SELECT AN OPTION: 1, 2, 3, OR 4. 


mob l A FEED=FORWARD DISTRIBUTION MATRIX 
Pot MATRIX eee: INPUT 2 


ae Gee es: One NO fs 


THIS OPTION DETERMINES THE CRITERIA FOR DECIDING WHEN A 
MARKOV PARAMETER IS ZERO-THE MARKOV PARAMETER INDICATES 
THE ORDER OF THE NUMERATOR POLYNOMIAL OF EACH TRANSFER 
FUNCTION. | 
ALL "N" ZEROS OF THIS POLYNOMIAL ARE PRINTED OUT AND 
THIS TEST TELLS HOW MANY EXTRA ROOTS EXIST AT Z = 0. 
IFSS THAN 10.0**{-IE} IS CONSIDERED ZERO. 
TEE DEFAULT VALUE OF THIS PARAMETER {IF} IS 6. 
IN CTHER WORDS, IE = 1.0E-6. 
IF YOU DESIRE A DIFFERENT MARKOV CRITERIA, 
TYEE THs INTEGER VALUE. 
IF YOU DESIRE THE DEFAULT VALUE, TYPE "0" {ZERO} 
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DO YCU DESIRE TO SYNTHESIZE A STAC ree eek REGULATOR eee 
DESTABILIZING THE ORIGINAL SiS ferr 


{NOTE:WORKS FCKk FILTER OR REGULATOR BUT NIT FCR BOTH 
IN THE SAME RUN.} 


TYPE YES SCra NO 


DC YOU DESIRE ITC PRIN] THe Coes Gh ots 
FRICR TO DECCMPOSIT ION {POR Mee CRrNG hE PROGES 


TYPE MYES Chae Or. 


i 


TE 
M} ? 


POWER “SPECTRAL DENSITY {PSDiv Olle i. 


OPTICN 1 -- COMBUTE THE PSD OF BENE VOU Pome) h Tee 
CONTROLS OF THE CONTROLLED SYSITEN VHEN PORGERS 
PROCESS AND MEASURE EN esse. eee BOTH A 
REGULATCK AND A FILT MUST BE RESIDENT IN Dae 
PROGRAM 10 USE THiS Clemens 

OPTION 2 ~~ SAME AS OPTION 1 ABOVE ee UO Noreen one 
RESIDUES OF EACH TRANSFER FUNCT fen 
USED IN THE PSD Cerne one 

OPTIGN 37=— NOD 2s ere 





SELECT AN OPTION: 17 zee r 3. 


3 

TO YOU DESIRE REGULATOR SYNTHESIS ONLY? 

TYPE "YES" CR "NO", 

Nn 

ENTER THE # OF STATES {NS} OF THE SYSTEM MAI RIX 
, {"F"—MATRIX}. 
Z 

ENTER THE # OF CONTROLS {NC} OF THE CONTROL SYSTEM MODEL 

, ("GY-MATRIX}. 
1 

FNIEFR THE # OF MEASUREMENTS OR OBSERVATIONS {NO} OF THE 
2 {"H"U-MATRIX}. 
, 

FNIFR THE # OF EFROCESS NOISE SOURCES {NG} OF THE 

{"GAMMA"—-MATRIX}. 


FIAG/PAFANETER SETIINGS FOR = lies 2 UN seen CL LO. 
IOL I¢ IR ISS IN “20 FV? se Peer DUE eG 
0 0 2 0 0 0 0 0 0 0 C 
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iS ee IDSTAR IPSD IYU INORM ITREG NS NC NOB NG 
0 0 0 0 0 0 2 1 1 1 
CROPR CE SYSTEM = Z 
NUMBER CF CONTROLS = 1 
Neeser OF OBSERVATICNS = 1 
NUMBER CF FROCESS NCISE SOURCES = 1 
pun (HE SYSTEM MATRE X orem Ae nL x} 
_DIMENSION = # STATES {NS} X # STATES {NS} 
, THE ELEMENT F( 1, 1)= 
0 
jf THE ELEMENT F( 1, 2)= 
1 
. Pero eNeNt F{ 2, 1)= 
0 
, fies leMeENr F{ 2, °2)= 
=. 6 
THE SYSTEM MATRIX {"F"-MATRIX}... 
0.0 1.00000 
0.0 -4.60000 
LO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
TYPE "YES" CR "NO", 
Ti 


Q 


Ord ad) 


0.0 


SeaneweoCr DY NANT CS MarR IX. se sec cee cee cess elon 


0.1000D+01 
- 0.46 00D+01 


Pipe ne one MPA URER ENT SCALING MATRIX {“H!-MATRIX} . 


DIMENSION = # OBSERVATIONS {NO} X # STATES {NS} 
THE ELEMENT H( 1, 1)= 


THE ELEMENT H( 1, 2)= 


TEE MEASUREMENT SCALING MATRIX {"H"-MAT RIX}... 
1.00000 O20 


DO YOU YISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
eee oe eo kh MeO W. 
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MEASUREMENT SCALING MAT RIX...-.---.- ere eeerererll es « 


0.10COLT+01 0.0 


ENTER THE OUTPUT MEASUREMENT COST MATRIX {"A"-MATRIX}. 
DIMENSION = # OESERVATIONS {NO} X # OBSERVATIONS {NO} 
THE ELEMENT A( 1, 1)= 


THE OUTPUT MEASUREMENT COST (Wage ae X} seer 
1.90000 


CO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELERENa 
TYPE ves Che tere 


CUIPUT COST MATRIX.......... ceccce eer cece che o 


0. 10C0D+01 


ENTER THE COCNTRCI DISTRIBUTEON VATRIX ('Go-' 2 eae 


DIMENSION = # STATES {NS} X # CONTROLS {NC} 
THE ELEMENT G( 1, 1)= 


THE ELEMENT G( 2,0 8) 


THE CONTROL DISTRIBUGION Hert Cee Go Y= MAT ED Xe 
0.0 
Oo 7e7 00 
DC YOU WISER TO CHANGE THE VALUE OF ANY OMATRD Ye ELEMEN 
TYPE “YES OR SNe. 


MATRIX upn oMad rx} 


ENTER THE CONTRCI COST ING 
DI Nera SutROLs inc 


WEL 
MENSION = # CONTROL 


‘< 
THE ELEMENT eee ee = 


THE CONTROL COST MATRIZ......2.-.-Beee 
1.00000 


DO YOU WISH TO CEANGE THE VALUE OF DANY SWATHS XE EM Ei 
TYPE "YES" ORs Neer. 
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Memeo tROL DISTFIBUTTON MATRIX... 5.5..... Ge. 
OU 
O.7870D+00 


PaO NRO. COS] (A ORin \ a ctsuape sie are c 6s © siepsue-0. Bee 
0.10C0D+01 


EIGENSYSTEM OF OPTIMAL REGULATOR. .....ce. 
C-LOCE CPTIMAL REG. E-VALUES...DET(SI-F+G¥#C).. 
—1.712C6D-01:-4.59681D+00: 


reer entond “FiGrNVEC TO RMATRIX.. 520 0. simile es 
Paes ooD le 2. 125703D-0 1 
-1.687503D-01 9.771458D-01 

meno Geet GENVECTCR MAT RIX ..< cece wece oso C¥M.. 

meee cl 1SPp-O01 3.957155D=-03 

See ere OF. neGs LEP B-VECT@GR MATRIX. -M- INV. . 
arse ae Zola SD OM 
lee co? SD— 1.062579D+00 


THE CPTIMAL FEEDBACK GAIN CONTROL MATRIX...C=BINV*¥GT*5S... 
reowcCUr+G0 =2.1349T-01 


TEE CLOSED LOOP DYNAMICS MATRIX ......F-G*C.. 


0.0 1.000000D+00 
-7.870000D-01 -4.768018D+00 


ENIER THE PROCESS NOISE DISTRIBUTION 
MATEIX {"GAMMA"—MATRIX}. 

DIMENSICN = # STATES a X # PROCESS NOISE SOURCES {NG} 
THE ELEMENT GAM( 1, 1) = 


THE FLEMENT GAM( 2, 1)= 


Ow Ord 


Loe PROCko so NOlo tml o TRIBUTION MATRIX 
PPGAMM AN =MATREX} . «. 


0.0 
0.10000 


DO YCU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
ii Sees OO NON. 
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TER THE PROCESS NOISE PSD WEIGHTING MATRIX 
ONMATRIX}. 


DIMENSION = 
THE ELEMENT 


EN 
it 


SES fag 


amd 04) 


THE PROCESS NOISE WEIGHTING SATE: ee 
1¢.00000 


DO YOU WISH TO CHANGE THE VALUE OF ANY MATRKEX ELEM 
meee "YES “Ven “NOM 
FROCESS NOISE DISTRIBUTION MATRIX... 3 GARI. 
O20 
0. 10C0L5+00 


FOWER SPECTRAL. DENS IY o—. PROBES Se 1S) eee eee 
0. 1000D+02 


ENTER TEE MEASUREMENT NOISE DISTRIBUTION MATRIX (“aA hee 


DIMENSION = # OFSERVATIONS {NO} X # OBSERVATIONS {NO} 
THE ELEMENT R( 1, 1)= 


2 
0.0000001 
THE MEASUREMENT NOISE DISTRIBUTION MATRIKe ss eee 


0. 00000 
DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEHENT? 
TAPES VYES CR NOiae 


Nn 
FCWER SPECTRAL DENSITY-MEASUREMENT NOISE..R.. 
0.1000D-06 
ENTER THE FEEDBACK GAIN ESTIMATOR MATPIX {"K"-MATRIX}. 
DIMENSION = # STATES {NS} X # OBSERVATIONS {NO}. 

, THE ELEMENT K( 1, 1)= 
$5.4 
, THE ELEMENT K( 2, 1)= 
4561 


THE FEEDBACK GAIN ES ot MATOR SMA Te a eee 
Ge. 0000 


Go 


4561. 00000 


DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
Dire “YES “OR Onis 


n 
fever olLOADY STATE GAR NSees «cc cco ec wocccs Kewe 
G$.5400C0D+01 
4.561000D+03 


TH 
—-9.540000D+01 1.0000C00D+00 
4.561000D+03 -4 .600000D+00 


Poh Noy Slee OF OPTIMAL “ESTIMATOR. « «cee ee 
See ee SUBOPTS EST. E=-VALUES...DET(SI-F+K*H).. 
porouOUODFO 1, &§.99984D+01: 


See moron EIGENVECTOR MADR sis .cc- Moss's 
pee I D- Os — 1.0962 16D-02 
1.00000CD+00 0.0 


eae cent EIGENVECTOR MATRIX.....H (BAR) *M.. 

Eee ao ODS — |. 962 16 D=0Z 
G>EQCP SubBOPT...¥ BPa BoVeerOR MATRIX«-N-INVs« 
O20 1. 
moet 22 9Z0+C le 9. 
BE COVAREA 


at 
De D0 5D—-06 
71090D-C4 


mat) 2% 


Ua. 0 
Zone 0 

PevomeveatUSS OF THE ESTIMATICN BRROR ww ca cccc ss 
2-674042D-03 1.0868C8D-01 


DO YOU WISH TO OBTAIN A 
OF THE SYSTEM YOU ARE 
(Y OR NW) 


NCTE: YOU MUST BE LOGGED ON AT A DUAL SCREEN 
ClLetsosleo [Fh A TL TOsUTILEZEe THIS SMODE. 


MoE F ({Srsaee4), G (CONTROL), 
A (CUTEUT cee AND B_ (COQ 
SAVE De rORSR EE st 


H (OBSERVABLES), GAM (NCISE), 
TROL COST) MATRICES WILL BE 
E 


N 
Poe eeor TSsy5 FROGRAN. 
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Me eH Re ee eK oe ee eR Re eo eK eK KK KK ek Kk ok & 


IF YOU ARE DISSATISFIED SWI Tei eee 
THUS FAR ANT WOULD LIKE TORE Git cee, 


“TYPE "'Y* TO 3hwE 


i 
OTHER INPUT 20 CON Un 
% % HK 


a 
* He Rt a ok kk HK eK KK KK KK 


A 
HK aK kk 


N 
ee 
LOADING OPTCALC. 
EXECUTION BEGENS 


DURING THIS SECTICN OF THE PROGRAM YOU WILL: 

SELECT THE TYPE CF SYSTEM RESPONSE TO PLOT 

OPEN LOOP, CLOSED LOOP OR FILTER/REGULATOR) 

ROVIDE START AND STOP TIME FOR PLOTTING CALCULATIONS 
CI THE TYPE CF DRLV ING gUFoE ue} (STEP OR RAMP) 
I riVING FUNCT IO (S) 
IDE DRIVING FUNCTION MAGNITUDE(S) 


CLEAR THE SCREEN® TOSGCONT DAU 


THE ETAT 


0.0 1.00000 
0.0 -4.60000 


THE. G MATREX 


0.0 
0.78700 
THE. -C. yagi ere 
= (UeOoe =O 21289 
2s. BAER 
1.900000 0.0 
THE K MATRIX 
92.40000 
4561.00000 


THE FOLLOWING PLOTTING OPTIONS PARE AVAL EA Se eee 
REQUIRED MATRICES WERE CALCU Ae OO eyo. 


1. OPEN LOOP TIME RESPONSE 
XDCT = fF}sxX + fOuFUC 


2. CLOSED LOOP TINE RESPONSE 
XDCT = {F—-G#Gher + Gonpetiom Mum Mercere x 
3. OPTIMIZED FILIER CLOSED LOOP SYSTEM RESPONSE. 
XDOT = {F]*X + {6} #0C, Z. H} *X 
XHDOT = {F}¥XH #{G}U + {K}*{Z - H*XH} 
4. FILTER + REGULATOR CLOSZD LOOP SYSTEM RESPONSE. 
XDCT = _[F4G¥C}#X 4 (G}*UC Z = H 3% 
XHDOT = {F}*XH +{G}U + {K}*{Z - H*¥XH}, OU = {C}*XH 
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SELECT oe, 2A, cones 


3 
THE (K¥H) MATRIX 
95.40000 0.0 
45€1.C00000 0.0 
THE CCMBINED SYSTEM F MATRIX (2*NS X 2*NS) 
0.0 1.00000 0.0 0.0 
0.0 -4.60000 0.0 0.0 
95.40000 0.0 -~-95.40000 1.00000 
4561.00000 0.0 -4561.00000 -4.,60000 
THE AUGMENTED G MATRIX (2*NS X NC} 
0.0 
0.78700 
0.0 
0. 7E700 
AT WHAT TIME DO YOU WANT TO START 
THE TIME RESPONSE CALCULATIONS? 
; INPUT START TIME IN SECONDS. (NORMALLY 0.9) 
3 
AT WHAT TIMF DO YOU WANT TO STOP 
THE TIME RESPONSE CALCULATIONS? 
: INPUT STCPE TIME IN SECONDS. 
0.3 
THIS PROGRAM DIVIDES THE TIME INTERVAL YOU HAVE 
JUST SPECIFIED INTO UP TO 500 SHALL INTERVALS FOR 
TEE INTEGRATION AND PLOTTING ROUTINES. IN CRDER 
TC SAVE COMPUTER TIME, UMBER OF POINTS CAN BE 
CAN BE REDUCED WITH SOME LOSS IN CURVE FIDELITY. 
; HOW MANY PCINTS DO YOU WANT TO CALCULAL E? 
500 


DOES THE SYSTEM UTILIZE A DRIVING FUNCTION (CONTROL INPUT)? 
(Y)ES OR (N)O 
TWO TYPES CE FUNCTIONS CAN BE USED AS DRIVERS. 
1. STEP INPUT 
2. KAMP INPUT 


ENTER VouUcm or Be lLON one ero nRe DRIVING FUNETION NUMBER 1 
5 


oni YOUsDE sing INPUT NUMSER t fT) START? 
ive Ole tye Stan. Stine DIN SSCONDsS- 


GZ 





0.0 
eC 
2 


& 
& 


Nie Ee 
W 


a 
ORTU 
OLLOWING AREAS. 


AXIMUM VALUE OF 
ION NUMBER 1 
(Y) ES OR (N) 0? 
PP 
2 *F 


M 
Gi 
IS THE INITIAL CONDITION FOR £RAT (2 


STARI@AND@S20P Tings 
DRIVING FUNCTIONS 
INITIAL CONDITIONS 


CONTINUE 


6 
WN 
A 
Zz, 
© 
O 
fr] 
WY) 
>, 
aa 
tx] 
= 
ro 
Sa 
Ay 
© 
Sa 
Wn 
wey 
Hs 
ES 
EH 
= 
Ay 
a 
FA 


oe 
4 
Se 
SELECT A NUMBER SB E2wW Ee eae. 


WHAT IS THE INITIAL CONDITION FOR ara) 
ae 


WHAT IS THE ENITIAL COND? Prone CR Aaa 
WHAT 2S THESENITI AL COND TICNe eRe ( 2) 


WHAT 


Cre 
Ay 
O 
{H 
Ww 
- 
ct 
A Sead 
AG 
fx) 
mM 
= 
ee | 
= 
EX 
— 
Ay 
= 
tH 
FY 
RG 
be 
WY 
ft 
A 
— 
O 
> 
O 
=) 
fg 
= 
tH 
[A 
| 
KG 
ee 
= 
be 
<< 


DOES THE SYSTEM STARI WITH ALL LNDTTLTAL CONDMEIONS 


Ow 


Wi 
DICATIONVOPSEROPER PRO 


Na ENG INFORMATION IS PR 


E 
A 


All CONTINOUS, 


fr AG 
C4 by 


BE PLOTTED. 
AGS) 


X (2) 


X (1) 


STATES ANDERS TATE ESTIMASES CAR 


U(1) 


ee 


— 
- 


= OOOOC®M 
DODDOOO 


ODOM rem 
TMM IFTMNO 
MNDMQHNOMM 
~Ororor-~ 
OVUNONAAN 
OF Or rere 
e © ¢ @ @ 8 o 8 


OOOOODOCOO 


OOQOo@oco 
ODdOXVdOoOVOO® 
+-+eee+et te + 
ARARAANA 
OFNDWOAO 
OF MNDOOO 
OMmOT Te" COW 
DONMHAOMNW 
ODT FOMNNs 
AMON OO 
WIFOMNNNAN eo 
es 2«¢ @ © &@ ee 8 @ 


ODOOOOO©O 


DOOCCOOCO 
Solelelelelolo 
tHe tetet 
AARARQAARAR 
DWCONOOOD) 
OOrMN~NO 
OTDAM™AWC 
CONM™OMM 
ONMNMHDOrN 
QOO°0Oorers 
eee 
eo ee © @ @ @ 8 8 


DOOODOO0OO 


AANNAANN 
DODOOCO 
++eetetee 
ARAAAAAQ 
DOODOOO 
DCOCD0O 


DOOOOOO 
OMe ee ee 
oeeee#ees® 
DINWDO ADO 
peebvrue 


mT ONNM St 
ODOoVOdO00O0MO 
ee ee 


COODdDDOOO 
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DODOO00 D009 QO CCCP SCT TST TTC CST CTT TTC NANNANST Cee eee 
SS DOOoOOOO OOO OOOO OO COC OOOO O OOOO OCOOOOOOOOOOOOO 
+Heeteeeeeeeteer Peer berrrrrtetr ee er rerbrrerbtrebrrebtoebreunre 
ARANRNRARAAAARAARAAAARARARARAARAARAARARAARARARRARAARANAARARARAR 
ONO FANDHOOMN DMNoQodQwNT MAW TTIMSTODMMNAWMIMOHNOOTOOAWNO 
T= TM DOD SCMHMNOM DMN SPE KH KOO TANNMNANHD&HANKNN DD AFANMOWOM To 
TNODNHAHDOODOODN™ MOM SF FOWONATACAIM OTOT ME OMA TEM OT AMO 
PMP TM OM OM DODO TM FATE KTODDMNN OUT NOW ODMNUIM KM ONO™Y 
THOON -ODOANONVOOM AFHANMMOVTANHWNSTSCKHKTAODTO Or Korwoaowy 
AAT TH THT OOOCOOCODNNANONHANTANDMAOMAMMEANOOHMAMOOMDAIONM 
TT NON NDD OOM ME OOM SS AMMNAN ST STEER KANANMM aN 
eeeeeeee8ffeeee?e¢®#eeeesFfe?e¢¢ @ #e «© © @ «© © ¢@ © © @ e¢ @ © @ ee eH 8 
DOTODDDODDODDVDOVDOVCDODVDOCOOCDO TOD DODO ODDOOOCOOOOCOOO0O 
titt tt btoee 

(SS SOCIO OOO OOOO OOOO OOOO OOOO OCOOCeer Ts rt 
DDO DVO DODO DODVDO DO DOCCODDOCCOD OCODOOCDODOOODOOO000 
EEL Lr eR Er He Ee EE HEH He Hr HEH HEH HEHE HEHE HEHEHE HH HEH HH HHH HOH HH + Ft 
DARRARADAARGRARARARAAAARAARAARAARARARAARARARAARARARAARARAA 
COT UW CON SOD FAS ODN ST TH PODMOOM HK FON WIM KHON asrHNODOON™ ANWO wormed 
DDE OVW DO ANNI TAM AMM E CODMONTHRONOE FSMMONM NINADOW MAAAWO 
WD FOAN ST DDODOOTMNST NT TNMNOMNANDOWIO MOM OMAN T AMM ODow 
M™ KH MNOMNN FF ADON DNNDANKH NOC DOr MOD AaMNNMD DMO ArNAOMNHA => 
DO™NADAMO DNA NAM ANODOKN KOM FrCOUCTMNDWN FOP E KE KM DONMWN 
TMMNDNWO FONOODMPM WME DMN ODT APHAINANSTORDeMNNY HOOCCOCOTTTN™ 
WO TOE STN AN OM & TF TMI WOO MOWOPE EE EEE CH WH CH DO DNDIDHDD NV 2 ff TET ee 
ee eee¢e¢?e?ee 8 @e¢e¢e? 8 ¢ #¢ @¢ @# @#@ @ @ ¢ ¢ ® @ @ @ @ @© @ @ @¢ © @¢ @ © &@ © @—hUmHmUlUMHCUMD 
DODD OVCDODODODO DDD OVDO OVD DOD OOD OODOODDVTODODVDOOOCOO 
i)retvtrtpretet tet tetretbettt)teretttetetvtb)_etrettt tt bebe 
See aoe OOOO OCT TS SCS CSS NNN OO ee 
DDDD VCO DDODVDODODO DDO DOOD DODD VDOVDOOVDOOOOOVOO00O 
+++eeeeeetetteteetiurttirtetbtte t_bbrebretv t bbe bbe bet t breed 
AA ARARAARAAAARAARARARAARARAARARARAARARAARARAARRAARAARAARAARAAAARA 
DOT OMNNDOMWN FHOMYO KH DOMNVHDNWMNMNNMIMNTHNOANANTMNANNOHASFOOCO 
OM DOTmMMMN NOMNDOT FNOAMHYANDOOST NOM OTT KDOANNNOMNAANE 
DNOMNAW AN QDODMODMMYMN KH HAPONMAMNOSTFHNOMNSFDOVOTKYOAMNMNOn 
Tm OADM FMW FORM FOP ANAM DOOMNANTOUWTNADNY™ OO FMM KONO WH 
MAMAN TOOMWOIMN TNO SFAHNMNNANOCOTAONS eRe erTtoTyOOor TOoOTWDOW 
FFF TH HTH KH OCD OODDNNADON DIT NOMODMOMENOOMMIAMOMAOM 
TTT TT NANO DOME OOM ST IMNYANANAT CACHE KNAMM AWM 
ee ee ee er rr rr rr 2 ee 
DOODDODO DDDOO DOOVOD OVD DODO DOVOTVDDODVDOOVODODODOOOOO00OD 
iiretee*!e!let 
NAANANANANANAAAANANANANAAAAAAANAAAAAAANNAAAINANAAAQAANAAG 
DODDODDVODDVVODVDOVDDODDOOVDOVDWVDOV OVD VOVODVOVOOVOOO0O 
++teeteteteteteteete tee eee eteeeteeteeeeeeteeteteteeteeterretet+ 
ARAARAARAARARAARARAARAARAAARAAAARAARAARARARADNAARAARAARARAA 
DDDD ODO DOD OVD ODVDDV DD VDVDOVDDOWDOVDVODOVDDODOVODVDOVODOOO000 
DODDODODD OD DOD OVODVDOVDDOOVOVDOOVDMDODDODODOOCOVOVDVOO0O00O 
CODODOVD DVO ODDD DODO OWDD ODD ION DOOD DVOVDOADDDMDOVOODVOOOCO 
DDODDOD DO DOC DVDDOD GCDO DDO DO VDOODVDDOODODODOVDDOVDOODVDOOOOO 
DODODVD DDODOCODVDOOC TDD ODDOV ODOC OVODVOODODOOVXVOODO0O 
DOVDODODODWDVDOODVDDOVWDD OVC VDVDOOVDVDODODVDDOOVVDODOD DODDOODVOOO00O 
Sie ae oe oe oe Coe ee Se Se eo ee eo eo SS ee ee eS Se fee es ee SS 
er 2 
OOOO OOOO OO OOOO DODD OO OOOO OOO OCOD OCOOOOOOOO0OSO 
rPrr)edtbitevwrt)ttttttettttett t tt bbit tt t tt tt tre’t!)tvt tet 


MNOS DO DONODOSCTAMM FS TWIOVOE (EH DNDN TAAMM aszmMnowl DOANANGE 
DDOMO0 O00 CTCF ST TT TT TOTO MMO NNN NANNANAANANANANANAANACING® 
*eeevereeeeeereeevreereesereeree tere eere ese ewereeeree ewe wee we we we we 
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THE PROGRAM PLOTS: 
74 


Om EX Lh 
= LHE STATES, 


aw 


aa 


‘yt 


(eS OTHER tNeGT=tO CONTINUE) 
FOO a kok kiki ok ak kak ok ak ke kok kok oak kok ko ok ok ok ok ak a kk ok ok took cok kak ak akc akc ak ake ok oe akc oka ok cok 


execution will soon follow ... 
O 


=e 
PXSCUTICN BEGINSsa. 


YOU ART Slims Cott ep NETH THE 


THUS FAR ANT WOULD LIKE TO EXilI TO CMS, 
Tonia OLtLan, program 1S HOW being loaded ... 


THiS Se OR Ee 


eK He eK eK KK kK eK oe oe 2K eae ake 2K ie aK ak ok 2K oe eK a a a a oe eK a ak eK OK 
ee 


eS 
Waar, 
S 
LCULATED. 


CA 
LABLE TO REVIEW ANY 
SLY SAVED AS DATA 


THE FOILOWING OPTIONS AKE AVAILABL 


Es 
1. PLCI THE DATA YOURS USie ceive 
2- PLCT A CURVE THAT YOU PREVI 


ENTER Wen 


ULATED. 
SUSEY SAVELe 


YOU MAY PLCIT UP 20 4 SYST Eiy ss sets Vo el tae. 
HOW MANY VARIABLES DO YOU WISH TO PLOT? 


@ 
“~ 


ae 


WHICH TYPE OF VARIAFPIE DO YOU WISH TO PLOT AS CURVE NUMER 


TE VARTABLE (lea x2 er) 
DBACK CONTRCL (If, 6 = =C#x) 

TROL INPUT (IE., U1, U2, ETC. +), 

TE ESTIMATE (OBSERVER) (IE., XHAT1 KHAT 2, Ene) 
ATE RECONSTRUCTION ERROR (IE., X1-XHaT1 


ENTER 1,2,3,% ORS 


int toh) 
6 6 6 e 6 


WHAT IS. THE, SUBSCRIPT SOE STARE STALE VARS Pear 
YOU WANT TC PLOT AS THE NUMBER T°CURVE VS Ie 


ond 


WHAT I5 THE CURVE LABEL FORM Aro eVAnRLASEE: 


NOTE: 1. 4O CHARACTERS MAX LENGTH 
. GREEK SYMBOLS WILL BE DRINTED FOR ANY LETTERS 
ENCLOSED IN PARENTHESES. 
IE. (A) => ALPHA 
E) => BETA 
F) => PHi 
< Cc) => THETA 
angular positions —st 
WHICH TYEE OF VARIABLE DO YOU WISH TO PLOT AS CURVE NUMBER 2? 
1. STATE VARIABLE (IE., X1, X2, ETC) 
2. FEEDBACK CCNTRCI (If., 6 = -c¥xX 
3. CCNTRCL INPUT (1E., Di, U2, ETC.) 
4. STATE ESTIMATE (OBSERVE fr) {IE., XHAT1, XHAI2, EC.) 
5. STATE RECONSTRUCTION FRROR (IE., X1-xXaT1, 
X2-XHAT2, ETC) 
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angular fosition estina 


Q 
filter only closed lcop 


exameE 


ENTER 1, 2,.3,0505 05 


Vit tomth Pe oleseniel sOr The STATE ESTIMATE THAT 
YOU WANT TC PLOT AS THE NUMBER 2 CURVE VS TIME? 


Hiei tS THE CURVE LABBE POR THIS VARIABLE? 


NOTE ia OO CHARACTERS Clink LENGI A 
Zee okenk  SYMEo Wo Seon Lie Se OR wANY LETTERS 
ENCLOSED IN ReNtoibsto. 
res WA HA 
E A 
a 
T 


IHhuUDYS 
Mn hittWdN 
—~ HHH oS 


A 


B 
al 
2 
> 
> 
> 
= x)el 


t 


BOUSNAY USE UP TO 3 HEADINGS. 
HOW MANY HEADINGS DO YOU DESIRE ON THIS GRAPH? 


OF ORS 


WHAT IS THE DESIRED HEADING NUMBER 1? 


NOTE: 4%. 40 CHARACTERS MAX LENGTH 
Pca SY NbOLo WIEieoE PRENTED FOR ANY LEITERS 
EPNGLOSED Th PAREN@HMSE S. 
It. A) => ALPHA 
E) => BETA 
F) => PHI 
=> THETA 


Wied to LHe DESTRED HEADING NUMBER 27 
De Iie Ac 4Q CHARACTERS MAX LENGTH 


REEK SYMBOLS WILL BE PRINTED FOR ANY LETTERS 
ee ee oN tA es 


A 
B ETA 
; 

de 4.1 


B 
Eee 
T 
oes THE DEBSwhED HEADING NUMBER 37 


Mrd@ 


VVVV 


HETA 








NOTE; 1.2 40 CHARACTERS MAX LENGTH 
Zeecuoune ott bolo Nii BE PRINTED FOR ANY LETTERS 
PReGLOscED IN UPAR ENTHESE >. 
IE. (A) => ALPHA 
E) => BETA 
F) => PHI 
Q) => THETA 
apt rtal COnebe! Syst ens 
iniGu A  PRE-~ALLOCAIEDSDA TASER PORSUNDISFTI/FO) 1. 
E FCLLOWING OPTICNS ARE AVAILABLE. 
FEGIN NER GRAPH OF OTHER CONTROLS, STATES, OR ESTIMATES. 
RFPLOT PREVIOUSLY SAVED GRAPH DATA. 
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SELECT A NUMBER BETWEEN 1 AND 5. 


tJ 


FAPH EDIT MENU 


FS OR ADD B&B CURVETON THE CURR ENDYP ler 
CM CURRENT PLOT. 


Ha 

tc} 
UMUHHQWO<4 
Nome mthm dcp 


SCALE. 
ES DEFAULT IS 8.5 X 6.0) 
HEIGHT. 

OF THE LEGEND. 


ONnWOMtY? wim HF 
to 


HOAAANAANHRHONON 
om We oe os i i oe 
bed Ss bs Br fo FE HY 
HA Sa Z2Sryqnia 


a 


G 
MBER EE 


tr} DOANE WNo 


@) 
ae 


Ie 
E 
C 
IE 
iE 
if 
re 
E 
fe 
N 
M 


BS HOAAMAGM 
4 Se ONE a ame 
= 
Ps 
H 
— 


EN 1 AND 10. 


tJ 
Ee anh 
= 


CN WHICH CURVE CC YOU WANT TO CHANGE THE Y-SCALE? 
ENTER CURVE NUM BER = ee eo cee 


— *4) 


WHAT IS THE NEW Y-MIN VALUE AT THE ORISIN? 
SO 
E WHAT IS THE NEW YeuaxX VaeeeeE : 
Oe le 


ty 


ADD A CURVE ON THE CURRENT PLOT. 
URRENT PLOT. 


HIG) 
> Cc} 
ton UHHH <4 
TT <SH OS 
rOht4 HK 
(J PHM NHI Pa 


E 
‘le omeS « Dik 0) 
Hive 

GEND. 


as 
OW OAINDMNIE Who 

6 666 66 eetseeesee °e6 
HOAAAONtIrMOION 

ey ales paket go Coe od Gee or 

PA Se bree Pt Bs Fs FF HH 
Hse A HA & 
HMA QAGQ 
2tdtdbthihitd Ohi 
Q0 Mate ht <p 
OMWWOtihy tito FH 


os EL 


GQ) 
2 mWOHnHerPoHtibt © 


SELECT A NUMBER BETWEEN 1 AND 10. 


10 
E FOLLOWING OPTIONS ARE AVAILABLE. 


N NEF GRAPH OF OTHER CONTROLS, STATES, OR ESTIMATES. 
CT EREVIOUSLY SAVED GRAPH DATA. 


GI 

PG 

Te Ee UR APH. 

an ON THE TEKOUG. 
at 


H 
2 
iE 
I 
U ArFLLEVOE Toe CURVES. 


z 
5 
ieee 
on E 
4. £F 
DO E 
ie 


J 
= 
a) 
mS FO 


ND 
REV BOU Si yas 
Se A NUMBER B 
mi 


EEN 1 AND 5. 
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UnredyM 


aa? JU) 


ti 


a 


Rj 


>) 


WN WN 


ub 


ci ot 


hd OOO WNOANNIE WN 


trYH OwoOd~1%0HnVNLE WH 


O 
am 


IFE FOLLOWING OFTIONS ARE AVAILABLE. 


PVE ohne Oh wOtitmeONTROLS, STATES, OF ESTIMATES. 
FREVIOCUSLY SAVED GRAPH DATA. 

HE CURRENT GRAPH. 

EVES EDeGE Ae eON Tne sl eRoOnG. 

Ne7OR MAKE eae Ee OF THE CURVES. 

HOUSTY SAVED. 


NUMBER BETWEEN 1 AND 5. 


rt ht 


POMEL TT NG 
PogOnee DD A, CURVE ONE THE SCURRENT PLOT. 
CM CURRENT PLOT. 


S (2h oR. 
wars Sep wiir 
(QEPAULT, Is 8.5 X 6.0) 


R TG 
OF THE LEGEND. 


Ce ee ee ee ee en ee ee | 
THQQNAANMMYOO 
=a 


BER BETWEEN 1 AND 10. 


Oo) 
3 


CN WHICH CURVE DC YOU WANT TO CHANGE THE Y-SCALE? 
EMER CURVESNUMBER= i, 2, 3, OR 4 


WHAT IS THE NEW Y-MIN VALUE AT THE ORIGIN? 
WHAT IS THE NEW Y-MAX VALUE? 


HE GRAPH EDIT MENU 


Peer OR ADDPASCURVE ON THE CURRENT PLOT. 
V aoe CURRENT PLOT. 


hIG? 


ese ee 68te6¢ @e @ 
HANINAAAMIBAOON 
=~ OM cre Onin 
PA roe Soe Be Dr Sy fg Fy Fd 
Soe eS ariinij’) 


HF MPH OH I ee 


GCG HOAMOaAaM 
aM MOHRH Ero 
rH HSN bhtiodtH 


EN 1 AND 10. 


= 


MBER EE 


E FCLLOWING OPTIONS ARE AVAILABLE. 


W GRAPH OF OTHER CONTROLS, STATES, OR ESTIMATES. 
LY SAVED GRAPH DATA. 
T GRAPH. 

AEH ON THE TEK618. 

E METAFILE OF THE CURVES. 


tho Ato 
GAhrOoriht om 
HOM 
FHWTCH 
NWAW 
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FREVIOUSLY SAVED. 








SELECT A NUMBER EETWEEN 1 AND 5. 
4 
TRE FOLLOWING OPTIONS ARE AVAILABLE. 
1. BEGIN NEW GRAPH OF OTHER CONDTEGIS, Stiwies, Of Sto 
2e« EEELCT EREVIOUSLY SAVED GRAPH DATA. 
3. EDIT THE CURRENT GRAPH. 
4G. ELCT REVISED GRAFH ON THE TEK618. 
>. CUIT AND/OR MAKE METAPI LESSOR Tae ecu ee. 
PREVIOUSLY SAVED. 
SELECT A NUMBER BETWEEN 1 AND 5. 
5 
DO YOU WANT TO SAVE THE CURRENT GRAPH DATA [IO 
EF USED LATER TO GENERATE A YUeIA PILE? 
Y OR N 
NOTE: A METAFILE IS See aee FOR SMOOTH VERSATEC PLOTS. 
THERE WILL BE AN CPPORTUNITY TO GENERATE A METAFILE 
JOST BEFORE EXITING THIS FROGERAM. 
filteron 
YOUR ANSWER MUST BE "YY" OR "N", 
DO YOU WANT TO SAVE THE CURRENT GRAPH DATA IO 
EF USED LATER TO GENERATE A METAFILE? 
Y OR N 
NOTES Ae META iS Sine FOR SMOOTH VEESATEC PLOTS. 
TEERE WILL BE AN CPPORTUNITY TO GENERATE A METAFILE 
JUST BEFOR& EXITING THIS PROGRAM. 
y 
WHAT FILE NAME IC YOU WANT THE CURVE DATA SITORED UNDER? 
(8 CHARACTERS MAX) 
filteron 
THE CURVE CATA IS BEING FILED UNDER FILTERON DATA 
END CF DISSFLA 9.0 -- 43644 VECTORS GENERATED IN 3 PLOT FEAKES 
FRCPRIETARY SOFTWARE PRODUG? OF 2SSCO, cA eoateGe, Cag 
8493 VIRTUAL STORAGE KREFEPENCES>s 9 READS; 0 WRITES. 
THE FOLLOWING OPTIONS AKE AVAILABLE: 
cles MAKE METAFILE OF PREVIOUSLY SAVED CURVE. 
Ze. OU Tie 
- ENTER 1 OR 2 
| 
WHAT FILE NAME IS THE DATA STORED UNDER? 
filteron 
THE CURVE DATA IS BEING LOADED FROM PILE FILTERON DATA 
>> JUSING A PRE-~ALLOCATED DATASET FOR UNIT FI18FO001. 


a 


TEE FOLLOWING OPTIONS ARE AVAILABLE: 
Le OTE Deere REV rOUSEY SAVED CURVE. 


ta ee Ohare 


2 

END OF DISSFLA 9.0 -- 14919 VECTORS GENERATED IN 1 PLOT FRAMES 
PROPRIETARY SOFTWARE PRODUCT OF ISSCO, SAN DIEGO, CA. 

267% VIRTUAL STORAGE REFERENCES: 9 READS; O WRITES. 

DAsvy WZ) DETACHED 

PeoDeb22 DETACHED 

DASD 120 DETACHED 


Be eH i ee ee ee eek eae ae ee oc ok ke a oe oe a oi ck oo ae ee oe oe fk oe ok ake ok ake ok ok ke ok ok ok ok ok 


DO YOU WANT A VRSTEC PLOTTER SMOOTH COPY OF THE 
@ek DISSPLA asa Orme JUST SCR EATERS 


Me Re Ae ee ee Re oe ee ee ee ee oe ok oe oe he a eof eo oe oe oe ok oe oo ok oe oe a ek eo i a aK i ook ote ok ok ok 


Ten k/C 

Pao eOO1 LINKED R/O; R/W BY MVS 

ZL ean aa oS 

Peoomeo! DETACHED 

Soo tinG NEW FILE: 

SrmalinG NEW FILE: 

Peer res 6910 TO MVS COPY 001 NOHOLD 
Paso) i120 DETACHED 


Me He oe ek ee ee oe ae eo ee ae ee ak oe ok oe oe oe a ee oe oe oe oe oe oe eo ee oe ea oe ee oe ee oe ok ok ok oe ok ok ok ok ok ok 


Memon nee (5) "CAN BE PICKED UP AT THE COMPUTER CENTER. 


THE GRAPH(S) WILL BE ADDRESSED TO "POP (USER ID)". 
Se oe co oe oe oo oe oe eo oe oe coe oe ca ke oko ok i ko oo i ke a oo ok oie ak co oko ako ok ok a ke shake oc oe ok 


AK oe eK ie oe ie ok ie ok oe ok ee oe oe ee oe ok oie oe ee oe oe ae oe ee oe oe oe ee ee ee ee ee i ke oe ee a ok ak 


DO YOU WANT TO 


tem WOR hSYSx AGAIN 
Zen uNebee ELOT PROGRAM USING THE SAME MATRICES? 
3 Cart ANODAGR EYPE OF VoYSiIEN (OPEN /CLOSED)) 


PONE Rep Oe 
co oii ick gok ok ck igiokgok ok digi ak dk dk dick iad dak i keke ak ke ie ak aka aki oc akc ok 


3 
SG GIGI ak go Gio igigioiok ddkok kok lokoigdoigiok dokukadok dik dk dado ak dak dake dk ok ak 


HAVE A GOOD DAY!! 


Re ke oe ok ok eke ok oe te ae fe oc oe fe oe oe ok oe fie oe ae ae eo oe oe oe oo ok oe ee ot oe oe ae ok ek ae ee ok ye Ko ee ok ok a ok 


Peep Woe o1/ 2. 24 21225209 
meecora oft 
END RECCFDING OF TERFINAL SESSION 
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The graphical output generated by this example follows 
as figure —-o. 
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D. FILTER AND REGULATOR CLOSED LOOP TIME RESPONSE 


The following filter and regulator example was’ taken 
from [Ref. 7 pp. 382 - 384]. 

The full terminal session is recorded below, with user 
input at the left margin in lower case letters or numbers 
below each "?", 


ad on 
RECORDING OF TERMINAL SESSION 
Op0n70. 02 09205 340 


tee tA @® 
NRHHO 


Cor 
EGIN 
tsy 


0 toh 


THE OFTSYS EXEC CONTROLS A TREO CO! WecrAna: 


OPTsY¥ok POR TRAN Sisfiu Cetin 
PTCALC FCRTRAN CALCULATE TIME eo te 
PTPLOT ECRIRAN DISSPLA PLOVTING ROUTINE) 


AM PASSES INFORMATION TO THE NEXT PeRee ran 
DATA FILE WRIITEN TO THE USERS DISK.) tie 
E FILES ARE "OPIMAT DATA" AND “OPTPLOl para 
F THESE FILES VARY WITH THE SYSTEM ORpe yee 
OUT 20% OF THE USERS DISK SPACES) Gee ees 

T SUFFICIENT DISK SPACE I50AVAELAPIL ES 


- TYPE "EW TO EXIT, ANY OTHER ENTRY TOSCO a. 


FIND 


O 
0 
‘GR 
A 
ES 
O 
AB 
HA 


iMOrm 
man & 
THN Oy Whe 


tt Se Tt fe rt 
Na@tynwdon 


THOAHOHA 
co 


MER I I eK ee ak ak oe ea ee a ak ek i oe ek ae ee ee oe ok ee oie ae ee i ea Ke i OK Ke KK Ke ok 


YOU HAVE A DATA FILE NAMED 'OPTMAT DATA ON YOUR 
DISK THAT WAS PREVIOUSLY GENERATED BY THE ~OPTS 
FROGRKAM AND CCNTAINS THE F, G, H, GAMMA, A AND 
MATRICES FROM THAT RUN. 


A 
Ys 
B 


IF YCU WOULD LIKE TO USE THESE SAME TATR IGE SI ECR 
THIS RUN, THE CPISYS PROGRAM WSL Shear 
TESIRED DATA AT THE APPROPRIALT ES Lone 


a 


IF YOU TYPE (Y) ES. 


ANY OTHER INPUT WILL RESULT IN THAT FILE BEEING BRASED: 


coo ok cok ok doko kok dk kk kkk dokkok dokokkk kook kok kokokokokokok dofok ok kkk KK RK kK KK 
y 


eK Xe ek a ie ie ak ek ee ok oe ok ie ie ak ek eae oe oe oe i a ee eo ie a OR eo Ke ek ok ok 


x 


DO YOU WANT THE NUMERICAL OUTPUYT FROM OPTSYSX TO GO 
10 YOUR TERMINAL 14 BED} ORZ TG SA D(@SKie Veil re 


We eK Re a eae ae ee oe ae ie a a eo a eK oe ae ee eK ee ok KK KK IK ok 


S 


eH oe oR ke oR ok ke oe oo eR a a KE Ke SK KK a Ko ka Ke Ke 
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OU TEU sete eee womeOoer OUR TERTGNAL SCREEN. 
Fe RR a kk ka ok gk kk a kk kk kok kok kok kk kk ak koe ak aka kek ak akc ok 


LOADING OPTSYS 


FXECUTION BEGINS... 

OPTSYSX IS A COMPLETELY INTERACTIVE OPTIMAL SYSTEMS CCNTROL 
PROGRAM. IT WILL SOLVE NUMEROUS CONTROL PROBLEMS ON THE 
FOLLOWING TYPES CF SYSTEMS CONTROL EQUATIONS: 

XDOT = [F}*X + {G}¥*U + {GAM} * (W+W9) 
MEASUREMENT EQUATION-- 
Z = {H}*X + {D}*U + V 
REGULATOR PERFORMANCE INDEX-- 
J = 1/2 * INTEGRAL (Y *{A}*Y + U *{B} *U) DT 
STATE FEEDBACK GAIN DEFINITION-- 
p= = {Cc} #X 
DO YOU WISH TO CONTINUE? ‘TYPE "YES" CR "NOM, 
: --DATA ENTRY-- 
ALTHOUGH OPTSYSX IS SPECIFICALLY DESIGNED TO READ 
ALL MATRIX DATA INTERACTIVELY, SEVERAL ALTERNATE 
METHODS ARE AVAILABLE TO USERS: 
METEOD 1--THE "F"oecnt AND "GAMMA" MATRICES 
MAY BE READ FROM SEPAEATE DATA FILES. 
METHOD 2--THE "F","G") AND "GAMMA" MATRICES MAY BE 
EXPLICITLY DEF{NED WITHIN SUBROUTINE "SETUP". 
{NOTE: IN EITEER CASE, THE USER SHOULD OBTAIN A COPY 
OF THE PROGRAM LISTING AND EXAMINE 
THE EXAMPLES CONTAINED IN sit "SETUP" 3 
DO YOU WISH TO CONTINUE? TYPE "YES" OR "NO", 
y 
DO YOU WISH TO INPUT THE "F", "G", AND "GAMYA" 
MATRICES FECM SUBROUTINE "SETUP" IAW THE 
METHOD DESCRIBED CN THE PREVIOUS SCREEN? 
TYPE "YES" OR "NOt, 
nN 
GENERAL OPTSYSX OPTIONS: 
OFTION 1 -- SYSTEM ANALYSIS WITHOUT 
OEEN-LOOP EIGENSYSTEM CALCULATIONS. 
CETION 2 -- SYSTEM ANALYSIS WITH OPEN-LOOP 
EIGENSYSTEM CALCULATIONS. 
CETION 3 —- OPEN-LOOP EIGENSYSTEM FOUND 
AND PROGRAM TERMINATES. 
("FUMATRIX ENTRY FOLLOWS IMMEDIATELY. } 
OFTICN 4 —— MCIAL DISTRIBUTION MATRICES COMPUTED 
WITHOUT FILTER OP REGULATOR SYNTHESIS 
OR STEADY-STATE ANALYSIS. 


=) *s) 


Shae T eae OPEL ON: 


1,2, 3, 


OR 4. 


DO YOU DESIRE RMS VALUES OF STATE eeoen ce 
A YES " 


yee 


OFTION 1 


OFTION, 2 


OFTION 3 


OFTION 4 


OR TNO 


CPISYSK LOR7ZGEA SS CROP riOns : 


j 
‘ 


OPTIN ALeRiL iar 
SYNTHESIS Vir 
MATRIX INPUT. 


OPTIMAL OE Lia 
OLN TH So 5 nee 
MATRIX INPUT. 
OP TOM Ad Pale ea 
SYNTHS iS nee 
MATRIX INPUT. 
OPTIMAL FIR 
SYNTHESIS Wer 
MATRIX INPUT. 


AND/OR REGULATOR 
NO BATERNAL “CY OR ag 


ae REGULATOR 
GALERNAL "CU 


m™ 


-_EGULATOR 


ten AND WKN 


PO YOU WISH TO LTETERMNINE JHE STEAD Y= 1A wero oP ONS 
FOR A CONSTANT DISTURBANCE? 


" YES " 


DO. YOU WaSH TO, LETERMENE TH ESM@DAL. DiS Tare non 
AND GAIN MATRICES? 


" YES " 


tyes 


ee 


CR aN OU. 


Ch “NOt: 


OPEN-LCCP TRANSFER FUNCTEOCN SOPs: 
NO OPEN=LOOP TRANSFER PUNCTIONS COMNPUDER 


GETION 7 


CFTION 2 
CFTION 3 
CETION 4 


SELECT AN OFZGOn: 


POrtoon 


RESIDUES, 


AND ZEROS COMPUTED. 


—-~- ONDY POLESMAND ZER@S CONPULED: 
== ONLY PORES Ay eRe UES eon Pure 


1, 2, 3, 


NOISE TRANSFER FONCTIIONSO Lawes. 
NOISE MTR ANSrae FUNCTIONS CoumPuLlEepe 


OET ICN aa 
OFTION 2 
OETLONE 
OFTICN 4 


NC 
Ee 


LES, RESIDUES, 


OR 4. 


AND ZEROS COMPUTED. 


ONLY POLES AND ZEROS COMPUTED. 
ONITY POLES AND RESIDUES COMPUTED. 
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eeecl AN OPTICNGgel,2, 3, OR &. 


? 
1 
COMEENSATCR TRANSFER FUNCTION OPTIONS: 
CFTICN 1 — NO COMP. TRANSFER FUNCTIONS COMPUTED. 
CETICN 2 — PCLES, RESIDUES, AND Z=ROS COMPUTED. 
OFTICN 3 -- ONLY POLES AND ZEROS COMPUTED. 
CETICN 4 — ONIY POLES AND RESIDUES COMPUTED. 
(NOTE: A COMPENSATOR TRANSFER FUNCTION CAN EE 
CCMPUTED ONLY IF BOTH A REGULATOR 
AND FILTER ARE SYNTHESIZED 
ANL/OR INPUT. } 
SuepeT ANG OETION@ 4, 2, 3, OR 4. 
? 
1 
WILL A FEED-FORWARD DISTRIBUTION MATRIX 
{"D" — MATRIX} EEF INPUT ? 
TYPE "YES" OR "NO", 
Nn 
THIS OFTICN DETERMINES THE CRITERIA FOR DECIDING WHEN A 
MARKOV PARAMETER IS ZERO-THE MARKOV PARAMETER INDICATES 
THE ORDER OF THE NUMERATOR POLYNOMIAL OF EACH TRANSFER 
ALL 'N" ZEROS OF THIS POLYNOMIAL ARE PRINTED OUT AND 
THIS TEST TELLS HOW MANY EXTRA ROOTS EXIST AT Z = 0. 
LESS THAN 10.0**{-IF} IS CONSIDERED ZERO. 
THE DEFAULT VALUE OF THIS PARAMETER {IE} IS 6. 
IN CTHEE WORDS, IE = 1.0E-6. 
IF YOU DESIRE A DIFFERENT MARKOV CRITERIA, 
TYPE THE INTEGER VALUE. 
, IF YCU DESIRE IHE DEFAULT VALUE, TYPE "0" {ZERO} 
0 


MomreeueNFolks TO SYNTHESIZE A STASLE FILTER {OR KEGULATOR} BY 
Peon etbestNG THESORIGINAL SYSTEM? 


{NCTE:WORKS FCR FILTER OR KEGULATOR BUT NOT FOR BOTH 
IN THE SAME RUN.} 


eer tes YE oes RNG tee 


DO YOU DESIRE TC PRINT THE EULER-LAGRANGE EIGENSYS‘TEM 
FRICR TO DECOMPOSITION {FOR CHECKING THE PROGRAM}? 


Pere YESS CR NOT. 


® 


MOWER SPECTRAL DENSITY {PSD} OPTION 7 ; 
Orimete i —— COMPUTE THE PSD OF THE OUTPUTS AND/OR THE 
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CONTROLS OF THE CONTROLITEDSS2 STEW vit N PORCH pe 
PROCESS AND MEASUREMENT NOISE. NOTE: BOTH A 
REGULATOR AND A FILTER MUST BE RESIDENT IN THE 
PROGRAM IC US& THIS OPTierey 
OPTION 2 -- SAME AS CETION 1 ABOVE BUT ONLY PRINT THE 
RESIDUES CF EACH TRANSFER FUNCTION 
USED IN THE PSD COMPU re: 
CPTITCN 3° =3 NOt DEST EE 
z SELECT AN CETEON: 1p. 2 pe@ibe oe 
3 
CO YOU DESTRE REGULATOR Siewert Ss ONG? 
IYPE “YEou CR Maier. 
Nn 
THE Ju Ne ue iG ne “GANAL eA it AND W Et YATRICES 


ErOl oy OUR PREV TOUS OPTSYS RUN WERE SAVED. 


THE FOLLOWING OPTIONS ARE AVAILABLE: 
1. USE ALE OF Tht eon teks AS AGN. 
Ze USE SELECTED (helo hc 
3. ENEUI ALL NEW OMATRT Crs 


ENDER tc po cee 


NOTE; EACH SAVED MATRIX WiLL BE REDISPea ea! 
THE PROPER INPUT SEQUENCE INTERVAL 
AND YOU WIIL HAVE THE OPTION OF CHANGING 


INDIVIDUAL MATRIX ELEMENTS. 


FLIAG/PARAMETER SETIINGS BORP THIS Rune RE ioe een. 
TOL I¢ IR 2SS GM TTR ears Bee De eee 


0 0 arauO 0 0 0 0 0 0 “0 
ISET IDSTAE IPSD IYU INORM IKEG NS NC NOB NG 
0 0 0 0 0 0 7 1 
ORDER OF SYSTEM = 2 : 
NUMBER CF CONTROLS = 1 
NUMBER CF CBSERVATICNS = 1 
NUMEER CF FROCESS NOISE SOURCES = 1 
THE SYSTEM MATEIX ["F"-MATRIX}... 
G84 88888 


TO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
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mie YES On meiOwt. 


nN 
eee e ORD Y NAMI CS MAG di Xceemeteterescie ele eels estes Fes 
0.0 Q0.1000D+01 
0.9 -0.4600L0+01 
THE MEASUREMENT SCALING MATRIX {"H"—MATRIX}... 
1. 90000 0.0 
DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
iveeiee LEO" CR UNO, 
n 


(eon w SCALING MAT RIX sc 00.0. c osc cees cle 
0. 10COL+01 OF0 


Bcc OND ROL DESIR IB UIrONeraAika Xx {"GUSKMATRIX} . .< 
0.0 
0.78700 


DC YOU WISH TO CEANGE THE VALUE OF ANY MATRIX ELEMENT? 
(ive ieee (CR "NOM. 


N 
ENTER THE FEEDBACK GAIN CONTROL MATRIX {"C"-MATRIX}. 
DIMENSION = # CONTROLS {NC} X # STATES {NS}. 
THE ELEMENT C( 1, 1)= 
=-254.1 
: THE ELEMENT C( 1, 2)= 
-19.57 
THE FEEDBACK GAIN CONTROL MATRIX {"C"-MATRIX} 
-254.1000C ~-19.57000 
DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
: TYPE "YES" OR "NO", 


eemechinOl (Dr SARTBUTTON MATRIikss cece ecce o Gan 
0.0 
O27 870D+00 


THE OPTIMAL FEEDBACK GAIN CONTROL MATRIX...C=BINV*GT*S... 
eo TODtO2 ~=1.95 705401 
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THE CLOSED LOOP DYNAMICS Mal Riera te, 
0.0 1.000000D+00 
—-1.9997670+4 025-77 0.00 eo 04 
C-1OCP SUBOPI. REG. E-VALURPSSe te - tG*Ce, 
—1,000CED40 1, 9. S93 CAre0C: 


C-LOCP RIGHT EIGENVECTOR MATRIX.32 2... .55Ms5.. 
-5. 00098 0p - 025 - FI 3 96GU2D U2 
1.000000D+CO0 0.0 


CCNTFCL EIGENVECTCR MAT Ritter rete tere ye eteeee ere 
-6.€62510D+00 1. 27033 3p 20a 
C-ICCP SUBOPI-REG. LEFT E-VECTOR OAT Ro hee Thy 
0.0 1.000000D+00 
-2.00015SD+C1 —-llGG0 2 Gnas 


THE PROCESS NOISE DISTRIBUTION MATRIX 
{NGAMMA"=MATRIX}... 


0.0 
Ueno ay 


DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEM Eee 
TYPE VES “Chern: 


N 
ENTER THE PROCESS NOISE PSD WEIGHTING MATRIX 
{"C"MATRIX}. 
DIMENSION = # PRCCESS NOISE SOURCES ie X 
#PROCESS NOISE SOURCES {NG 
‘ THEeELEMENDT 0 ( ied j= 
10 
THE PROCESS NOISE WEIGHTING MATRIX....<Qe- 
10.00000 
DOC YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
TYPE "YES" CR "NO", 
nN 


FROCESS NOISE DISTRIBUTION MATRIX.....GAMMA.. 


On9 
0. 10C0D+00 
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Eemeke SPECTRAL DENSITY e-=PROCESS NOTSE. ...Qa. 
0.10C00LD+02 
iio ec PASUREMENT ROUSE DISTRIBUTION MATRIX {"R"MAT RIX}. 


DIMENSION = # OEFSERVATIONS {NO} X # OBSERVATIONS {NO} 
THE ELEMENT R( 1, 1)= 
6. 000000 1 


fPteeood RE MENT NOISE DISTRIBUTION MATRIX»... ...Re< 
C.900000 


POmwoURW2 oh) 1O CHANGE THE VALUE OF ANY MATRIX ELEMENT? 


Ce Eee YS. SC Re NO tee 


Nh 
POWER SPECTRAL DENSITY-MEASUREMENT NOISE..R.. 
0. 10COD-06 
FNIER THE FEEDBACK GAIN ESTIMATOR MATRIX {"K"-MATRIX}. 
DIMENSION = # STATES {NS} X # OBSERVATIONS {NO}. 
THE ELEMENT K( 1, 1)= 
95.4 
: THE ELEMENT K( 2, 1)= 
456 1 
THE FEEDBACK GAIN ESTIMATOR MATRIX {"K"-MATRIX} 
4567200000 
DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
; TYPE "YES" CR "NO", 


FIILZER STEADY STATE GAINS. ccc ce eee ce cee een Kees 


S.54000Q0D+01 
woo ro CUDtT O03 


Gl@sSEDSLCOP FELTER DYNAMICS MATEEXR  I5.2.. 


D+#+01 1.000000D+00 


Eee 
=o oe OOO 
561000D+03 -4.600000D+00 


ie 


EIGENSYSTEM OF OFTIMAL ESTIMATOR. ....e-e- 


Groce SUSOP Ta ot S- VALUES. -~ .DET(Si-Ftk). . 
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—5. OOUCOT Oat ees io cen rs 


C-LCOF RIGHT BSIGENVECTOR MATRIX. Oe eM... 
9.953595 /D= Cem 10 Joop 
1.00000CD+00 0.0 


MEASUREMENT EIGENVECTOR MATRIX.....H (EAR) *M.. 
9299595 70-03 ==". 0 oop 0 2 


C-ICCP SUBOPT. FILTER LEFT BeVieren Bie een oe 
Ee 1.0000C0D+00 
=9. 12229207 0V “92060 Zoo 


THE COVARIANCE OF THE ESTEMATUONe Sr MO teers one 
© 150503D-06" 2527 eer 
5: Z2INQJVQOD=Ca™ Tete rep ez 


RMS VALUES OF THE ESTIM AT ICN SEG eee < oseeereere 
2.674042ZD-03 I.0868Gap-—01 


DO YOU WISH TO OBTAIN A TINE VRESeONSE 
OF THE SYSTEM YOU ARE EVALUATING? 


(Y OR N) 


NCTE: YOU MUST BE LOGGED (ON AT AS DUAL SACRE 
(TEK 618) TERMINAL ITO UTILIZE THiS @iepE. 
TE ee he ee at G (COMES ae COCR 3AM a OISE), 
A (QUTEUT See AN B_ (CONTE L_ COST) MAT fCES WILL BE 
SA WE eOer KF ENTRY TQ THE MAIN OPTSYS PROGRAM. 


He Fe eo A KK ea a Ka oe He ee oe ee ke eae ee eae Ke eK a ok ok ke oe a eK eke ee aK ie ok ok a ek ok KK 


IF YOU ARE DISSATISFIED Viti ete reo LS 
THUS FAR AND WCUDD LIKE” 20 Gxt) so] elo, 


=T YP.En * Y Ses TORE eld = 


ANY OTHER INPUT TO CONTINUE 
eK KK 


He MK I He AK te ee I IK ok oe ee a ee kK ok Le, BG cna eee 


OPT CALC... 
ON BEGINS... 


bit 
~O 

by 

OW +# 
CH 4 
HI 
Ha + 


- SELECT THE TYPE CF SYS) Sher hoeens lh Cm. 
(CrEN LCOP, CLOSEL LOOP TOE she ER? ne oeraenon 
- PROVIDE START AND STOP TIME FOR PLOTTING CAL Be aan 
- SELECT THE TYPE Creep R ie Creu fC) ee Hoe eae By 
- PROVIDE START AND STOr EINES yeCky Oita em PUNCTIO (S) 
- FROVIDE DRIVING FUNCTION MAGNITUD4®(S). 


CLEARB@EEeSsenREEN TO CONTINUE 


THE F MATRIX 
ee 1.00000 
0.0 -4.60000 
THE G MATRIX 
O20 
ve 7 07 00 
THE C MATRIX 
=204, 10000 migs.5 70CH 
THE H MATRIX 
1.900000 0.0 
DE OR Ao ek 
$5.40000 
a One O00 0 


meee OLuLOn TING PLOIDING CPTTONS BARES AVAIVABLE IF THE 
PaO Pee bei AT RICE S@HERE CALCULATED @IN GCPTS YSX; 


1.  OFEN LOOP TIME RESPONSE 
XDOT = {F}*X + {G}*UC 
2. CLOSED LOOP TIME RESPONSE 
RDOT = {F-G*C}*X + {G}*UC, U= {C} *x 
3. OFTINIZED FILTER CLOSED LOOP SYSTEM RESPONSE. 
XDOT (HL +E TP ioNs UC, Z, H} *X 
XEDCT = } *XB Shy 0 + {K}*{Z - H*®XH} 
FILTER + REGULATOR CLOSED LOOP SYSTEM RESPONSE. 
XDOT =_(F+GeC}#X + {G) * Z = ta} 3% 
XHDCT = {F}*X $ {6} NEL * (2 - H*¥XH}, U = {C}*XH 
SELECT 1, 2, 3 OR 4. 
THE (G*C) MATRIX 
0.0 0 
-199.67670 -15.40159 , 
THE (K*H) MATRIX 
$5.40000 O10 
4561: C0000 : 
THE COMBINED SYSTEM F MATRIX (2*NS X 2*NS) 
0.0 1.00000 0.0 0.0 
0:0 -4260000 -199:97670 -15.40159 


oi 


500 


5.-40000 0.0 =95-40000 1.00000 
tT. C0000 0.0 =a 700s oon aenes OO a9 
THE AUGMENTED G MATRIX (2*NS X NC) 

0. 0 
0.78700 
Cad 
Cee 0.0 

Al WHAT TIME DO YOU WAN?) TOs Se27 

THE TIME RESPONSE (GAECULATIGCNE 

INPUT START TIME IN SECONDS. (NORMALLY 0.0) 

AT WHAT TIME DO YOU WANT TO STOP 

TEE TIME RESPONSE CATCUMRA ets 

INPUT SIGE TINE iN steeiys- 


THIS PROGRAM DIVIDES THE TIME INTERVAL YOU HAVE 
JUST SPECIFIED INTO UP IC SO0RSUALT Se LN i eave ee 
IHE INTEGRATION AND PLOTIING ROUTINES. IN ORDER 
TO SAVE COMPUTER TIWE, THE NUMBER OF PORNTS (ean BE 
CAN BE REDUCED WITH SOME LOSS IN CURVE FIDELITY. 


HCW MANY POINTS DO YOU WANT TO CALCULATE? 


DOES THE SYSTEM UTILIZE A DARLVING = BUNCTION SeenreoL INP Ue 


Nn 


DOES JHE SYSTEM START WHITH ALL EPNiirAL Conve 


Oed 


(Y)ES OR (N)O 


It 
© 
e 
© 

v 


(Y)ES OR {N)0? 


WHAT IS THE INITIAL CONDE vO Nr Orme 1) eee 
WHAT IS THE INDTDAL SCOR DET LONSY OR coal (ae 
WHAT IS THE INITIAL CONDITION FOR zie 
WHAT IS THE INITEAL CONDI TGON 1 OR ane een 2 


THiS IS YOUR VLAST Reel On itm Tyee 
MAKE CHANGES IN THE FOLLOWING AKEAS. 


ie SELECT ANOQTHER PYVREwOr SiS TE tem ro: 


Bis 


FILTER OR FILTER/R EGULATOR) 


CLOSED, 


(OPEN, 


DL AR Tease Ore NES 
DRIVING FUNCTIONS 


2 


oe 


INITIAL CCNDITIONS 


CONTINUE 


4. 


Se 
SELECT A NUMBER BETWEEN 1 AND 5. 


CATICN OF PROPER PROGRAM OPERATION. 
DIATE AND Bolla thes TIMATES CA N@EE, PLOTTED. 


L 
I 


HE FOLLOWING INFORMATION IS PROVIDED ONLY 
Ck AN INDI 
AT ieCCNIROLS, 


Ei hy 


X (3) 


X (2) 


X (1) 


U(1) 


TIE, 


TOT MMT TANYA NAAQS SCT TSM TK TK ONNNAANNAAANNAMMAMMNMAMNMAMMAMMMAMMNN 
DODDOVD CDODOVOV OC DOWDVD OVO CODD TODVDOOVDODVODO VOODOO VOVOOOCOO 
t{+tretlrtitereerett_|toprtetretrttt'etettbtrbrerte)!et¢t)rtrytrereetas 
LaanunaanaganaaanaeaaqgaAaaannAanAaanaAaaaa aa aamaaoaaomaasd 
DTI MNITAADA™ OW THN W ODOM N STH NHNVOANAOKYNMNOOO TAK TOrHNON 
DO TNO DOMNHODNT KHDDOMME KKH OMIT KH HOO FONODONOFTOMRMETNMNOOMm mmr ON 
MOF OND MOK ST ANT OMNAHDOCO ANTOWAN ANTNKNCODNOODO—-DTOMW 
ONDA SOF FON FS TMHFRMYHADNNANS WIS DOM NMORMLM FANDOM DEHDOHANAM 
OFM NONINES — DNINAT DNOMOMOMEM OMNOMDAINTADONMEMODOKONDNArMWIN 
OMT DOR NNVN KH ANgerNOWOO NANA TFONMNATSKMNNMMNMNAMNONAWMNOrP—K—NCwoO Ass 
Be Te On Oo aa, a cen 
Cr ee er er ee eee2e 2«§ &©& ee 8 @ @ @ @ 
SSCS COCO OOOO OO OOOO OO OOOO OOOO OC OOOC OOOO OOOO OSDCOOSOK 
{ tfeutetetoue tfteetet_tet’ertetet 
DODDODOCOSO DODO DDDD DODO CODCOD CDODDDDDOOVOCODVDVODVOCOOCOOCOCOCOCOOCO 
++r++eteteee eee eee i_ UC trrrrreeee de erretderedeett tbe!e td etetvbptert 
AARAARAAARAARAARAARAANRAARAARAARAAAARARARAAAAAAANRAARARARAARAARAARNA 
MONTNMNAMNONTODOHNDOOMMNOMNTOVWODAANST SFAPFOODOAMNVHAMANMNAMmomMmM=st st 
TNT DMNDNOMNODM FNDOOMN OM KH OF DODMNHONTHDMNHMNOMNE™ SMOMKe MNO HOO 
TNO MNAFMNM FAM HDDMNDANOTKHKMNONAN FHWODNSTHRHTHADNAOM ANWDHANOWOKMODONM 
D™MNNOMNEAQNDOS FOCSCANDKHANOOS AHTMOM KK NHODHANMNANIM™AMADODO 
WTAMNWOM OPK ANNO ADNAN ANANT THOM STVDANHOOHNODMNHMNOATPTONNOCAOOMGA 
DAMN YANO VWOSFTEKNAST FTOMNDVODSOCO OT FOPANK—MNOOARHDRTOOT™MNoPYwors 
DO FDNT THT TOON ANAT OMN HK FOODHDNSCTKKNNOM OM AMNMNACSCEKR SOT rnnstst 
Cr er er 
SMG a LET ee oe ae Manin, 
tfrrrtetrertets 
DOH MONA NGS ST TTS TT THO NANAANANAANANANANMMOAMMAMNMANMNMAMMMNMAMNEH 
DO DO DMO DODOVSN DOOD DODD OD ODOC D CODD OVDOVOVOOVDVOCVDCCOCOOC0O0O 
+frtrtirrertirttetrtrrt)ererrer)_tverrr+)rtteterertttrerrttetretetrerrgs 
AAARAAARARAANRNAANRANAARAANAANAANAANAAANAANAARAANAAARAANRAANAAANAAANAARAANAAANAAAQAANQAN 
DTANDHAMNOWN ST FONDS VDT DONO +ONOM ANEAMNMNANNONSTMNODODOSATTPOKRVHON 
ODWINWIIADNT ONNAAM DOMAMOM DOTMOMNONNOHODNDOAOREKNMNOO? Hl ON 
DOTmONANADNSAATANO KODMNKH CODON K-HNOOMNMNADODAMNIFN KH MNKNDODOWDNODDM—OOMNW 
OVDNANANOOlM OTFOME OM OD FOMAANATNET DOM AMON FHANHOFOMHOANSAM 
OTF OTFHNT(NQDOMAOAN ST DNWMOM OME MNONOMNDWNNIMTNONEKITMODOCKTONNAMWLN 
ODO DY IDNA MHOWOON FIND YON FMNATWDNMNNAT SK MNDMOMNODHNOTMPrNCweO AOS 
TNO NFAT OF PDS HT TMK TT THT TH NOM ON ANAT COAT KHMNNOOPE EM OWMWN 
re er ee | 
ae ee rT ae Games TT | TPT TMi ee 
i t | \ i 


DODO DOD DODO DOOD DCO DODVDOV DD DVDODDDVDDDODOVOV CDDODODVDVCOCOCOO 
e* @ 8 e*ee?ee ®e eee @ @ @ @ eeee#ee?e##e#e8t eee @ 


—olelololelololelelslelololololelolereio(aiclololololorslelolovororo) DODOCODTOOCCO0O 


Mm ATWO WOW NMAFAOlP DNONAM SFUIVWVDNOPeNTIOKM WOT NAM STOMP WHAODNMAWNO 
DOOO OCOD CO SCT K TTT TTHANANNAAANANAMMMMNMAMNMN MY = sa as Toss aT AT NMIUNMWMW 
i rr er ee | 


CO DODD DODOODOOOD DODOVDOVD ODDO DODD OOO DOSS DOOCVCVODOOCO0O 
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0.58 0.0 Q.4904551D-03-0.4634245D-02 0.4904551D-03 
Uo 0.0 0.4347377D-03-0.462959 (pe 0 4 347s 
0.60 0.0 0.3798064D-03-0. 4aUerech- v2 053 798050n-03 


eR HE Fe KF ee Fe Fe ee ae eae he oe ie ee a eI OK a 2 OK oe Kok ok a ea a ke oe oc ok a a oe a ek Ok OK 2K 


Ir YOU ARE DISSATISFIED eo eee eouLTs 
THUS FAR AND WOULD LIKE TO EXIT TO CMS, 


=TYPE *Y* 7O Exai- 


(ANY OTHER INPUT TO CONTINUE) 
CCC CCSCLCSSE CSCS TCS SS LSS CLS TCLS CCC CCS Se SSS CT e Ce ee ee ee ee se se eS 


B (120) E/0O 
C 1321 R/O 
EM 122) “B7© 


-es. Your Fortran programers nor ber1ng wi oaded ee 
ee. Execution will soon follow ... 
EXECUIZON BEGINS. 


Thao er ORE i ae a FROGRAM PLOTS: 


= FRELEACK CCNT ECL INPUTS, 

= STATE ESTE ie omeaD 

= RECCNSTRUGC TION ERAORS 
FROM THE DAIA THAT Yours! CALCULATED. 
ITFE CAPABILITY IS ALSO AVAILABLE TO REVIEW ANY 
GRAPHS THAT YOU HAD PREVIOUSLY SAVED AS DATA 


FILES CN YOUR Disk. 
CLEAR THE SCREEN TO CONTINUE. 
THE FOLLOWING OPTIONS ARE AVAILABLE: 


1. PLOT THE DATAVYOUR2UST pear Ccuree 
2. PLOT A CURVE THAT YOU PREVIOUSL 


ENCLE eel Orn ae 


ED. 
¥ SAVED: 


YCU MAY PLOT UP TO SSySTem VAREABLES@VS TEME- 
HOW MANY VARIABLES DO YOU WISH TO PLOT? 


3 


WHICH TYEE CF VARIABLE DO YOUUWHSH TOSPLOT AS CURVE NUMER 


1. STATE VARIABLE (IE; X1 X2y ETC) 

2: FFEDBACK CCNTRC f, 4 6G = 

3. CONTROL INPUT (ig, 1 Woe Bme) -) 

4. STATE ESTIMATE a véR) (IE., XHA aT XHAT2, ETC.) 

5. STATE BE CONS TEUCTION ERROR (IE., X1-XHAT1, 
X2-XHAT2, ETC) 


ENTER 1,2,3,4 OR 5 
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met Vis thes tpoenie tor stot STATE VAST ABLE THAT 
YCU WANT TC PLOT AS THE NUMBER 1 CURVE VS TIMMS? 


(tio THE CURVE DABE Reon THIS VARIABLE? 


NOTE: 1. 40 CHARACTERS MAX LENGTH 
2. GREEK SYMBOLS WILL BE PRINTED FOR ANY LETTERS 
ENCEOSED IN Beet Bee 
HA 


EE. A) => A 
E) => BETA 
F) => PHI 
A C) => THETA 
angular FrOoSition = (x)1 


WHICH TYPE CF VARIABIE DO YOU WISH TO PLOT AS CURVE NUMBER 2? 
1. STATE VARIABLE (IF., X1, X2, ETC) 
2. FEFDBACK CONTRCL (If., 6 = =c*X) 
3. CCNIROL INPUT (IE., U1, U2, ETC.) 
4. STATE ESTIMATE (OBSERVER) (IE-, XHAT1, XHAI2, ETC.) 
5. STATE RECONSTRUCTION ERROR (IE., X1-XHaT1, 
X2-XHAT2, ETC) 
; ENTER 1,2,3,4 OR 5 
2 


oes UP oOe hemmed coe GALE BOTT MATE THAT 
PeUR Nant PC PLOI MAS TRE RVUUBER. 2 CURVE VS TIME? 


WHAT IS THE CURVE LABEL FOR THIS VARIABLE? 


Nett ete 40 CHA RAGELERS MAX LENGTA 
Z-On tok eo OCLonhLG BE PRINTED FOR ANY LETTERS 
PMChOSEy. IN PARENTAES ES. 
IE. {A) => ALPHA 
B) => BETA 
E> PHS 
_. Q) => THETA 
angular position estimate - (x)el 


WHICH TYEE CF VARIABLE DO YOU WISH TO PLOT AS CURVE NUMBER 3? 
1. STATE VARIABLE {IE., X1, X2, ETC) 
2. FEFDBACK CONTRCL (if. = =C*X) 
DeGONTRCL LNPUT (Tm, 01, U2, ETC.) 
4. STATE ESTIMATE (OBSERVER) (IE., XHAT1, XHAI2, ETC.) 
5. STATE RECONSTRUCTION ERROR (IE., X1-xfaT1, 
X¥2-XHAT2, ) 


ENTER 172,2,4 0R 5 


Noe 


io oun sino hie Ob Or win S FEEDBACK CONTROL THAT 
VOU iiieTOmehOT AS THE NUMBER 3 CURVE VS TIME? 
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WHAT IS THE CURVE LABEL FOR THIS eae eee 
NOTE: 1. 40 CHARACTERS MAX DENG 
2e GREEK SYMBOLS WILL BE PRENGED FORSANY La toes 
ENCLOSED IN PARENTHESES 
IE. (A) => ALPHA 
B) => BETA 
F) => PHI 
; Q) => THETA 
input voltage = 94 


YCU MAY USE UP TO SyHEaDiiGS- 
HOW MANY HEADINGS DO YOU DESIRE ON THIS GRAPH? 


OF ee na 


WHAT IS THE DESDTRED AEADING Sere 


NCTE: 1. 40 CHARACTERS MAX LENGTH 
2e GREEK SYMBOLS WELL GBE Se Re Nt De Okeer th eri Bee 
ENCLOSED IN PARENTHESES. 
IE. (A) => ALPHA 
B) => BETA 
F) => PHI 
0) =>. TH 
filter + regulator closed pecs system 
WHAT IS THE DESIRED BEADING Oiees. 2e 
NCTEs 1. 40 CHARACTERS MAX LENGTH 
2- GREEK SYMBOLS WELL Be CORINDTED FORSANY Lettie: ss 
ENCLOSED: IN EARENTHES ES. 
IE. (A) => ALFHA 
B) => BETA 
F) => PHI 
QO) => THETA 


example 5.1 
WHAT IS THE PESITRED HEADING SN Ue reo: 
NCTE. ag ae CHARACTERS MAX LENGTH 





- GREEK SYMBOLS WEEE BE BREINTED FOR ANY LET @aa 
ENCLOSEL IN PARENTHESES. 
IE. (A) => ALPHA 
B) => BETA 
EF) => Pe 
Q)_=> eps 
linear optital cConmtrel syste 
2> USING A PRE-~ALLOCATED DATASET POR LUND FLIP OoM: 


THE FOLLOWING OFTIONS ARE AVAILABLE. 


NEW GRAPH OF CTHER CCNTROLS, STATES, OR ESTIMATES. 
PREVIOUSLY SAVED GRAFH DATA. 


Un GIN 
® @ 6 


8 
PO "OD bd bo to 


Sly SAV noe 


H 
Ma ° 
ae OR MAKE METAFILE ORS tie etry eS. 
NUMBER BETWEEN 1 AND 5. 


oe 





THE GRAPH EDIT MENU 
1. CHANGE VARIABLES OR ADD A CURVE ON THE CURRENT ELOT. 
2. DELETE CURVE FROM CURRENT PLOT. 
3. ELIT CURVE TITLE(S). 
4, EDIT PAGE HEADING(S) . 
5. CEANGE THE Y-AXIS SCALE. 
6. CHANGE THE TIME AXIS SCALE. 
7. CHANGE PLOT SIZE. (DEFAULT IS 8.5 X 6.9) 
8. CHANCE THE LETTERING HEIGHT. 
9. CHANGE POSITICN OF THE LEGEND. 
10. EDITINS COMPLETE. 
SELECT A NUMBER BETWEEN 1 AND 10. 
S 
HOW MANY INCHES IN THE X DIRECTION 
LEFT OR RIGHT), DO YOU WANT TO MOVE 
OVE THE LEGEND BCX FROM ITS PRESENT POSITION 
NOTE: 1. DEFAULT PLOT SIZE IS &.5 X 6.0 
2. LEFT IS NEGATIVE 
: . RIGHT IS POSITIVE 
3 
HOW MANY INCHES IN THE Y DIRECTION 
(up OR Dow) DO YOU WANT TO MOVE 
OVE THE LEGEND BOX FROM ITS PRESENT POSITION 
NOTE: 1. DEFAULT PAGE SIZE IS §&.5 X 6.0 
~ DOWN IS NEGATIVE 
; 3, Ur PS POSITIVE 
2 
THE GRAPH EDIT MENU 
1. CHANGE VARIABLES OR ADD A CURVE ON THE CURRENT PLOT. 
2. DELETE CURVE FROM CURRENT PLOT. 
3. ELIT CURVE TITLE (S) 
4. EDIT PAGE HEADIN (3) 
S. CHANGE THE Y-AXIS SCALE. 
SeMGHANIGE GZHE TIME AXIS SCALE. 
7. CHANGE PLOT SIZE. (DEFAULT fas Gsm. 62:0) ; 
8. CHANGE THE LETIERING HEIGHT. 
9. CHANGE POSITICN OF THE LEGEND. 
10. EDITING COMPLETE. 
SELECT A NUMBER BETWEEN 1 AND 10. 
? 
ON WHICH CURVE [TO YOU WANT TO CHANGE THE Y- SCALE? 
: SNTER CURVE NUMBER- 1, 2, 3, OR 4 
1 
. WHAT IS THE NEW Y-MIN VALUE AT THE ORIGIN? 
eS 
i WHAT IS THE NEW Y-MAX VALUE? 
0.125 


THEEGhArPHe SDT MENU 
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VE SONS THE CURRENT Pict: 


OR ADD ASU 
Tt PLOT. 


FIG) 
26M Oth 


GQ) 


~~ Cc} 
WHO HHQaWdO< 
QOmry math Aaay 


g {3 ALE. 

AXIS SCALE. 
Ae ace Lor oe aus 6 Cie, 
G HEIGHT. 

F THE LEGEND. 


OwWO~IJIO01E Wh 
@e*e6see6@4e0ee@e 6 
Orin tr rOowohty 
bi > es fo OH Hit 4 
F222 242nHm=m 
Wr AKA 
SS PHM WH I eH © 


aad 
HRaADAAMNM 


SELECT NUMBER BET N 1 AND 10. 


E FOLLOWING CPTIONS ARE AVAILABLE. 


F OTHER CONTROLS, STLTES, OR ESTIMATES. 
VED GRAPH DATA. 


A 

APH 

ON THE TEK618. 

TAFILLE OF THESCURVESS 


LOWING OFTIONS ARE AVAILABLE. 


GRAPH OF . OTHER CONTROLS, STATES, OR ESTIMATES. 
VIOUSLY SAVED GRAPH DATA. 
URRENT GRAPH. 

ED GRAFH ON THE TEK618. 

R WAKE METAFILE OF THE CURVES 


Un4= GaNu 


5 
DO YOU WANT TO SAVE THE CURRENT GRAPH DATA TO 
EE USED LATER TO GENER ADE ee ee 
Y OR N 
NOTE: A METAFILE IS REQUIRED BOR S1002 TH VSRoA fe eCrc le ea. 
THERE WILL BE AN CPPORTUNITY TO GENERATE A METAFILE 
JUST BEFORE EXITING THIS PROGRAM. 
Y 
WHAT FILE NAME [LC YOU WANT THE CURVE DATA STORED UNDeae 
' (8 CHARACTERS MAX) 
filtereg 
THE CURVE [ATA IS BEENG FILED UNDER FILTEREG DATA 
END OF DISSELA 9.0 -— 26332 \VEGTORS, GENERATED IN 2 PLOT PRA 


FROPRIETARY SOFTWARE PRODUCT OF ISSCO, SAN DIEGO, CA. 
9020 VIRTUAL STORAGE REPRR ENG Eo ee a) ore er eet or 


THE FOLLOWING OPTIONS ARE AVAILABLE; 
1. MAKE METAPILE VOPR SP REV ICUSD SSA el eGu ae. 
Jae (O10) ab 


a9 : 


BND Ree az 


@ 

a 

WHAT FIBE NAME IS THE DATA STORED UNDER? 
filtereg 
THE CURVE DATA IS BEING LOADED FROM FILE FILTERES DATA 
2? USING A PRE-ALLOCATED DATASET FOR UNIT FT18F001. 


THE FOLLOWING OPTIONS AKE AVAILABLE: 


leek METAR GGeOn PREVIOUSLY SAVED CURVE. 
eee wee Weer 


ENTER 1 OR 2 


= 

2 

Poemeeume tool iA 930 ——- 13201 VECTORS GENERATED IN 1 PLOT FRAMES 
Pee nie tARY SOFTWARE PRODUCT OF ISSCO, SAN DIEGO, CA. 

Moen TUAL STORAGE REFERENCES; 9 READS; 0 WRITES. 

Daob12\ DETACHED 

DASD 122 DETACHED 

DASD 120 DETACHED 


FEA IK A I A a a eK A i a a eK eK ee oe oe ee te ee KK KK RK KK Kk 


BOmrou WANT TA VRSTEC PLOTTER SNOOTH COPY OF THE 
mes DISSPLA AOE Ny MOUmoUon CheArE yD ? 


EA IKK ee eK i eK AK ee a i a a a Ke eK eK KK eK ok ok 


PA 7/0 

CAS D Ob 1 DPMS ween /Osan/WaBY “MVS; R/O BY O700P 
Z ey RC. OS 

DASD 001 DETACHED 

Creat inG NEW FILE: 

CREATING NEW FILE: 

BUNefite 831/ TO MVS COPY 001 NOHOLD 
DASD 120 DETACHED 


Jo koko ok kfok kok koko doko kok dik kk kok kk kkk kkk doko kk kok kk dk ook ok ok ok 


YOUR GRAPH(S) CAN BE PICKED UP AT THE COMPUTER CENTER. 
TEE GRAPH (S) WILL BE ADDRESSED TO "POP (USER ID)". 
oe eK ok kee ee eK oe KK KK i ee eK KK KK i ee i ee KK KK KK KK KKK Ke 


He eK KK eK ee ee ee eK ee ie oe oe ee EK Ka a Ke OK 


DO YOU WANT TO 


ono NOR Lo. So keaGa LN 
Cemenon tHee FLOR PR OG 


RAM U 
(TO PLOT ANOTHER TYPE 
QUIT 


Sil 
OF 


ENTER, 2 ORS 
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We He HK I i eo Ae ie a i eK oe eK eo i ee i ee ee ee ae a ek oe Xe eo oe KK ok ok aK eo ok OK 
3 


IO kk a kokok dkok lok ok loki ak kkk kk kook kak dai ak de ak ak ak ak ak ak ok ata ak ak ok 


HAVE S25 GCOD eA 


FOI GOK oko kk io a kolo ick loki doi kok kok dak kodak oko kok akok ok ak ak ak ake ake ok ak otek ake a ok 
R; T1959. 0072 2530S = 7 se 
record off 
END RECORDING OF TERMINAL SSSION 
The grarhical output generated by this exanple follows 


as figure 3.4. 
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FIRVER + REGUMATOR CLOSED ECO resis lial 


EXAMPLE 0.1 
LINEAR OPTIMAL CONTROL SYSTEMS 


Seto 


Geto 


0°0¢ 


Oof'O 


o'oT 
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1? - NOILISOd UVTINDNV 
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Ta? —- ALVWILSA NOILISOd YVINDNV 


0°0 O'OI- O0'OZ- O'OE- O'0%=- 0'09- 
A- dOVLTOA LOdNI 


2 


0°09- 0'02- 


0.3 
TIME — SEC 


Figure 3.4 


Filter plus Regulator Closed-loop Time Response 


IV. CONCLUSIONS AND RECOMMENDATIONS 


A. CCONCIUSIONS 


As an ultimate evaluation or the computational abilities 
of OFTSYSX, the program was tested using an 82 X 82 matrix 
of aircraft longitudinal motion equations for the i-29A 
experirental forward-swept wing Fighter aircraft prototyfe, 
Provided by NASA-Edwards. 

For a system of equations of this magnitude, all progran 
arrays were re-dimensioned, and a 2-mejgabyte virtual machine 
size was required. 

The graphical time response curves generated irom the 
X~29A longitudinal system matrix follow as figures @. ive 
4.3 and 4.4. The accuracy of these time response curves is 
mixed. All of the states shown have the correct waveforms, 
but differ ina scale factor of approximately times 10.0. 
Unfortunately the data supplied by NASA was not explicit 
regarding hcw the control input was applied, and whether any 
additional gains were used in their simulation of the 
system. Time constrains did not allow the clarification of 
these prceblem areas, however the results of the X-29A longi- 
tudinal system are encouraging and should be a topic of 


further research. 
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It is hoped that Sencrol System instructors will 
encouraye their students to use this interactive graphical 
time response pregram for all applicable class projects; and 
that its enhanced capabilities will stimulate both interest 
in and research on basic systems control problems, as well 


as more advanced designs. 


Fe. RECCHEMENDATIONS 


Based on the results of this thesis, three areas emerged 


as possibilities for further research and study: 
1. Erogram 


The use of OPTSYSK and Similar design programs 
should be encouraged in all undergraduate and graduate level 
courses involved in tke analysis and design of control 
systems. Toward this end, it is recommended that the OPTSYS 
time response programs be placed in the non-IMSL likrary of 
Subroutines, making it easily available to all pctential 


users. 


2. Frogram Memory Requirements 


When configured for large matrix operations (98X98) 
the CETSYSX program requires over 2 megabytes of virtual 
memory. Virtual machines of. this size are not normally 
available tc a user. The memory usage for matrix storage is 
a possiktle area of improvement in the efficiency of the 
CPTSYSXK program design. All matrices calculated in OPTSYSX 
{except DUMMY matrices) are still available when the run is 
finished. This simplifies program operation but uses an 
excesSive amount of memory. Memory usage should re studied 
and program modifications should be made to reduce the 


excessive memory reguirements. 
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3. PUEther Modi ticarrons 


Program sequencing during optimal filt2r synthesis 
Should be modified. At the present time a regulator must be 
designed or supplied when a filter is designed. Various 
test runs indicate that this problem can be overridden if 
the number cf controls (Nc) 1S given as zero, but this is 
not a viakle soluticn for systems which use a driving 
function. 
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APPENDIX AW OPISYS EXEC LISTING 


SIR aGe OFF 

FO ok doko ok tok kk kkk dk kk dk dk dak dk dkokok dk doa teak dk kok akote cake 
x 

“ PROVOET Sis (2 5G 

x CONTROLS TE RgOPTSS hye ORECALC ANDmOPTPLOT 

é TC DETERMINE THE TIME RESPONSE OF A SYSTEM. 

“ BY H. A. DIEL 

Fs VERSION 1.0 16 JULY 19384 


Me He ee eK Ke ea ee ee he ee ae fe ate he eee ae oe oe ak ee ae ake ole oe ae Coe oe i ee ie ake ok ke fe oe ok oe oe eke teak o& 
We A He He ee ee eee ie ee i ae eo ae 2 2 eo he a 2K fe ee ok oe ae oe oe ake kc 2 ie ke ak oe oe ke ok ek ok 


* * 
a CHECK FOR USER'S VM SIZE = > THAN 1 MEGBYTE : 
FR RO ak lok kk tk kk kkk kok doko kak ak ok de skate keke ook ok ok ok aoe ook ok 
Wiis ZF 

SIF &kC GE 1024 §&GOTC -THO 

CLESCEN 


pemeayPh -~ENDTHREE 
AOR I OR kk kk dik ok kk sk dk kok kok ak ok ok ke ake ok ek 


YOU MUST HAVE A 1M CR LARGER VIRTUAL MACHINE 
Om@RUN THISMOPTSYS PROGRAM 


HOS DEN INE AIM ViIRTUAS MACHINES 


DEFINE STORAGE 1M PKESS ENTER 

i CS PRESS ENTER 

OPTSYS PRtoSs ENTER 

FOR SYSTEMS LARGER THAN 32 X 32 
OETAIN A LISTING OF THE OPTSYS PROGRAM 
AND FOLLOW INSTRUCTIONS CONTAINED IN THE LISTING. 

$k KO GR kk kok doiok kok lok kkk doi iki ok dokkok dokoiiokk kodkok dk kok iok dodo kkk dk kok ok te ok 
-ENDTHREE 
SEXIT &RC 
= IE. 
CUR s Cc 


RN 
cpg. ie 8 ~ENDZERO 


ieee ls Youn “eC "CC NTROLS AQTRLO OF SPROGRAMS: 

WaemeOrTSYS x FORTRAN So. bu Pipeaninne Gee 

2- OPICALC FORTRAN CALCULATE TIME RESPONSE 

3. OFTPLOT FORTRAN Piss oA PLOTTING “.OULINE) 
PmenerRGGRAM PASSES INBORMATION TO,THE NEXT PROGRAM _ 
MoameuGHheA DATA FILE WRITZIEN TO THE USERS DISK. IN THIS 
eer teron Flees ARE "OFTMAT DATA AND “OPTPLOT DATA". 
inmeeores Of THESE FILES VARY WETH THE SYSTEM ORDER, ANT 
Gwe uok ABOUT 20% Cr THE USERS DISK SPACE. THEREFORE 
MyouRE THAT SUFFICIENT DISK SPACE IS AVAILABLE. 


- TYPE "E" TO EXIT, ANY OTHER ENTRY TO CONTINUE - 
-ENDZERO 
GREAD VAES GANS 
SIF .oANS ere wemor 1 eeore 
KK He ee Kk IK IK 5 45G HSE eRET EG eae ee ee eee 
o x 
* ALLOW THE USE OF AN OLD "OFIMAT DATA A1™ 


FRC GG Sok tok io kk ka kkk iia tok dak kk ak dak oe koe ok ok 
RENAME OFTMAT DATA Al OPTSYS DATA Al 

GIf Chen haere COMO metho. 

RENAME OFTSYS DATA Al OPTMAT DATA At 


CLRSCRN 
SBEGIYPE -ENDONE 


HEI ie eK a ae ee ee ee ea ie ie ie ie ie He Ke ie a ie ee Sie eo oo i ik ke ek Kk KK OK 


YOU HAVE A DATA FILE NAMED ‘OPTMAT DATA! ON YOUR A 
DISK THAT WAS PREVIOUSLY GENERABEDESY The (OR tsis 
PROGRAM AND CONTAINS THE F, G, H, GAMMA, A AND B 
MATRICES FROMM GUAT SCUNG 


a 
R 
A 
F YOU WOULD ILKE 70 USE THESE SAC wh fae fo eo 
HIS...RUN, THE OPTSYS Pe OG Fr eee erga ee erdee 
ESIRED DATA AT THE APPROPRIATE rs 


Lf YOU TYRE { jeto- 


ANY OTEER INPUT WILL RESULT IN THAT FILE SEiNe aS 


SOR RR OK ok ok ka kk ka kkk kk oka ok kit oie oe ic ke kat ake ok 2 
~iNDONE 

SREAD VAES GANS 

SIF .S6ANS EC .~.Y &GOTC -ONE 


=F Te Om 
eA HE He HH eK 3 ee ie ee HK ok ee ie a 2k ee eo A aie ee a ie ie ae oe ae ie ee ie a oe ok i a ie ke ok eK AK 


re 


x * 
*. ERASE THE OLD MOPINATODA Tim eet re 7% 
* PLACE "O00 0" TNSTEE NEW “ORT ee x 
i TO ACT AS A FLAG FCR OPTS YSX AND OPTCALC : 
HK KR RRR RR ORR RR a I kk a aka Ka kkk ak ok ke ok ca 
FRASE OPIMAT DATA At 

SSTACK 000 0 

apo aen OPTMAT DATA A Ugh ec Ua 

SH Pn 

CLRSCEN 

SBLGTYP Es ENDER OUR 


eK 3K IK IK eK i eH He ee eK I ea A i oe ie a ee ae i ee ee KK KK KK ee eK A ee ok KK KK 


DC YOU WANT THE NUMERICAL OUTPUT FROM CPTSYSX TO GO 
™C YOUR TERMINAL S (CRE EN) CR TO A D(ISK) FIL#? 
(S OR D 
Me He Ke He i eK ie ee ee oe ee ie ie he oie ie ae ek ee Se ae i oe ie oe ae ke i oe i KK i ek ie eo ok ok oe ok ok ek kok 
-ENDFCUR 


SREAD VAKS €EANS 
SIF .&ANS EQ 2S €GOTC —-eure a 
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oo 


Ee er 


meeeecanes £C .D &GOTC —FIFI A 
Chno GaN 
SBEGTYPz -ENDFIVE 


eK AE Re oe oe ee oe ee oe oe ee oe eo i ae tee oe ake ake oe oo ok ke oe oe 2 ok oe oe oe ok ok ke ok kK ok ok kK KKK KK 
YOU MUST ANSWER S(CREEN) OR D(ISK). 


He Be He RE ee ie eae ee a ae he oe oe ie oe ee oe KK RK OK KK KK KKK KK KKK 


whenor ive 
Emo 


He Fe HE Ek ae eee eee eae ee ee 2 he 2 oe ie A ee ie oe ee eK ke Ok oe ee a KK Ke ORK KKK 


OUTPUT WILL COME TO YOUR TERMINAL SCREEN. 
eR ok oe tok ko kk kkk kok ie dk ok ek ok dk kok kk ok kk ok & 


-ENDSIX 
CP SLEEP 1 SEC 

ETYPE LOADING OPTSYS.... 

FILEDEF 06 TERM (RECEM FA BLKSIZE 133 

FILEDEF 8 DISK OPTPLOT DATA AT (PER 

FILEDEF S DISK OPTMAT DATA Al (PERM 

GLOBAL TXTLIB FORTMOD2 MOD2EEH IMSLDP WNONIMSL 
LOAD OPTSYSX (START 

SGOTC -FIVE 

-FIFTH 

CLRSCEN 

EBEGIYPE -ENDSEVEN 


He He He ee ie Ke oe he oe eee ee eee ee ee ie ie aye ee A oe He 2c ke ae KK oe ee eK ee KK Ae eK KO KK 


Comevur NebL GO Te DISK Tite *OUDEUTX LISTING Al 
He ek ee ke kak kk ke ko ak kk a ook ke ke kek oo oe ko ke a oo kook ke ok ke ok oe ok ek ook ok 


=EniDio BV EN 

Ceo erEP | SEC 

Ceeree LOADING OPTSYS.... 

Weeks OG DISK OUTPUTX LISTING Al 

Pereoee 8 DISK OPIPLCT DATA Al eae 

FILEDEF 9 DISK OPTMAT DATA Al (PERM ; 
GLOBAL TXTLIB FORTMOD2Z MODZEEH IMSLDP NONIMSL 
ae (START 
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EO ENE GH 


HK He He He He eK eK a ae ea ae ee ie ee ee ee A eR KK KE OE RK ok eke ok 


IF YOU ARESDUSSATISI I EDR Tee ceo oS 
THUS FAR AND WOULD LIKE TOSr tere Oser., 


Sol PE we Orn kr 


ANY OTHER INPUT £O CONTiMUE) 
eR ee a a RR a RR oe Ra a oo ao a eo Rk eo 


-ENDEIGHT 

EREAD VARS ans 

CIF = SAN ee Shien eRe 

He He ee A ek ee ee Hee he Hee he a oe oe ok ee ae ee ee ee ee ee oe ok ee ek ee a a KOK EK ok ek 
* 4 
“ CHECK FOR DATA IN THE FILE "OPTMAT DATA " * 
‘ BEFORE LOADING OPTCALC + 


$k I ok tok kk Rk kK RR do kk ok eK 
FIND on ak Stes eae Ai cO1 0 LIMOO2 ALL GiOCUP1 


SREAD VARS &A1 GA 

SIF .&A2 EQ .0 SEXIT ESRC 

STYPE LCADING OPTCALC.... 

FILEDEF 8 DISK OPTPLCT DATA 44 (BER i 

FILEDEF 9 DISK OPTMAT (PERM 

GLOBAL TXTLIB FORTMODZ "OD2 EEE IMSLDP NONIMSL 
LOAD CETCALC (START 

-SEVENTH 

CLRSCEN 

SBEGTYPE -ENDNINE 


RR oR ko ko dk fk ak ok dk kk tok 
If yOU ARE DISSATISFIED Weel oo UTS 
THOS FAR AND WOULD LIKE 20 Exi® TOeene 
~T YEE Tome ell 


(ANY OTHER INPUT TO CONTINUE) 
Me HK IK KK Ke He Ke ee eee a ee ee i he he he oe eK ee ee eK ok eK KK ee KK i kK eek Kok & 


-ENDNINE 
READ VARS §ANS 

SIF .SANS EQ .Y SEXIT &RC 

FILEDEF 8 DISK OPTPLCT DATA A1 (PERM 

EXEC ,DISSPIA OPTELOT 

MK HH KK ee ee ee ee ae ee a ee ee ee ee ee EK ee RE ee Ke KK ok eK KK 
ss x 
* CHECK FOR FILE "DISSPLA METAFILE A4" ON + 
+ THE USER'S [ISK BEFORE GOING TO DISSPOP x 


Re I RR Ra RR RR RR a a a Ra a oo eo oe ok eo ok Ke KK 
RENAME DISSELA METAFILE AY OPTSYS METAFILE A4 

SlF ERC KE 0 €COTO] -—TEiea 

RENAME OFTSYS METAFILE AY DISSPLA METAFILE A4 
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CLRSCRN 
Spare! —~ENDTEN 


Fe He ee Re a ae ee ee a ee a a KK ko oko ok 


PemeoOU WANT A VRSTEGC PLOTIER SMOOTH COPY OF THE 
ie too Pe MOR Shh LOUIS. CREATED 2 


BR He He a ee ee a a tee eae ee ae he ie ie ee ea i he oe ke ok ook ok ok ote eo ok oe ete oko ok 


=ENDIEN 

GSREAL VARS SANS 

LF «SANS EQ .Y &GOTC -NINI A 
Gum «CANS E@e.N EGORO -TENTH 
GLERSCEN 

pooG ©YPE —~ENDELEVEN 


FA FE A I CK I I Ke i A a ke ae ee eC ae oe A oe eK ok eK ok RK KK 


YOU MUST ANSWER Y(ES) OR N(O). 
Me We ee ee he ee ee i ee ae eee ee he ee i oe ee ee ae ok oie oe oe ie ake Ce oe a ee ke ee ek ke ke kk ok OK KK 


~ENDELEVEN 

Ce SLEEP 4 SEC 
&EGOTC -EIGHTH 

enti tH 

x EC et VAS DEC 


CLRSCRN 
GBEGTYEFE -ENDIWEILVE 


AE He ee ee a ee ee ee ee ee ee ee ee ee a a Ke i KK eK Ne eK a Ko 


YOUR GRAPH(S) CAN BE PICKED UP AT THE COMPUTER CENTER. 
PheeonaArH (S) Witt BE ADDRESSED TO "POP (USER ID)". 


Me He ee oe eo he ee ee ee eo ee eke ae ke ie oe a ee ee i ee Ke ae eK ee eK KK 


-~ENDIWELVE 

Crmoe th? 5. SEC 

eS. hu 1 o 

CLRSCEN 

SebGaiyPs ~ENDTHLIRTEEN 


Be a eke eo ee oe a a oe eK KK EK KK KK KK KE 
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DO YOU WANT TO 
Tomek UN OR TSS nahn 
Ze RUN THE PLOT PROGRAM USING Dawes a kee, 
Coit ANOTHER TYPE OF SYSTEM Cran CEOSED 


BN TPES, 2 Ones 
FO ak dR doko ok kkk kkk ak kde ok ok ak kk kkk dk ok oko ak ok ok dk kkk ak ok ek dk ook dott ok 


—~ENDTRae wean 
ne SANS 


aie 


SIF .€ANS PCa cGClOG —THine® 
SIF SANS BE SZ CGOTe =—S iki 
CLRSCRN 

SREGITYPE -~ENDGOODBY 


He Fe HK oe I ee ee eK oie eK ee IK eK ae Fe oie eae ee oe I eK oe ae ie oe ae XS oe ie eo ok i ok ok oe Ke ek ek & 


HAVE A GOOD DAY!! 


eH HK ie Re ee Ke oe oe A i ee ea ee ee ee ie i ie ie OK ee ie ee ee ee ke oe oe ie ke kok 


=~EbN DGGODEY 
CP SEEREVs ore 


CER E Tene 


LTSWISSSITPar sor STOr /I0dd/ NOWWOD 


CUL2E) On Sr hae tae 
TIAH SETS TITZSM) SCCTSTIANS OTS TOTTMD FC CTS TIMDS CTS TOTTIM) SADNAIVAINGSA 


PATI SAVLsSar*’ 3arséMSQsST SES LIS Zair st 
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